RADIATIVE PROCESSES IN HIGH 
ENERGY ASTROPHYSICS 



Gabriele Ghisellini 

INAF - Osservatorio Astronomico di Brera 
February 28, 2012 



2 



Contents 



1 Some Fundamental definitions 

1.1 Luminosity 

1.2 Flux 

1.3 Intensity 

1.4 Emissivity 

1.5 Radiative energy density 

1.6 How to go from L to u 

1.6.1 We are outside the source 

1.6.2 We are inside a uniformly emitting source . . . . 

1.6.3 We are inside a uniformly emitting shell 

1.7 How to go from I to F 

1.7.1 Flux from a thick spherical source 

1.8 Radiative transport: basics 

1.9 Einstein coefficients 

1.9.1 Emissivity, absorption and Einstein coefficients . 

1.10 Mean free path 

1.10.1 Random walk 

1.11 Thermal and non thermal plasmas 

1.11.1 Energy exchange and thermal plasmas 

1.11.2 Non-thermal plasmas 

1.12 Coulomb Collisions 

1.12.1 Cross section 

1.13 The electric field of a moving charge 

1.13.1 Total emitted power: Larmor formula 

1.13.2 Pattern 

2 Bremsstrahlung and black body 

2.1 Bremsstrahlung 

2.1.1 Free-free absorption 

2.1.2 From bremsstrahlung to black body 

2.2 Black body 



Beaming 

3.1 Rulers and clocks 

3.2 Photographs and light curves 

3.2.1 The moving bar 

3.2.2 The moving square 

3.2.3 Rotation, not contraction 

3.2.4 Time 

3.2.5 Aberration 

3.2.6 Intensity 

3.2.7 Luminosity and flux 

3.2.8 Emissivity 

3.2.9 Brightness Temperature 

3.2.10 Moving in an homogeneous radiation field 

3.3 Superluminal motion 

3.4 A question 

Synchrotron emission and absorption 

4.1 Introduction 

4.2 Total losses 

4.2.1 Synchrotron cooling time 

4.3 Spectrum emitted by the single electron 

4.3.1 Basics 

4.3.2 The real stuff 

4.3.3 Limits of validity 

4.3.4 From cyclotron to synchrotron emission . 

4.4 Emission from many electrons 

4.5 Synchrotron absorption: photons 

4.5.1 From thick to thin 

4.6 Synchrotron absorption: electrons 

4.7 Appendix: Useful Formulae 

4.7.1 Emissivity 

4.7.2 Absorption coefficient 

4.7.3 Specific intensity 

4.7.4 Self-absorption frequency 

4.7.5 Synchrotron peak 

Compton scattering 

5.1 Introduction 

5.2 The Thomson cross section 

5.2.1 Why the peanut shape? 

5.3 Direct Compton scattering 

5.4 The Klein-Nishina cross section 

5.4.1 Another limit 

5.4.2 Pause 



CONTENTS 



5 



5.5 Inverse Compton scattering 

5.5.1 Thomson regime 

5.5.2 Typical frequencies 

5.5.3 Cooling time and compactness 

5.5.4 Single particle spectrum 

5.6 Emission from many electrons 

5.6.1 Non monochromatic seed photons 

5.7 Thermal Comptonization 

5.7.1 Average number of scatterings 

5.7.2 Average gain per scattering 

5.7.3 Comptonization spectra: basics 

Synchrotron Self Compton 

6.1 SSC emissivity 

6.2 Diagnostic 

6.3 Why it works 

Pairs 

7.1 Thermal pairs 

7.1.1 First important result 

7.1.2 Second important result 

7.2 Non-thermal pairs 

7.2.1 Optical depth and compactness 

7.2.2 Photon-photon absorption 

7.2.3 An illustrative case: saturated pair production . 

7.2.4 Importance in astrophysics 

Active Galactic Nuclei 

8.1 Introduction 

8.2 The discovery 

8.2.1 Marteen Schmidt discovers the redshift of 3C 273 

8.3 Basic components 

8.4 The supermassive black hole 

8.5 Accretion disk: Luminosities and spectra 

8.5.1 Spectrum 

8.6 Emission lines 

8.6.1 Broad emission lines 

8.6.2 Narrow emission lines 

8.7 The X-ray corona 

8.7.1 The power law: Thermal Comptonization .... 

8.7.2 Compton Reflection 

8.7.3 The Iron Line 

8.8 The torus and the Seyfert unification scheme 

8.8.1 The X-ray background 



6 



CONTENTS 



.9 The jet 144 

8.9.1 Flat and steep radio AGNs 144 

8.9.2 Blazars 146 

8.9.3 The power of the jet 152 

8.9.4 Lobes as energy reservoirs 154 

.10 Open issues Il56 



Chapter 1 



Some Fundamental 
definitions 



1.1 Luminosity 

By luminosity we mean the quantity of energy irradiated per second [erg s~^] . 
The luminosity is not defined per unit of solid angle. The monochromatic 
luminosity L{l') is the luminosity per unit of frequency v (i.e. per Hz). The 
bolometric luminosity is integrated over fequency: 



Often we can define a luminosity integrated in a given energy (or frequency) 
range, as, for instance, the 2-10 keV luminosity. In general we have: 



Examples: 

• Sun Luminosity: Lq = 4 x 10^^ erg 

• Luminosity of a typical galaxy: Lg^i ~ IO^^Lq 

• Luminosity of the human body, assuming that we emit as a black- 
body at a temperature of (273+36) K and that our skin has a surface 
of approximately S = 2 m^: 



This is not what we loose, since we absorb from the ambient a power 
L = SaT^^^ ~ 8.3 x 10^ erg/s if the ambient temperature is 20 C 
(=273+20 K). 




(LI) 




(1-2) 



Lbody = SaT^ ~ 10^° erg/s ~ 10^ W 



(1-3) 
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1.2 Flux 

The flux [erg cm~^ s~^] is the energy passing a surface of 1 cm^ in one 
second. If a body emits a luminosity L and is located at a distance R, the 
flux is 

1.3 Intensity 

The intensity / is the energy per unit time passing through a unit surface 
located perpendicularly to the arrival direction of photons, per unit of solid 
angle. The solid angle appears: [erg cm^^ s^^ sterad"^]. The monochro- 
matic intensity /(i/) has units [erg cm~^ s~"^ Hz~"^ sterad"^]. It always obeys 
the Lorentz transformation: 



Z/3 (j,/)3 



invariant (1-5) 



where primed and unprimed quantities refer to two different reference frames. 
The intensity does not depend upon distance. It is the measure of the irra- 
diated energy along a light ray. 



1.4 Emissivity 

The emissivity j is the quantity of energy emitted by a unit volume, in one 
unit of time, for a unit solid angle 

' = dvdtd^^ ^'-'^ 

If the source is transparent, there is a simple relation between j and /: 

I = jR (optically thin source) (1.7) 



1.5 Radiative energy density 

We can define it as the energy per unit volume produced by a luminous 
source, but we have to specify if it is per unit solid angle or not. For 
simplicity consider the bolometric intensity /. Along the light ray, construct 
the volume dV = cdtdA where dA (i.e. one cm^) is the base of the little 
cylinder of height cdt. The energy contained in this cylinder is 



dE = IcdtdAdn 



(1.8) 
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In that cylinder I find the hght coming from a given direction, I can also say 
that: 

dE = u{n)cdtdAdn (1.9) 

Therefore 

u{n) = ^ (1.10) 

If I want the total u (i.e. summing the light coming from all directions) I 
must integrate over the entire solid angle. 



1.6 How to go from L to u 

There at least 3 possible ways, according if we are outside the source, at a 
large distance with respect to the size of the source, or if we are inside a 
source, that in turn can emit uniformly or only in a shell. 



1.6.1 We £ire outside the source 

Assume to be at a distance D from the source. Consider the shell of surface 
AttD'^ and height cdt. The volume of this shell is dV = AirD'^cdt. In the 
time dt the source has emitted an energy dE = Ldt. This very same energy 
is contained in the spherical shell. Therefore the energy density (integrated 
over the solid angle) is: 

^ n in 



1.6.2 We cire inside a uniformly emitting source 

Assume that the source is optically thin and homogeneous. In this case 
the energy density will depend upon the mean escape time from the source. 
For a given fixed luminosity we have that the longer the escape time, the 
larger the energy density (more photons accumulate inside the source before 
escaping). Therefore: 

u = ^(W) (1.12) 

For a source of size R one can think that the average escape time is (tcsc) ~ 
R/c, but we should consider that this is true for a photon created at the 
center. If the photon is created at the border of the sphere, it will have a 
probability ~ 1/2 to escape immediately, and a probability less than 1/2 to 
pass the entire 2R diameter of the sphere. Therefore the mean escape time 
will be less that R/c. Indeed, for a sphere, it is (tesc) ~ (3/4)-R/c Therefore: 

_ 3L 3R _9 L 

47ri?3 4c 4 AttR'^c ^ ^ 
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Note that this is more (by a 9/4 factor) than what could be estimated by 
Eq. 11.111 setting D = R. We stress that one can use Eq. ll.l3l onlv for a thin 
source. 



Figure 1.1: A source of radius R emits a luminosity L. The way to calculate 
the energy density of the emitted radiation is different if we are outside the 
source (at the distance D) or inside the source. 



1.6.3 We are inside a uniformly emitting shell 

Assume that a spherical shell of radius R emits uniformly a luminosity L In 
this case the radiation energy density returns to be 

^ (1.14) 



47riZ2c 

and it is equal in any location inside the shell. 

1.7 How to go from I to F 

The relation between the intensity / and the total flux F must account 
for the fact that, in general, each element of the surface is seen under a 
different angle 6. See Fig. 11.21 We should then consider the projected area, 
and introduce a cos 9 term in the integral: 

IcosOdn (1.15) 
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Figure 1.2: This figure illustrates the fact that each element dA of the surface 
is seen under a different angle with respect to the normal n. 

1.7.1 Flux from a thick spherical source 

If we see only the surface of a sphere, namely the source is optically thick, 
and if the intensity is constant along the surface we have symmetry along 
the (j) directions and the solid angle dil. = 2tt sva.9d9 . See Fig. 11.31 We get: 



At the surface, R = r, and we have F = irl. 

Very far away, r ^ R, and we have F = nO'^I = /resource- 

1.8 Radiative transport: basics 

Once the radiation is produced in a given location inside the source, we have 
to see how much of this can leave the source and reach the observer. To 
calculate that, we must introduce the absorption coefficient a^, whose 
dimension is [length"^]. It is defined by the following equation, describing 
the decrement of lu when passing trough an infinitesimal path of length ds: 




(1.16) 



Here R = r sin 6c- Therefore we have 




(1.17) 



dL 



a^Fds 



(1.18) 
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Figure 1.3: The observer P sees the spherical source as a disk of angular 
radius 9c- Each element of the sphere has a different projection. 

The absorption coefficient can be thought as the product of a density of 
"absorbers" times the cross section of the absorbing process: 

Ui, = nuy (1-19) 

But inside a source, besides absorption, there is a contribution to Ii, coming 
from the emitters distributed along ds. The increment of I^ is 

dl^ = j^ds (1.20) 

Therefore, combining emission and absorption, we have the basic equation 
of radiative transport: 

^ = -a,I^+j, (1.21) 
ds 

We solve it in some specific cases: 

1. Emission only: 

dl 

^ = = lufi + juds (1.22) 

where S is the total emitting path. 

2. Absorption only: 

^ = -aylu ^ = -ayds (1-23) 

ds ly 

Note that the form of this relation immediately implies that an expo- 
nential is involved: 

h{s) = lAso)e-^^>-'^'^'' (1.24) 

Before passing the layer, the intensity was 1^^. While passing the layer 
of length s, the intensity decreases exponentially. 
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3. Emission plus absorption: In this case it is convenient to introduce 
the optical depth Tj^: 

dTi, = ttyds = noyds (1.25) 
The transport equation then becomes: 

/ 1/ H > 3— = -ly H (1-26) 



a^ds ay dTy 

We can call source function the quantity Sy defined as 

Su = — source function (1-27) 

dy 

Then the formal solution of Eq. 11.261 is 

h{Ty) = Iy,oe-^^ + r e-^^'^-^^'^Sy{T'y)dTl (1.28) 
Jo 

here Ty is the final value of (i.e. when the intensity has travelled 
the entire distance s). If Sy is constant: 

Ury) = Iy,o e-"" + 5, (1 - e-^^) (1.29) 

IfIu,o = 0- 

Iu{Ty) = ^{l-e-^^) (1.30) 
ay 

Now, a trick: multiply and divide the RHS hy s = R (the dimension 
of the source) to obtain: 

In this form it is immediately clear that when the source is optically 
thin (and Tj, ^ 1), we have 1 — e""^" — )• 1 — 1 + r^, and therefore 

Iu{Ty) = jyR, {Ty < 1) (1.32) 

When instead the source becomes optically thick, and Ty ^ 1, then; 

Iu{Ty) = {Ty » 1) (1.33) 

Ty 

Usually, this happens at low frequencies. The above equation explicitly 
shows that the intensity we see from a thick source comes from a layer 
of width R/Ty, i.e. the layer that is optically thin. In other words we 
always collect radiation from a layer of the source, down to the depth 
at which the radiation can escape without being absorbed (riayer = !)• 
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1.9 Einstein coefficients 

The Einstein coefficients concern intrinsic properties of tfie emitting par- 
ticles. As such, they involve more general concepts than the emission or 

the absorption coefficients. They concern the probability that a particle 
spontaneously emits a photon, the probability to absorb a photon, and the 
probability to emit a photon under the influence of another incoming pho- 
ton. The latter concept, at first sight, is weird. In the current world we are 

surrounded by electronic devices using lasers, and therefore wc arc used to 
the concept of stimulated emission. This was however introduced by Ein- 
stein (in 1917) because, without it, he could not recover the black-body 
formula. 

Let us define a system in which the emitting electron can be in one of 
several energy levels. Each level has a statistical weight gi. The simple case 
that comes into mind is an atom. 



3 g3=4 



2 §8=2 



1 gi=l 



Figure 1.4: The shown 3 energy levels have different statistical weights gi. 

Spontaneous emission: A21 — This process corresponds to a jump of 
the electrons from energy level 2 to level 1 by emitting a photon of energy j^, 
corresponding to the jump in energy E2 — E1. The subscripts 1 and 2 indicate 
energy levels 1 and 2, but are not important, we can equivalently talk of a 
jump between the energy levels 5 and 4... and so on. The corresponding 
Einstein coefficient is A21. Therefore 

A21 = transitionprobability for spont. emission [s~^] (1-34) 

Absorption: B12 — This occurs when the electron is in level 1 and absorb a 
photon of energy hv that corresponds to the energy difference between levels 
1 and 2. The probability of one electron to make this transition depends on 
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Figure 1.5: Illustrative sketch for the processes of spontaneous emission, 
absorption and stimulated emission. 

how many photon there are, and therefore to the intensity Ji, of the radiation 
field. The latter is averaged over the solid angle: = (l/47r) / I^dCl. 

Bi2Ju = transitionprobability for absorption per unit time (1.35) 

Stimulated emission: B21 — This occurs when the electron is in level 
2, and the arrival of a photon of energy hv corresponding to the energy 
difference between levels 2 and 1 makes the electron jump to level 1 by 
emitting a photon. The energy of the emitted photon is the same of the 
energy of the incoming one. It can be shown that also the direction and the 
phase of the two photons are the same. We then have a coherent emission, 
at the base of all laser and maser devices. The probability of one electron 
to make the 2 — >• 1 transition depends on how many photons there are, and 
therefore on the radiation intensity: 

B2iJu = transition probability for stimulated emission per unit time (1.36) 

The Einstein coefficients are completely general, and are valid for any radi- 
ation field (i.e. even out of equilibrium). They are valid even when we have 
a free particle, (i.e. not an atom), because we can define even in that case 
the energy state of the particle. 

In particular, we can find the relations between the Einstein coefficients 
in the particular case of equilibrium, that is when the transitions 2 — >• 1 are 
equal to the transitions 1 — )• 2. Assume then that ni and 712 are the number 
density of electrons in the levels 1 and 2. We must have: 

niBu'Ju = n2A2i + n2B2iJu (1-37) 
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Solving for J^: 



T - A21/B21 

- {m/n2){Bu/B2i) - 1 ^'•^''^ 



At equilibrium the ratio ni/n2 must be: 

ni ^ 51 e-^/'^ ^ 9i hu/kT 

where gi and 52 are the statistical weights of the levels 1 and 2. Inserting 
in Eq. 11.401 we have 

J = ^21/^21 , . 

" {9iB^2/92B2i)eh-/kT - 1 ^ • ^ 

We know that at equilibrium Jy must be equal to the black-body intensity 
By. For this to happen it must be that: 

91B12 = 92B21 

2/11/3 

A21 = -^B2i 1.41 

These are the relations among the Einstein coefficients. Although we have 
found them in one particular case, these relations are always valid, also out 
of equilibrium, because they concern intrinsic properties of the system (note 
for instance that they do no depend on temperature). The relations between 
the Einstein coefficients imply that we need to know only one of them to 
find the others. 



1.9.1 Emissivity, absorption and Einstein coefficients 

The emissivity [erg s~^ cm^^ Hz~^ sterad"^] is related to the Einstein coef- 
ficient A21. If n2 is the density of electrons in the state 2, we must have 

jy = ^n2^2i (1.42) 

The reason of the 47r is because n2A2i is the probability to emit a photon 
(for spontaneous emission) in any direction, while the emissivity jy is for 
unit of solid angle. 

Also the stimulated emission produces photons. Why do we not include 
it the calculus for the emissivity? Because, to have stimulated emission, 
we must have the incoming radiation, and it is therefore more convenient to 
think to stimulated emission as negative absorption. For the total absorption, 
consider the quantity a^Iy. This is the amount of energy absorbed for unit 
volume, time, frequency and solid angle. We also have that niBi2lu is the 
number of the 1—7-2 transitions per unit time, volume and frequency. We 
must multiply by hu to have the energy (per second, Hz, cm'^) and divide 
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by 47r to have the same quantity per unit sohd angle. 1^ drops out and we 
have 

hv 

ay = -7- [niBi2 - n2B2i) 
47r 

/ nm \ „ 

= n2 B21 

\ 91 ) 

,2 



' • (1.43) 



^ f nig2 _ 

\n2g1 J 2hv'^ 

We conclude that absorption and emission are intimately linked, and know- 
ing one process is enough to derive the other one, from first principles. If a 
particle emits, it can also absorb. 

1.10 Mean free path 

The mean free path (for a photon) is the average distance H. travelled by a 
photon without interacting. It corresponds to a distance for which Ty = 1: 

Ty = 1 ^ gMv = 1 ^lu = — (1.44) 

na,. 



If a source has radius R and total optical depth > 1 we have: 

1 R R 



1.10.1 Random walk 



(1.45) 



Assume that a photon interacts through scattering with electrons inside a 
source of radius R and optical depth r > 1. It is a good approximation to 
assume that the scattering angle is random, so the direction after a collision 
will be random. The photon will travel "bouncing" from one electron to the 
next, each time traveling for a mean free path. We then ask: i) how many 
scatterings does it have to do before escaping? ii) How much time does it 
takes? 

Fig. 11.61 illustrates the random motion of a photon inside the source. 
Each segment corresponds to a free path, whose average length is i. The 
total displacement after scatterings will be i?, but if we calculate R simply 
summing all the vectors fi, r2j ■■■■ "Hv, we obtain 0. 

To calculate the distance travelled by the photon we calculate the square 
of R: 

{R') = (rl) + (4) + {rD 

+ 2{fi-f2) +2{fi-f3) + .... (1.46) 
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Figure 1.6: Random walk of a photon inside a source. 



The cross products, on average, vanish, because the cos 9 between the vectors 
vanishes, on average. Therefore we remain with the squares only, and each 
term rj is, on average, long as a mean free path. We have N of these terms. 
Therefore 



(R^) = N{ff) = Nt y = VA^£ (1.47) 

If we want to know how many scatterings, on average, a photon undergoes 
before escaping a source of radius i?, we write: 

y/N = J = R(^n = N = t'^ (1.48) 

Remember: this is valid only if r > 1 (if r < 1 then the majority of photon 
escape the source without interacting...). 

The other question is how long does it take to escape. If ti = i/c is the 
average time for one scattering, we have a total time ttot 



U o R _Rt^ _ R 
Ranc c T c 



hot = Nti = T^- = — = r- (1.49) 



So, on average, the photons makes scatterings before escaping a source 
with a scattering optical depth r > 1. Since each mean free path is long 
R/t, the total distance travelled by the photon is tR and of course the total 
time to escape is tR/c. Beware that this is true on average. As we shall see, 
the formation of the high energy spectrum through Comptonization involves 
those photons that make more than scatterings, because in this way they 
can reach higher energies. 



1.11. THERMAL AND NON THERMAL PLASMAS 



19 



1.11 Thermal and non thermal plasmas 

By definition, a thermal plasma is characterized by a Maxwellian distribution 
of particles. Therefore a non thermal plasma is anything else. The non 
relativistic Maxwellian distribution is: 




(1.50) 



In this form, F{v) is normalized, i.e. F{v)dv = 1. This can be seen 
changing variable of integration, from v to x = mv^ /{2kT), and remember- 
ing that /q°° y/xe~^dx = y/TT/2. 

It is worth to stress that the physics is in the exponential term, while 
the A-KV^dv term is simply equal to dvxdvydvz (in three dimensions). This 
then suggests the questions: it is possible to have a "Maxwellian-like" dis- 
tribution in 2 dimensions? And in one dimension? 

Instead of the velocities we may consider the momenta p of the particles. 
If we write 

p = 7^mc (1-51) 

where 7 = 1/(1 — the above relation is valid both for non-relativistic 

and for relativistic velocities. The Maxwellian momenta distribution be- 
comes (setting 6 = mc^): 

where (1/©) is the modified Bessel function of the second kind. 
Note the following: 



• To define a temperature, the distribution of velocities must be isotropic. 

• Written in terms of momenta, the Maxwellian distribution has the 
same form in the relativistic limit. 

• The Maxwellian distribution is very general, it is a result of statistical 
mechanics. But to achieve this distribution, it is necessary that the 
particles exchange energy between themselves. 

• If competing processes are present (i.e. cooling) it is possible that 
one has a Maxwellian distribution only in some interval of veloci- 
ties/momenta (for instance for low velocities). 

• The exponential term contains the physics, the term p^ is 
simply due to dp^dpydpz = A-Kp^dp. 
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1.11.1 Energy exchange and thermal plasmas 

There are two main ways in which particles can exchange energy among 
themselves: 

1. Collisions 

2. Exchange of photons (emission and absorption; scattering) 

Traditionally, one thinks to collisions as the main driver, and to Coulomb 
collisions as the main mechanism. Of course, the more the collisions, the 
faster the energy exchange and the faster the relaxation towards the thermal 
equilibrium. We have 

# of collisions of a single particle , 
oc density n 

time 

total # of collisions n /-, 
oc n (1.53) 

time 

However, the density n is not the only factor. The other factor is the energy 
of the particle: the cross section decreases with the velocity of the particle 
(and then with its energy, or the temperature of the plasma) . This means 
that it is difficult to have a Maxwellian in hot and rarefied plasmas. 

Ask yourselves: what does it happen if I put particles all of the same 
energy in a box with reflecting and elastic walls? Why should the energy of 
a single particle change? 



1.11.2 Non— thermal plasmas 

In rarefied and hot plasmas, the relaxation time required to go to equi- 
librium, and allow for sufficient energy exchange among particles is long, 
compared to the typical timescales of other processes, as acceleration, cool- 
ing and escape. The particle distribution responsible for the radiation we 
see is then shaped by these other, more efficient, processes. 

The queen of the non-thermal particle distribution is a power law: 

N{E) = NqE-p (1.54) 

When all particles are relativistic we can equivalently write 

N{j) = K-f-P (1.55) 

Usually, N is the density [cm~^] of the particles, but sometimes it can in- 
dicate the total number. Furthermore, one can also specify if N{'y) is per 
unit of solid angle (when one has a distribution that is not isotropic), or 
not. One must also specify in what energy range A^(7) is valid, so in general 
A''(7) has a low and an high energy cutoff (7min and 7max)- Within these 



1.12. COULOMB COLLISIONS 
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two limits, a power-law particle distribution has no preferred energy (there 
is no peak or break). 

Often one deals with a broken power-law distribution, defined by two 
power-laws of slopes pi and p2 joining at some 7break- 

The main acceleration process leading to a power-law distribution (but 
not the only one) is shock acceleration. In this process a particle gains energy 
each time it crosses the shock, and there is a probability, each passage, that 
the particle escape. This is sufficient to yield a power-law distribution. 

1.12 Coulomb Collisions 

Fig. 11.71 shows the trajectory of a particle of mass mi and charge qi when 
it passes near a particle of mass m2 and charge q2- For the example in the 
figure the charges are of the same sign, and m2 3> mi. In general, we can 
introduce the reduced mass fj, defined as 



mim2 






mi + 7712 




Figure 1.7: Sketch of the trajectory of the charge qi (of mass mi <^ m2) 
that is deflected by the charge 52 ■ 
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The impact parameter b corresponds to the minimum distance of the 
unperturbed trajectory from the mass m2 ■ The deflection angle is the angle 
between the incoming and outgoing velocity vectors (both measured at large 
distances from the location of m2). 

One also defines the quantity 



tan6'/2 



(1.57) 



where 6o(u) corresponds to the impact parameter for deflection angles of 
90° (namely, when tan^/2 = 1. In this case 6o(^) = b). bo{v) is related to 
the electrostatic force. One makes the approximation that the interaction 
between the two particles is important only when the particle 1 is close to 
the particle 2. One can then say that the interaction time is of the order of 
t ba/v. The acceleration can then be approximated by a ~ v/t ~ /bo. 
On the other hand, the force responsible for this acceleration \s F = (/i 92/^0' 
and the acceleration is then a = (?i(?2/(m^o)- Equating the two expression 
for the acceleration we have: 



<M2 



Therefore 



b = 



boiv) 

bojv) ^ 

tan6'/2 ~ ^uv^ tan(6'/2) 



(1.58) 



(1.59) 



1.12.1 Cross section 

It is easy to understand that particles with a large impact parameter will 
be deflected much less than particles that could pass close to 52- Also, we 
can understand that particles with larger velocities, thus interacting for a 
shorter time, will be deflected less. We can simply write the infinitesimal 
cross section da as 



da = 2Tr b\db\ 

,db , 
= 2Trb—de 
od 



(1.60) 



The differential cross section (i.e. as a function of deflection angle) is 

27r6 



da 

'dn 



2iT sin 



dbde 

rndB 



b 


d 


\ Mv) 1 


2sin(6'/2)cos(6'/2) 


de 


_tan(6'/2)J 



(^1^2) 



4sin4(e/2)/x2^;4 



(1.61) 
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For electron-proton interactions of non-relativistic particles: 
da 1 



dn Amlv^ sin^(0/2) ' 



v<^c (1.62) 



In the relativistic case, accounting also for the magnetic moment of the 
electron: 

da 1 l-/32sin2(0/2) 



dO. 4mlv^ 72 sin^(e/2) 
Note the following: 

1. 



(1.63) 



da _4 



The number of collisions in the unit of time is proportional to (av) oc 
(a faster electron encounters more protons in the unit of time). 
For a Maxwellian, {v) « T^/^ Then, for the non-relativistic case: 



# of collisions 



time 

Collisions are then less frequent in hotter plasma. 

Collisions with small deflections arc largely favorite. This is because 
it is more probable to have distant particles, rather than close ones. 

As a consequence of 2), the field due to the spin (magnetic moment of 
the electron) is negligible with respect to the Coulomb one. Therefore, 
in general, the term (1 — P'^ sin 6/2) is negligible. 

If we have electron-proton collisions, 

da tq 

oc 



dO, rrigU^ m|/3^ 



where tq = /mc(? is the classical radius of the electron. Wc can see 
that for non-relativistic particles the cross section can be much larger 
than the scattering cross section a^ = (87r/3)ro. 



5. In the relativistic case: 

da 

oc 



dO, m?^"^ 

This simply flags the fact the relativistic mass is larger (by a factor 7) : 
the particle has more inertia, and deflects less, yielding a lower cross 
section. 
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1.13 The electric field of a moving charge 



The treatment follows Rybicky & Lightman (p. 80) and Jackson (§14.1). 
The electric and magnetic fields produced by a moving charge is: 



E{f,t) 
B(f,t) 



k'^R'^ 
n X E 



+ 



ck^R 



n X 



(n - /?) X /3 



(1.65) 



P 




Figure 1.8: A charge is moving along a trajectory. At the time t we measure 
the electric field produced by the charge. To find it out, we have to calculate 
its position, velocity and acceleration at the retarded time tret- 

All the quantities must be evaluated at the retarded time (it is called 
retarded, but is really a earlier time) : at the retarded time the particle was 
at a distance R such that tret + R/c is the time t (the time t is the time of 
the observer that measures the electric field). 

We have the following definitions: 

_ Rjtret) 
'■ret — ^ 

C 

k = 1-n-p 

/? = 

c 

7 = (l-/32)-i/2 (1.66) 
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We see that Eq. 11.651 for the electric field is made by two pieces. 

1. The first is called the velocity field: acceleration does not appear, and 
E is proportional to R^^. 

2. The second is called acceleration field: it contains /3, i.e. the accelera- 
tion, and is proportional to R^^. It is perpendicular to ft and to the 

acceleration /3. 




Figure 1.9: The electric field produced by a charge initially in uniform 
rectilinear motion that is suddenly stopped. At large distances, the electric 
field points to where the charge would be if it had not been stopped. At 
closer distances, the electric field had time to "adjust" and points to where 
the charge is. There is then a region of space where the electric field has to 
change direction. This corresponds to radiation. The width of this region 
is cAt, where At is the time needed to stop the particle, namely, the At 
during which the charge was decelerated. 

For a charge moving at a constant speed (3c, we have only the velocity field. 
There is something apparently weird associated to this field: it points at 
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the location of the particle at time t, not at the retarded time tret- This 
is due to the (n — /3) term. If the particle does not move, one recovers the 
usual E = qn/R'^ Coulomb law. But if does move, E does not point in the 
n direction, even if we have to calculate the position of the particle at the 
retarded time. The (n — /3) term makes the E vector to point somewhat 
along the /3 direction. 



Fig. 11.91 illustrates in a graphic way what the radiation is. Suppose that a 
charge, in uniform motion along a rectilinear path, is suddenly stopped. Very 
far way, the observer "does not know" that the particle was stopped, and 
measure an electric field assuming that the particle has instead continued its 
motion. The electric vector then points to a point that would be occupied 
by the particle, if it had not been stopped (dashed lines in Fig. II. 9p . In a 
region closer to the particle, instead, the information: "the particle is now at 
rest" has been propagated, and the E points radially to the actual position 
of the particle. In between these two regions, E must change direction 

(and value). The change is more pronounced perpendicularly to /3, and is 

negligible along /3 (that here coincides with the direction of the velocity /?). 
The region where E changes direction propagates in time (it goes at c), but 
its width is always cAt. 



1.13.1 Total emitted power: Larmor formula 

For simplicity, we specialize to the non relativistic case, consider only the 
radiative field, and set: 

P < 1 

k = 1 -n-/3 ^ 
n-/3 n (1.67) 

With these approximations we have: 



E(r,t) = — n X (n X P) 
cR 



tret 



B{f,t) 



fix E 



(1.68) 



To calculate the power per unit solid angle carried by this electromagnetic 
field let consider the Poynting vector S (in cgs units): 



S = — E X B = —\E\^n = r fix (fix B) fi 

47r 47r Att c^R^ L lu^t 



The power crossing a surface dA = R?dQ. is 



dP 



dP = SdA = SR'dn ^ — = SR 

ail 



(1.69) 



^1.70) 
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Therefore 



dP 



2 -.2 

n X (n X /3) 



Q 

Attc 

n2 ■ 2 

- — p sm I 
47rc 



2 • 2 

a sm I 



tret 



(1.71) 



This is the Larmor formula for non relativistic particles. The angle 6 is 
the angle between n and the acceleration. The power is null along the 
acceleration, it is maximum perpendicularly to it. Integrating over the entire 
solid angle we have the total power: 



P 



dP 



27r 



sin^ 6 d(cos 6) 



3^" 



(1.72) 



1.13.2 Pattern 





Figure 1.10: The pattern of the radiation emitted by a charge with its 
acceleration parallel to the velocity (top) and acceleration perpendicular to 
the velocity (bottom). 

Eq. 1 1.71 1 shows that the pattern of the emitted radiation has a maximum 
perpendicular to the acceleration, and vanishes in the directions parallel to 
it. Fig. 14.21 shows two examples: in the top part the particle has /3 || /3 while 
in the bottom part /3 _L /3. The particle is non relativistic. 
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Just for exercise, consider an antenna consisting of a linear piece of metal. 
It is called a dipole antenna. What will be the pattern of the emitted radi- 
ation? See: http://e11.wikipedia.0rg/wiki/File:Felder_u111_Dipol.jpg 

What will be the pattern of an oscillating electron? 

We will soon see the modification occurring when the particle becomes rela- 
tivistic, both for the total emitted power and for the pattern of the emitted 
radiation. 



Chapter 2 

Bremsstrahlung and black 
body 



2.1 Bremsstrahlung 

We will follow an approximate derivation. For a more complete treatment 
see Rybicki & Lightman (1979) and Blumental & Gould (1970). We will 
consider an electron— proton plasma. 
Definitions: 



• b: impact parameter 

• v: velocity of the electron 

• rie number density of the electrons 

• Up-, number density of the protons 

• T: temperature of the plasma: mv^ ^ kT ^ v ^ {kT/mY/"^. 

We here calculate the total power and also the spectrum of bremsstrahlung 
radiation. We divide the procedure into a few steps: 

1. We consider the interaction between the electron and the proton only 
when the electron passes close to the proton. The characteristic time 
T is 

b , ^ 

T K - 2.1) 

V 

2. During the interaction we assume that the acceleration is constant and 
equal to 
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3. Prom the Larmor formula we get 

2eV e2 _ 
3c3 ~ c3 m2c364 - ,^2^354 (2.3) 

Note that we have dropped the 2/3 factor, since in this simphfied 
treatment we neglect all the numerical factors of order unity. Later 
we will give the exact result. 

4. Since there is a characteristic time, there is also a characteristic fre- 
quency, namely t~^: 

uj ^ - = - (2.4) 

r b 



5. Therefore 



^(^) ^ - = ^4t3 (2-5) 



6. We can estimate the impact factor b from the density of protons: 

b ^ n-^/^ ^ b^ = — (2.6) 



7. The emissivity j(i^) will be the power emitted by a single electron mul- 
tiplied by the number density of electrons. If the emission is isotropic 
we have also to divide by 47r, since the emissivity is for unit solid angle: 



■'max- 



We integrate j{uj) over frequency. The integral will depend upon 
What should we use for Wmax? One possibility is to set ^max = kT. 
This would mean that an electron cannot emit a photon of energy 
larger than the typical energy of the electron. Seems reasonable, but 
we are forgetting all the electrons (and the frequencies) that have 
energies larger than kT. In this way: 



f^'^"^ , , , neUp /mc\ 1/2 kT 



(2.8) 



We suspect that in the exact results there will be the contribution 
of electrons with energy larger than kT: since they belong to the 
exponential part of the Maxwellian, we suspect that in the exact result 
there will be an exponential... 
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9. The exact result, considering also that z/ = a;/(2-7r), is 



3 J mlc^ \kT J 

The Gaunt factor depends on the minimum impact factor which in 
turn determines the maximum frequency. Details are comphcated, but see 
Rybicki & Lightman (p. 158-161) for a more detailed discussion. 

We have treated the case of an electron-proton plasma. In the more 
general case, the plasma will be composed by nuclei with atomic number 
Z and number density n. The emissivity will then be proportional to Z^. 
This is because the acceleration of the electron will be a = Z / {nicb'^) (see 
point 2), and we have to square the acceleration to get the power from the 
Larmor formula. In cgs units we have: 

= 5.4 X 10-39 zVniT-i/V'^'^/'^^^ 
j = 1.13 X 10-2^ Z^neUiT^/^^ (2.10) 

2.1.1 Free— free absorption 

If the underlying particle distribution is a Maxwellian, we can use the Kir- 
choff law to find out the absorption coefficient. If is the intensity of black 
body emission, we must have 

- -^'^ - R _2hu^ 

In these cases it is very simple to find once we know ji,. Remember: this 
can be done only if we have a Maxwellian. If the particle distribution is 
non-thermal, we cannot use the Kirchoff law and we have to go back to a 
more fundamental level, namely to the Einstein coefficients. Using Eg. 12.111 
we have: 

= ir = o ^ u^2.2 -T^ 9s (2.12) 



Bi, 3 \ 3 y hm'^c^ \ kT J v'^ 
In cgs units [cm~^] we have 

al = 3.7 X 10^ -3 gs (2.13) 

When hu ^ kT (Raleigh- Jeans regime) this simplifies to 
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Log u 

Figure 2.1: The bremsstrahlung intensity from a source of radius R = 10^^ 
cm, density He = np = W^^ cm~^ and varying temperature. The Gaunt 
factor is set to unity for simplicity. At smaller temperatures the thin part of 
Ijy is larger (oc T^^/^), even if the frequency integrated I is smaller (oc T^/^). 

Fig. 12.11 shows the bremsstrahlung intensity from a source of radius R = 
10^^ cm and Ue = np = 10^^ cm~^. The three spectra correspond to different 
temperatures. Note that for smaller temperatures the thin part of is larger 
{Ii, oc r~^/^). On the other hand, at larger T the spectrum extends to larger 
frequencies, making the frequency integrated intensity to be larger for larger 
T (/ oc r^/^). Note also the self-absorbed part, whose slope is proportional 
to u'^. This part ends when the optical depth t = a^R ^ 1. 

2.1.2 Prom bremsstrahlung to black body 

As any other radiation process, the bremsstrahlung emission has a self- 
absorbed part, clearly visible in Fig. 12.11 This corresponds to optical depths 
^ 1. The term in the absorption coefficient ai, ensures that the ab- 
sorption takes place preferentially at low frequencies. By increasing the den- 
sity of the emitting (and absorbing) particles, the spectrum is self-absorbed 
up to larger and larger frequencies. When all the spectrum is self absorbed 
(i.e. Ti, > 1 for all v), and the particles belong to a Maxwellian, then we 
have a black-body. This is illustrated in Fig. 12. 2t all spectra are calculated 
for the same source size (R = 10^^ cm), same temperature {T = 10'' K), 
and what varies is the density of electrons and protons (by a factor 10) from 




10 12 14 16 18 20 

Log u 



Figure 2.2: The bremsstrahlung intensity from a source of radius R = 10 
cm, temperature T = 10^ K. The Gaunt factor is set to unity for simphcity. 
The density Tig = Up varies from 10^'' cm~^ (bottom curve) to 10^® cm~^ 
(top curve), increasing by a factor 10 for each curve. Note the self-absorbed 
part (oc i^^), the flat and the exponential parts. As the density increases, the 
optical depth also increases, and the spectrum approaches the black-body 
one. 

He = Up = 10^" cm~'^ to 10^® cm~'^. As can be seen, the bremsstrahlung 
intensity becomes more and more self-absorbed as the density increases, 
until it becomes a black-body. At this point increasing the density does 
not increase the intensity any longer. This is because we receive radiation 
from a layer of unity optical depth. The width of this layer decreases as we 
increase the densities, but the emissivity increases, so that 

= 1^ oc '^"'^P^ constant (r^ > 1) (2.15) 
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2.2 Black body 

A black body occurs when "the body is black" : it is the perfect absorber. 
But this means that it is also the "perfect" emitter, since absorption and 
emission are linked. The black body intensity is given by 

= I ^f,^^ (2.16) 
Expressed in terms of the wavelength A this is equivalent to: 

Note the following: 

• The black body intensity has a peak. The value of it is different if we 
ask for the peak of or the peak of vBi,. 

The first is at /ifpcak = 2.82 kT. 
The second is at hvpe^k = 3-93 /cT. 

• If T2 > Ti, then: ^^(Ta) > 5^(Ti) for aU frequencies. 

• When hv <C kT we can expand the exponential term: e^'^/^'^ 1 _)- 
hv/kT..., and therefor we obtain the Raleigh— Jeans law: 

/RJ = — kT (2.18) 

• When hu S> kT we have e^"/'''^ - 1 ^ e^"'/'''^ and we obtain the Wien 
law: 

jW ^ ^^-hu/kT (2.19) 



• The integral over frequencies is: 

fB... = ^r^. = (2.0, 

The constant (Jmb is called Maxwell-Boltzmann constant. 

• The energy density u of black body radiation is 

u=^rB^du = aT\ a = ^ (2.21) 
c Jo c 

The two constants {gmb and a) have the values: 

o"MB = 5.67x10"^ ergcm"^deg~"^s"^ 

a = 7.65x10-^^ ergcm-^deg-'' (2.22) 




Figure 2.3: The black body intensity compared with the Raleigh- Jeans 
and the Wien law. 



• The brightness temperature is defined using the Raleigh-Jeans law, 
since I^'^ = {2^"^ /(?)kT we have 

^2 rRJ 

Tb = ^ (2.23) 

• A black body is the most efficient radiator, for thermal plasmas and 
incoherent radiation (we can have coherent processes that are even 
more efficient). For a given surface and temperature, it is not possible 
to overtake the luminosity of the black body, at any frequency, for any 
emission process. 

• Let us try to find the temperature of the surface of the Sun. We 
know its radius (700,000 km) and luminosity {Lq = 4 x 10^^ erg s~^). 
Therefore, from 

Lq = ttAttR^ / B^du = A-KR^aMBT'^ (2.24) 
Jo 



10 

we get: 
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A spherical source emits black body radiation. We know its distance, 
but not its radius. Find it. Suppose we do not know its distance. Can 
we predict its angular size? And, if we then observe it, can we then 
get the distance? 



Chapter 3 

Beaming 



3.1 Rulers and clocks 

Special relativity taught us two basic notions: comparing dimensions and 
flow of times in two different reference frames, we find out that they differ. If 
we measure a ruler at rest, and then measure the same ruler when is moving, 
we find that, when moving, the ruler is shorter. If we syncronize two clocks 
at rest, and then let one move, we see that the moving clock is delaying. 
Let us see how this can be derived by using the Lorentz transformations, 
connecting the two reference frames K (that sees the ruler and the clock 
moving) and K' (that sees the ruler and the clock at rest). For semplicity, 
but without loss of generality, consider a a motion along the x axis, with 
velocity v = fic corresponding to the Lorentz factor V. Primed quantities 
are measured in K' . We have: 

V{x - vt) 

y 

r{t-fi^) (3.1) 

with the inverse relations given by 

X = T{x' + vt') 

y = y' 

z = z' 

t = T(^t' + /3^y (3.2) 

The length of a moving ruler has to be measured through the position of its 
extremes at the same time t. Therefore, as At = 0, we have 

X2 — x'l = T{x2 — xi) — FvAt = T{x2 — xi) (3.3) 
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I.e. 



Ax 



Ax' 



(3.4) 



r 



contraction 



Similarly, in order to determine a time interval a (lab) clock has to be 
compared with one in the comoving frame, which has, in this frame, the 
same position x' . Then 



An easy way to remember the transformations is to think to mesons pro- 
duced in collisions of cosmic rays in the high atmosphere, which can be 
detected even if their lifetime (in the comoving frame) is much shorter than 
the time needed to reach the earth's surface. For us, on ground, relativistic 
mesons live longer (for the meson's point of view, instead, it is the length of 
the travelled distance which is shorter). 

All this is correct if we measure lengths by comparing rulers (at the same 
time in K) and by comparing clocks (at rest in K') - the meson lifetime is 
a clock. In other words, if we do not use photons for the measurement 
process. 

3.2 Photographs and hght curves 

If we have an extended moving object and if the information (about position 
and time) are carried by photons, we must take into account their (different) 
travel paths. When we take a picture, we detect photons arriving at the same 
time to our camera: if the moving body which emitted them is extended, 
we must consider that these photons have been emitted at different times, 
when the moving object occupied different locations in space. This may 
seem quite obvious. And it is. Nevertheless these facts were pointed out in 
1959 (Terrel 1959; Penrose 1959), more than 50 years after the publication 
of the theory of special relativity. 

3.2.1 The moving bar 

Let us consider a moving bar, of proper dimension i' , moving in the direction 
of its length at velocity /3c and at an angle 6 with respect to the line of sight 
(see Fig. 13. ip . The length of the bar in the frame K (according to relativity 
"without photons") is £ = i'/T. The photon emitted in Ai reaches the 
point H in the time interval Ate- After Ate the extreme Bi has reached 
the position B2, and by this time, photons emitted by the other extreme 
of the bar can reach the observer simultaneously with the photons emitted 
by Ai, since the travel paths are equal. The length B1B2 = PcAte, while 
AiH = cAte- Therefore 




(3.5) 



c 




^'cos^ 



(3.6) 



3.2. PHOTOGRAPHS AND LIGHT CURVES 



39 



A, 




H 

to observer 



Figure 3.1: A bar moving with velocity /3c in the direction of its length. The 
path of the photons emitted by the extreme A is longer than the path of 
photons emitted by B. When we make a picture (or a map) of the bar, we 
collect photons reaching the detector simultaneously. Therefore the photons 
from A have to be emitted before those from B, when the bar occupied 
another position. 

Note the appearance of the term 5 = l/[r(l — /3cos0)] in the transformation: 
this accounts for both the relativistic length contraction (l/F), and the 
Doppler effect [1/(1 - /3 cos (9)] (see below, Eq. [XT5|) . The length A1B2 is 
then given by 

^1^2 = ^ = ' = Si'. (3.7) 

cos 6^ r(l — pcosd) 

In a real picture, we would see the projection of A1B2, i.e.: 

sm 

HB2 = ^1^2 sine = £'— = i' 6 sine, (3.8) 

r(i — pcosO) 

The observed length depends on the viewing angle, and reaches the maxi- 
mum (equal to i') for cos 9 = j3. 

3.2.2 The moving square 

Now consider a square of size in the comoving frame, moving at 90° to the 
line of sight (Fig. 13. 2p . Photons emitted in A, B, C and D have to arrive 
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Figure 3.2: Left: A square moving with velocity /3c seen at 90°. The observer 
can see the left side (segment CA). Light rays are assumed to be parallel, 
i.e. the square is assumed to be at large distance from the observer. Right: 
The moving square is seen as rotated by an angle a given by cos a = f3. 



to the film plate at the same time. But the paths of photons from C and 
D are longer — )• they have to be emitted earlier than photons from A and 
B: when photons from C and D were emitted, the square was in another 
position. 

The interval of time between emission from C and from A is i' /c. During 
this time the square moves by f3i', i.e. the length CA. Photons from A and 
B are emitted and received at the same time and therefore AB = i' /T. The 
total observed length is given by 

i' 

CB = CA + AB = -{1 + r/3). (3.9) 

As /3 increases, the observer sees the side AB increasingly shortened by 
the Lorentz contraction, but at the same time the length of the side CA in- 
creases. The maximum total length is observed for /3 = l/\/2, corresponding 
to r = and to CB = I' i.e. equal to the diagonal of the square. Note 
that we have considered the square (and the bar in the previous section) to 
be at large distances from the observer, so that the emitted light rays are 
all parallel. If the object is near to the observer, we must take into account 
that different points of one side of the square (e.g. the side AB in Fig. 
13. 2p have different travel paths to reach the observer, producing additional 
distortions. See the book by Mook and Vargish (1991) for some interesting 
illustrations. 
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Figure 3.3: An observer that sees the object at rest at a viewing angle given 
by sin a' = 6 sin a, will take the same picture as the observer that sees the 
object moving and making an angle a with his/her line of sight. Note that 
sin a' = sin(7r — a'). 

3.2.3 Rotation, not contraction 

The net result (taking into account both the length contraction and the 
different paths) is an apparent rotation of the square, as shown in Fig. 13.21 
(right panel). The rotation angle a can be simply derived (even geometri- 
cally) and is given by 

cos a = /3 (3.10) 

A few considerations follow: 

• If you rotate a sphere you still get a sphere: you do not observe a 
contracted sphere. 

• The total length of the projected square, appearing on the film, is 

+ l/F). It is maximum when the "rotation angle" a = 45° — >• /3 = 
l/\/2 — )• r = \/2. This corresponds to the diagonal. 

• The appearance of the square is the same as what seen in a comoving 
frame for a line of sight making an angle a' with respect to the velocity 
vector, where a' is the aberrated angle given by 

• ' ^^^^ X ■ 

sma = — 7 = osma 3.11 

r(l-/3cosa) ^ ^ 

See Fig. 13.31 for a schematic illustration. 
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Figure 3.4: Difference between the proper time and the photons arrival 
time. A lamp, moving with a velocity f3c, emits photons for a time inter- 
val Atg in its frame K'. The corresponding time interval measured by an 
observed at an angle 6, who receives the photons produced by the lamp is 
Aia = At'J6. 

The last point is particularly important, because it introduces a great sim- 
plification in calculating not only the appearance of bodies with a complex 
shape but also the light curves of varying objects. 

3.2.4 Time 

Consider a lamp moving with velocity ?; = /3c at an angle 9 from the line 
of sight. In K' , the lamp remains on for a time At'^. According to special 
relativity ("without photons") the measured time in frame K should be 
Ate = TAtg (time dilation). However, if we use photons to measure the time 
interval, we once again must consider that the first and the last photons have 
been emitted in different location, and their travel path lengths are different. 
To find out Ata, the time interval between the arrival of the first and last 
photon, consider Fig. 13.41 The first photon is emitted in A, the last in B. 
If these points are measured in frame K, then the path AB is 

AB = PcAte = rpcAt'^ (3.12) 

While the lamp moved from Aio B, the photon emitted when the lamp was 
in A has travelled a distance AC = cAtg, and is now in point D. Along the 
direction of the line of sight, the first and the last photons (the ones emitted 
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in A and in B) are separated by CD. The corresponding time interval, 
CD/c, is the interval of time Ata between the arrival of the first and the 
last photon: 

CD AD -AC 

Aia = = = Ate - /3Ate COS 6 

c c 
= Ate(l -/3cos6i) 

= At'j:{l- PcosO) 
At; 
6 



(3.13) 



If 9 is small and the velocity is relativistic, then 5 > 1, and Ata < Atg, i-e. 
we measure a time contraction instead of time dilation. Note also that we 
recover the usual time dilation (i.e. Ata = TAt^) if 9 = 90°, because in this 
case all photons have to travel the same distance to reach us. 
Since a frequency is the inverse of time, it will transform as 

1^ = u'5 (3.14) 

It is because of this that the factor 6 is called the relativistic Doppler factor. 
Its definition is then 

^ = T{l-lcos9) ^^-''^ 

Note the two terms: 

• The term l/T: this corresponds to the usual special relativity term. 

• The term 1/(1 — /? cos 9): this corresponds to the usual Doppler effect. 

The 6 factor is the result of the competition of these two terms: for 9 = 90° 
the usual Doppler term is unity, and only "special relativity" remains: 6 = 
l/T. For small 9 the term 1/(1 — (3 cos 9) becomes very large, more than 
compensating for the l/T factor. For cos 9 = (3 (i.e. sm.9 = l/T) we have 
5 = T. Tot 9 = 0° we have 5 = r(l + /3). 



3.2.5 Aberration 

Another very important effect happening when a source is moving is the 
aberration of light. It is rather simple to understand, if one looks at Fig. 13.51 
A source of photons is located perpendilarly to the right wall of a lift. If the 
lift is not moving, and there is a hole in its right wall, then the ligth ray enters 
in A and ends its travel in 5. If the lift is not moving, A and B are at the 
same heigth. If the lift is moving with a constant velocity v to the top, when 
the photon smashes the left wall it has a different location, and the point B 
will have, for a comoving observer, a smaller height than A. The light ray 
path now appears oblique, tilted. Of course, the greater v, the more tilted 
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the light ray path appears. This immediately stimulate the question: what 
happens if the lift, instead to move with a constant velocity, is accelerating? 
With this example one can easily convince him/herself that the "trajectory" 
of the photon would appear curved. Since, by the equivalence principle, the 
accelerating lift cannot tell if there is an engine pulling him up or if there is 
a planet underneath it, we can then say that gravity bends the light rays, 
and make the space curved. 




Figure 3.5: The relativistic lift, to explain relativistic aberration of light. 
Assume first a non-moving lift, with a hole on the right wall. A light ray, 
coming perpendicularly to the right wall, enter through the wall in A and 
ends its travel in B. If the lift is moving with a constant velocity v to the 
top, its position is changed when the photon arrives to the left wall. For the 
comoving observer, therefore, it appears that the light path is tilted, since 
the point B where the photon smashes into the left wall is below the point 
A. What happens if the lift, instead to move with a constant velocity, is 
accelerating? 

This helps to understand why angles, between two inertial frames, change. 
Calling 6 the angle between the direction of the emitted photon and the 
source velocity vector, we have: 

sin 6*' , sin^ 

smt^ = — — — ; smf 



r(l + /3cos0')' r(l-/3cos( 

cose' + f3 , cosO- 13 . 

cosa = — ; cos 6^ = (3.16) 

1 + /3COS0' 1-/3COS0 ^ ' 

Note that, if 9' = 90°, then sin 6' = l/F and cos 9 = /3. Consider a source 
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emitting isotropically in K'. Half of its photons are emitted in one emi- 
sphere, namely, with 9' < 90°. Then, in K, the same source will appear to 
emit half of its photons into a cone of semiaperture l/F. 

Assuming symmetry around the angle (j), the transformation of the solid 
angle is 

dn = 2^dcos9= r2(i + ^eosgO^ = - ^cos^)^ = ^ (3.17) 

3.2.6 Intensity 

We now have all the ingredients necessary to calculate the transformation 
of the specific (i.e. monochromatic) and bolometric intensity. The specific 
intensity has the unit of energy per unit surface, time, frequency and solid 
angle. In cgs, the units are [erg cm~^ Hz"-*^ ster~^]. We can then write 
the specific intensity as 

l(y) = hu 



= Shu' 



dt dv d^ dA 

dN' 



(dt'/S) 5dv' {dn'/S^)dA' 
= 5^l'{v') = 5^l'{v/5) (3.18) 

Note that dN = dN' because it is a number, and that dA = dA' . If we inte- 
grate over frequencies we obtain the bolometric intensity which transforms 
as 

I = S^I' (3.19) 

The fourth power of 5 can be understood in a simple way: one power comes 
from the transformation of the frequencies, one for the time, and two for 
the solid angle. They all add up. This transformation is at the base of our 
understanding of relativistic sources, namely radio-loud AGNs, gamma-ray 
bursts and galactic superluminal sources. 



3.2.7 Luminosity and flux 

The transformation of fluxes and luminosities from the comoving to the 
observer frames is not trivial. The most used formula is L = S'^L', but this 
assumes that we are dealing with a single, spherical blob. It can be simply 
derived by noting that L = Andf^F, where F is the observed flux, and by 
considering that the flux, for a distance source, is F oc IdQ,. Since is 
the source solid angle, which is the same in the two K and K' frames, we 
have that F transforms like /, and so does L. But the emission from jets 
may come not only by a single spherical blob, but by, for instance, many 
blobs, or even by a continuous distribution of emitting particles flowing in 
the jet. If we assume that the walls of the jet are fixed, then the concept of 



46 



CHAPTERS. BEAMING 



"comoving" frame is somewhat misleading, because if we are comoving with 
the flowing plasma, then we see the walls of the jet which are moving. 

A further complication exists if the velocity is not uni-directional, but 
radial, like in gamma-ray bursts. In this case, assume that the plasma is 
contained in a conical narrow shell (width smaller than the distance of the 
shell from the apex of the cone). The observer which is moving together 
with a portion of the plasma, (the nearest case of a "comoving observer") 
will see the plasma close to her going away from her, and more so for more 
distant portions of the plasma. Indeed, there could be a limiting distance 
beyond which the two portions of the shells are causally disconnected. 

Useful references are Lind & Blandford (1985) and Sikora et al. (1997). 

3.2.8 Emissivity 

The (frequency integrated) emissivity j is the energy emitted per unit time, 
solid angle and volume. We generally have that the intensity, for an optically 
thin source, is I = J/^jijdr, where AR is the length of the region containing 
the emitting particles. The emissivity transforms like j = namely with 
one power of 5 less than the intensity. 




Figure 3.6: Due to aberration of light, the travel path of the a light ray is 
different in the two frames K and K' 

To understand why, consider a slab with plasma flowing with a velocity 
parallel to the walls of the slab, as in Fig. 13.61 The observer in K will 
measure a certain Ai? which depends on her viewing angle. In K' the same 
path has a different length, because of the aberration of light. The height 
of the slab h' = h, since it is perpendicular to the velocity. The light ray 
travels a distance AR = h/ sinO in K, and the same light ray travels a 
distance AR' = h'/sind' in K'. Since sin 6*' = SsinB, then AR' = AR/5. 
Therefore the column of plasma contributing to the emission, for 5 > 1, 
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is less than what the observer in K would guess by measuring AR. For 
semplicity, assume that the plasma is homogenous, allowing to simply write 
/ = jAR. In this case: 

/ = jAR = 6^1' = 6^j'AR' ^ j = 6^j' (3.20) 

And the corresponding transformation for the specific emissivity is j(z^) = 



10 9 8 7 6 



OOOOOOO^OO 



AR 




Figure 3.7: Due to the difi^erences in light travel time, the number of 
blobs that can be observed simultaneously at any given time depends on the 
viewing angle and the velocity of the blobs. In the top panel the viewing 
angle is 9 = 90° and all the blobs contained within a certain distance R can 
be seen. For smaller viewing angles, less blobs are seen. This is because the 
photons emitted by the rear blobs have more distance to travel, and therefore 
they have to be emitted before the photons produced by the front blob. 
Decreasing the viewing angle 9 we see less blobs (3 for the case illustrated 
in the bottom panel). 

Fig. 13.71 illustrates another interesting example, taken from the work 
of Sikora et al. (1997). Consider that within a distance R from the apex 
of a jet {R measured in K), at any given time there are N blobs (10 on 
the specific example of Fig. 13. 7p . moving with a velocity v = /3c along the 
jet. To fix the ideas, let assume that beyond R they switch off. If the 
viewing angle is 9 = 90°, the photons emitted by each blob travel the same 
distance to reach the observer, who will see all the 10 blobs. But if < 90°, 
the photons produced by the rear blobs must travel for a longer distance in 
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order to reach the observer, and therefore they have to be emitted before the 
photons produced by the front blob. The observer will then see less blobs. 
To be more quantitative, consider a viewing angle 6 < 90°. Photons emitted 
by blob numer 3 to reach blobs number 1 when it produces its last photon 
(before to switch off) were emitted when the blobs itself was just born (it was 
crossing point A). They travelled a distance RcosO in a time At. During 
the same time, the blob number 3 travelled a distance AR = c/SAt in the 
forward direction. The fraction / of blobs that can be seen is then 

R-AR cPAt 
f = ^ — = 1 ^ = 1 - /3 cos 6' (3.21) 

Where we have used the fact that At = {R/c)cos6. This is the usual 
Doppler factor. We may multiply and divide by F to obtain 

/ = ^ (3-22) 

The bottom line is the following: even if the flux from a single blob is boosted 
by (5^, if the jet is made by many (N) equal blobs, the total flux is not just 
boosted by NS'^ times the intrinsic flux of a blob, because the observer will 
see less blobs if 6* < 90°. 

3.2.9 Brightness Temperature 

The brightness temperature is a quantity used especially in radio astronomy, 
and it is defined by 

where we have assumed that the solid angle subtended by the source is 
Af^s ~ vr^s, and that the received flux is F{v) = AQ.sI{v). There are 2 
ways to measure 

1. from VLBI observations, one can often resolve the source and hence 
directly measure the angular size. In this case the relation between 
the brightness temperature measured in the K and K' frames is 

2. If the source is varying, we can estimate its size by requiring that the 
observed variability time-scale Atvar is longer than the light travel 
time -R/c, where R is the typical radius of the emission region. In this 
case 

C 



where c^a is the angular distance, related to the luminosity distance 
dL hy dA = dL/{l + zf. 
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There is a particular class of extragalactic radio sources, called Intra- 
Day Variable (IDV) sources, showing variability time-scales of hours in the 
radio band. For them, the corresponding observed brightess temperature 
can exceed 10^^ K, a value much larger than the theoretical limit for an 
incoherent synchrotron source, which is between 10^^ and 10^^ K. If the 
variability is indeed intrinsic, namely not produced by interstellar scintilla- 
tion, then one would derive a limit on the beamig factor 5, which should be 
larger than about 100. 

3.2.10 Moving in an homogeneous radiation field 

Jets in AGNs often moves in an external radiation field, produced by, e.g. 
the accretion disk, or by the Broad Line Region (BLR) which intercepts a 
fraction of the radiation produced by the disk and re-emits it in the form 
of emission lines. It it therefore interesting to calculate what is the energy 
density seen by a an observer which is comoving with the jet plasma. 




Figure 3.8: A real case: a relativistic bob is moving within the Broad Line 
Region of a radio loud AGN, with Lorentz factor F. In the rest frame K' of 
the blob the photons coming from 90° in frame K are seen to come at an 
angle 1/F. The energy density as seen by the blob is enhanced by a factor 

~ i . 
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To make a specific example, as illustrated by Fig. 13.81 assume that 
a portion of the jet is moving with a bulk Lorentz factor T, velocity /3c 
and that it is surrounded by a shell of broad line clouds. For simplicity, 
assume that the broad line photons are produced by the surface of a sphere 
of radius R and that the jet is within it. Assume also that the radiation is 
monochromatic at some frequency fo (in frame K). The comoving (in frame 
K') observer will see photons coming from a cone of semi-aperture 1 /F (the 
other half may be hidden by the accretion disk): photons coming from the 
forward direction are seen blue-shifted by a factor (1 + /3)F, while photons 
that the observer in K sees as coming from the side (i.e. 90° degrees) will 
be observed in K' as coming from an angle given by sin^' = l/F (and 
cos 6' = (3) and will be blue-shifted by a factor F. As seen in K\ each 
element of the BLR surface is moving in the opposite direction of the actual 
jet velocity, and the photons emitted by this element form an angle 6' with 
respect the element velocity. The Doppler factor used by K' is then 

The intensity coming from each element is seen boosted as (cfr Eq. I3.'2.10p : 

/' = (3.27) 

The radiation energy density is the integral over the solid angle of the in- 
tensity, divided by c: 

[/' = ^ C I'd cos e' 



2^ 



1 + /3 + 



F4(l -/3cos0')^ 



d cos 6' 



l + /3+y^F2C7 (3.28) 

Note that the limits of the integral correspond to the angles 0' and 90° in 
frame K. The radiation energy density, in frame K' , is then boosted by a 
factor (7/3)F2 when /? ~ 1. Doing the same calculation for a sphere, one 
would obtain U' = V'^U . 

Furthermore a (monochromatic) flux in K is seen, in K\ at different 
frequencies, between Tvq and (1 + /3)Ffo, with a slope F'{u') oc z^'^. Why 
the slope u'"^! This can be derived as follows: we already know that I'{i^') = 
5'^I{v) = {u' /u)^I{u). The flux at a specific frequency is 

F'{iy') = 2tt J"' dfi' (^^^ I{u) (3.29) 
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where = cos9', and the integral is over those fi' contributing at i/'. Since 

1 ,^1 



we have 

*' = -^ (3.31) 

Therefore, if the intensity is monochromatic in frame i.e. I{v) = Iq5{v — 
1^0 )) the flux density in the comoving frame is 



IU2 

271 L 



r/3 v^^ 



%'2 



Tz^o < i^' < (1 + /3)rz^o (3.32) 



where the frequency limits corresponds to photons produced in an emysphere 
in frame K, and between 0° and sin0' = l/F in frame K' . Integrating Eq. 
13.321 over frequency, one obtains 

F' = 27r/or2 (^1 + /3 + =Y^[l + p + ^^ F (3.33) 

in agreement with Eq. 13.281 



3.3 Superluminal motion 

In 1971 the Very Long Baseline Interferometry began, linking different radio- 
telescopes that where distant even thousands of km. The resolving power of 
a telescope is of the order of 

<t>- ^ (3.34) 

where A is the wavelength to be observed, and D is the diameter of the tele- 
scope or the distance of two connected telescopes. Observing at 1 cm, with 
two telescopes separated by 1000 km (i.e. 10^ cm), means that we can observe 
details of the source down to the milli-arcsec level (m.a.s.). The first obser- 
vations of the inner jet of radio-loud quasars revealed that the jet structure 
was not continuous, but blobby, with several radio knots. Repeating the ob- 
servations allowed us to discover that the blobs were not stationary, but were 
moving. Comparing radio maps taken at different times one could measure 
the angular displacement between the position of the blob. Knowing 
the distance d, one could then transform is a linear size: Ai? = 
Dividing by the time interval Ata between the two radio maps, one obtains 
a velocity 

dM 

^app = -r— (3.35) 




Figure 3.9: Top: The apparent velocity /3app as a function of the viewing 
angle 9 for different values of T, as labelled. Bottom: the amplification S'^ 
as a function of the viewing angle, for the same T as in the top panel. 



With some surprise, in several objects this turned out to be larger than 
the light speed c. Therefore these sources were called superluminal. The 
explanation of this apparent violation of special relativity is in Fig. 13. 4t the 
time interval Ata can be much shorter than the emission time Ate- With 
reference to Fig. 13.41 what the observer measures in the two radio maps is 
the position of the blob in point A (first map) and B (second map), projected 
in the plane of the sky. The observed displacement is then: 

CB = /3cAteSin(9 (3.36) 

Dividing by Ata = Atc(l — /3cos0) we have the measured apparent velocity 
as 

_ (3c Ate sin (9 _ (3 sin 9 

"""PP " Ata ^'^P " 1-/3 cos ^ ' 
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Ask yourself: F does not appear. Is it ok? At 0° the apparent velocity is zero. 
Is is ok? At what angle /3app is maximized? What is its maximum value? 
Fig. 13.91 shows /3app as a function of the viewing angle (angle between the 
line of sight and the velocity) for different F. Is the apparent superluminal 
speed given by a real motion of the emitting material? Can it be something 
else? If there are other possibilities, how to discriminate among them? 



u = v'6 


frequency 


t = t'/6 


time 


V = V'5 


volume 


sva.9 = sin 9' / 6 


sine 


cos 9 = (cos 0' + /?)/(! + /3 cos 9') 


cosine 




specific intensity 


I = 6^r 


total intensity 




specific emissivity 




absorption coefficient 


Tb = %5 


brightn. temp, (size directly measured) 


Tb = r^s^ 


brightn. temp, (size from variability) 


U' = (1 + /3 + /3V3)F2[7 


Radiation energy density within an emisphere 



Table 3.1: Useful relativistic transformations 



3.4 A question 

Suppose that some optically thin plasma of mass m is falling onto a central 
object with a velocity v and bulk Lorentz factor F. The central object 
has mass M and produces a luminosity L. Assume that the interaction is 
through Thomson scattering and that there are no electron-positron pairs. 

a) What is the radiation force acting on the electrons? 

b) What is the gravity force acting on the protons? 

c) What definition of limiting ( "Eddington" ) luminosity would you give 
in this case? 

d) What happens if the plasma is instead going outward? 



54 



CHAPTERS. BEAMING 



References 

Lind K.R. & Blandford R.D., 1985, ApJ, 295, 358 

Mook D.E. & Vargish T., 1991, Inside Relativity, Princeton Univ. Press 
Penrose R., 1959, Proc. of the Cambridge Philos. Soc, 55, p. 137 
Sikora M., Madejski G., Modcrski R. k Poutanen J., 1997, ApJ, 484, 108 
Terrel J., 1959, Phys. Rev., 116, 1041 



Chapter 4 

Synchrotron emission and 
absorption 

4.1 Introduction 

We now know for sure that many astrophysical sources are magnetized and 
have rclativistic Icptons. Magnetic field and rclativistic particles are the two 
ingredients to have synchrotron radiation. What is responsible for this kind 
of radiation is the Lorentz force, making the particle to gyrate around the 
magnetic field lines. Curiously enough, this force does not work, but makes 
the particles to accelerate even if their velocity modulus hardly changes. 
The outline of this section is: 

1. We will derive the total power emitted by the single electron. Total 
means integrated over frequency and over emission angles. This will 
require to generalize the Larmor formula to the rclativistic case; 

2. We will then outline the basics of the spectrum emitted by the single 
electron. This is treated in several text-books, so we will concentrate 
on the basic concepts; 

3. Spectrum from an ensemble of electrons. Again, only the basics; 

4. Synchrotron self absorption. We will try to discuss things from the 
point of view of a photon, that wants to calculate its survival proba- 
bility, and also the point of view of the electron, that wants to calculate 
the probability to absorb the photon, and then increase its energy and 
momentum. 

4.2 Total losses 

To calculate the total (=integrated over frequencies and emission angles) 
synchrotron losses we go into the frame that is instantaneously at rest with 
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the particle (in this frame is zero, but not the acceleration!). This is 
because we will use the fact that the emitted power is Lorentz invariant: 



2p2 9^2 

P - P' - — a'2 - — 



ail + a I 



(4.1) 



where the subscript "e" stands for "emitted". The fact that the power is 
invariant sounds natural, since after all, power is energy over time, and both 
energy and time transforms the same way (in special relativity with rulers 
and clocks). But be aware that this does not mean that the emitted and 
received power are the same. They are not! 

The problem is now to find how the parallel (to the velocity vector) and 
perpendicular components of the acceleration Lorentz transform. This is 
done in text books, so we report the results: 

aj| = T^ay 

a'^ = 7'a± (4.2) 

where 7 is the particle Lorentz factor. One easy way to understand and re- 
member these transformations is to recall that the acceleration is the second 
derivative of space with respect to time. The perpendicular component of 
the displacement is invariant, so we have only to transform (twice) the time 
(factor 7^). The parallel displacement instead transforms like 7, hence the 
7^ factor. 

The generalization of the Larmor formula is then: 



2e2 

P - P' - — 



12 I 12 

an + a_|_ 



2e^ 4 
3c3^ 



2 2, 2 
7 «|| + «± 



(4.3) 



Don't be fooled by the 7^ factor in front of a|... this component of the power 
is hardly important: since the velocity, for relativistic particles, is always 
close to c, it implies that one can get very very small acceleration in the same 
direction of the velocity. This is why linear accelerators minimize radiation 
losses. For synchrotron machines, instead, the losses due to radiation can be 
the limiting factor, and they are of course due to a±: changing the direction 
of the velocity means large accelerations, even without any change in the 
velocity modulus. To go further, we have to calculate the two components 
of the acceleration for an electron moving in a magnetic field. Its trajectory, 
in general, will have an helical shape of radius tl (the Larmor radius). The 
angle that the velocity vector makes with the magnetic field line is called 
pitch angle. Let us denote it with 9. We can anticipate that, in the absence 
of electric field and for a homogeneous magnetic field, the modulus of the 
velocity will not change: the magnetic field does not work, and so there is 
no change of energy, except for the losses due to the synchrotron radiation 
itself. So one important assumption is that at least during one gyration, the 
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Figure 4.1: A particle gyrates along the magnetic field lines. Its trajectory 
has an helicoidal shape, with Larmor radius tl and pitch angle 6. 



losses are not important. This is almost always satisfied in astrophysical 
settings, but there are indeed some cases where this is not true. 

When there is no electric field the only acting force is the (relativistic) 
Lorentz force: 

d e 
Fi^ = — (7mv) = - V X B (4.4) 
at c 

The parallel and perpendicular components are 

Fl|| = ev\\B = «|| = 

Fl^ = im— = e —B ^ a^ = (4.5) 

at c 'jmc 

We can also derive the Larmor radius tl by setting a± = v'^/ri,, and so 

v'^i 'ymc'^/3 sinO 

TL = — = ^ 4.6) 

a± eB 



The fundamental frequency is the inverse of the time occurring to complete 
one orbit (gyration frequency), so vb = c/3 sin 6 /(lirri^), giving 
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where z^l is the Larmor frequency, namely the gyration frequency for sub- 
relativistic particles. Larger B means smaller tl, hence greater gyration 
frequencies. Vice-versa, larger 7 means larger inertia, thus larger tl, and 
smaller gyration frequencies. Substituting a± given in Eq. 14.51 in the gener- 
alized Larmor formula (Eq. 14. 3p we get: 

= ^^^V/3'sin2 (4.8) 

We can make it nicer (for future use) by recalling that: 

• The magnetic energy density is Ub = B^/{8-k) 

• the quantity / {rrieCp), in the case of electrons, is the classical electron 
radius rg 

• the square of the electron radius is proportional to the Thomson scat- 
tering cross section ctt, i.e. ctt = Stttq/S = 6.65 x 10~^^ cm^. 

Making these substitutions, we have that the synchrotron power emitted by 
a single electron of given pitch angle is: 

Ps(f) = 2aTcUB7^/3^sm'^e (4.9) 

In the case of an isotropic distribution of pitch angles we can average the 
term sin^ 6 over the solid angle. The result is 2/3, giving 

(^'s) = ^(ttcUbJ^P^ (4.10) 

Now pause, and ask yourself: 

• Is Ps valid only for relativistic particles, or does it describe correctly 
the radiative losses also for sub-relativistic ones? 

• In the relativistic case the losses are proportional to the square of the 
electron energy. Do you understand why? And for sub-relativistic 
particles? 

• What happens of we have protons, instead of electrons? 

• What happens for — )• 0? Are you sure? (that losses vanishes..). Ok, 
but what happens to the received power when you have the lines of 
the magnetic field along the line of sight, and a beam of particles, all 
with a small pitch angles, shooting at you? 

• Why on earth there is the scattering cross section? Is this a coincidence 
or does it hide a deeper fact? 
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4.2.1 Synchrotron cooling time 

When you want to estimate a timescale of a quantity A, you can always 
write t = A/ A. In our case A is the energy of the particle. For electrons 
with an isotropic pitch angle distribution we have 

_ _E_ _ -fniec^ 7.75 x 10^ _ 24.57 

In the vicinity of a supermassive AGN black hole we can have B = W^B-^ 
Gauss and 7 = 10^73, yielding tgyn = 0.75/(i?|73) s. The same electron, in 
the radio lobes of a radio loud quasar with B = 10~^i?_5 Gauss, cools in 
tsyn = 246 million years. 



4.3 Spectrum emitted by the single electron 
4.3.1 Basics 

There exists a typical frequency associated to the synchrotron process. This 
is related to the inverse of a typical time. If the electron is relativistic, this 
is not the revolution period. Instead, it is the fraction of the time, for each 
orbit, during which the observer receives some radiation. To simplify, con- 
sider an electron with a pitch angle of 90°, and look at Fig. 14.21 illustrating 
the typical patterns of the produced radiation for sub-relativistic electrons 
moving with a velocity parallel (top panel) or perpendicular (mid panel) to 
the acceleration. In the bottom panel we see the pattern for a relativistic 
electron (with v _L a): it is strongly beamed in the forward direction. This is 
the direct consequence of the aberration of light, making half of the photons 
be emitted in a cone of semi-aperture angle I/7 (which is called the beaming 
angle). Note that this does not mean that half of the power is emitted within 
1/7, because the photons inside the beaming cone are more energetic than 
those outside, and are more tightly packed (do you remember the 5^ factor 
when studying beaming?). 1 

To go further, recall what we do when we study a time series and we 
want to find the power spectrum: we Fourier transform it. In this case we 
must do the same. Therefore if there is a typical timescale during which we 
receive most of the signal, we can say that most of the power is emitted at 
a frequency that is the inverse of that time. 

Look at Fig. 14.31 the relativistic electron emits photons all along its 
orbit, but it will "shoot" in a particular direction only for the time 

AB 1 -f 1 2mcv 2 1 2 1 , , 

At, = -r^l = - = = 4.12 

V V 2 V eB 2tt z^l 2tt ^v-q 

This is the emitting time during which the electron emits photons that will 
reach the observer. We can approximate the arc AB with a straight segment 
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Figure 4.2: Radiation patterns for a non relativistic particle with the velocity 
parallel (top) or perpendicular (mid) to the acceleration. When the particle 
is relativistic, the pattern strongly changes due to the aberration of light, 
and is strongly beamed in the forward direction. 



if the electron is relativistic, and the observed will then measure an arrival 
time At A that is shorter than Atg: 

AU = At. (1 - « = A,. ~ ^ = 2^ (4.13) 

The inverse of this time is the typical synchrotron (angular) frequency = 
27ri/s, so that: 

^ 7^1^ = = 7^ ^ (4.14) 



27rAtA 27rmeC 

This is a factor 7^ greater than the fundamental frequency, and a factor 7^ 
greater than the Larmor frequency, defined as the typical frequency of non- 
relativistic particles. We expect that the particle emits most of its power at 
this frequency. 



4.3.2 The real stuff 

One can look at any text book for a detailed discussion of the procedure 
to calculate the spectrum emitted by the single particle. Here we report 
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AB = 2rjy 



Atp =AB/c 



line of sight 
► 



Figure 4.3: A relativistic electron is gyrating along a magnetic field line 
with pitch angle 90°. Its trajectory is then a circle of radius tl- Due to 
aberration, an observer will "see it" (i.e. will measure an electric field) 
when the beaming cone of total aperture angle 2/7 is pointing at him. 



the results: the power per unit frequency emitted by an electron of given 
Lorentz factor and pitch angle is: 



— sm ( 
2 



K5/3iy)dy 



(4.15) 



This is the power integrated over the emission pattern. K^i^[y) is the 
modified Bessel function of order 5/3. The dependence upon frequency 
is contained in F{v/vc), that is plotted in Fig. 14. 4i This function peaks at 
u ~ 0.29z^c5 therefore very close to what we have estimated before, in our 
very approximate treatment. The low frequency part is well approximated 
by a power law of slope 1/3: 



47r f u 
^ ' ' ^/3r(l/3) \2vJ ^ 



(4.16) 



Ki V ^ Vc the function decays exponentially, and can be approximated by: 

1/2 



7r\i/2 
2. 



— e 



(4.17) 
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Figure 4.4: Top panel: The function F{u/uc) describing the synchrotron 
spectrum emitted by the single electron. Bottom panel: F{v/vc) is com- 
pared with some approximating formulae, as labeled. We have defined 
X = V jvc- 
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Another approximation valid for all frequency, but overestimating F around 
the peak, is: 

47r / V^^ / 

Fiiy/i^c) 7= e-"/"^ (4.18) 

^ ' ' V3r(l/3) \2vJ ^ ' 

4.3.3 Limits of validity 

One limit can be obtained by requiring that, during one orbit, the emitted 
energy is much smaller than the electron energy. If not, the orbit is modified, 
and our calculations are no more valid. For non-relativistic electrons this 
translates in demanding that 

huB < m^c^ B < — j-^ = Be (4.19) 

where -Be ~ 4.4 x 10^^ Gauss is the critical magnetic field, around and above 
which quantum effects appears (i.e. quantized orbits. Landau levels and so 
on). 

For relativistic particles we demand that the energy emitted during one 
orbit does not exceed the energy of the particle. 

Ps 9 e/fjT 7.22 X 10^^ ^ 

— < -fniec^ ^ B < „ . 2^ Gauss 4.20 

z/B Y sm 9 Y sm 9 

Therefore for large 7 we reach the validity limit even if the magnetic field is 
sub-critical. 

For very small pitch angles beware that the spectrum is not described 
by F[u/uc), but consists of a blue-shifted cyclotron line. This is because, in 
the gyroframe, the particle is sub-relativistic, and so it emits only one (or 
very few) harmonics, that the observer sees blueshifted. 



4.3.4 Prom cyclotron to synchrotron emission 

A look at Fig. 14.51 helps to understand the difference between cyclotron 
and synchrotron emission. When the particle is very sub-relativistic, the 
observed electric field is sinusoidal in time. Correspondingly, the Fourier 
transform of E[t) gives only one frequency, the first harmonic. Increasing 
somewhat the velocity (say, /3 ~ 0.01) the emission pattern starts to be 
asymmetric (for light aberration) and as a consequence E{t) must be de- 
scribed by more than just one sinusoid, and higher order harmonics appear. 
In these cases the ratio of the power contained in successive harmonics goes 
as /32. 

Finally, for relativistic (i.e. 7^1) particles, the pattern is so asym- 
metric that the observers sees only spikes of electric field. They repeat 
themselves with the gyration period, but all the power is concentrated into 
AtA- To reproduce E{t) in this case with sinusoids requires a large number 
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Figure 4.5: Prom cyclo to synchro: if the emitting particle has a very small 
velocity, the observer sees a sinusoidal (in time) electric field E{t). Increas- 
ing the velocity the pattern becomes asymmetric, and the second harmonic 
appears. For < /3 <C 1 the power in the second harmonic is a factor less 
than the power in the first. For relativistic particles, the pattern becomes 
strongly beamed, the emission is concentrated in the time At a- As a con- 
sequence the Fourier transformation of E{t) must contain many harmonics, 
and the power is concentrated in the harmonics of frequencies z/ ~ 1/AtA- 
Broadening of the harmonics due to several effects ensures that the spec- 
trum in this case becomes continuous. Note that the fundamental harmonic 
becomes sma//er increasing 7 (since ub oc I/7). 
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of them, with frequencies going at least up to l/AtA- In this case the har- 
monics are many, guaranteeing that the spectrum becomes continuous with 
any reasonable line broadening effect, and the power is concentrated at high 
frequencies. 



4.4 Emission from many electrons 

Again, this problem is treated in several text books, so we repeat the basic 
results using some approximations, tricks and shortcuts. 

The queen of the particle energy distributions in high energy astrophysics 
is the power law distribution: 

dE 

Ar(7) = K-f-P = N{E) — ; j^^^j < (4.21) 

(17 

Now, assuming that the distribution of pitch angles is the same at low and 
high 7, we want to obtain the synchrotron emissivity produced by these par- 
ticles. Beware that the emissivity is the power per unit solid angle produced 
within 1 cm'^. The specific emissivity is also per unit of frequency. So, if Eq. 
14.211 represents a density, we should integrate over 7 the power produced by 
the single electron (of a given 7) times A^(7), and divide all it by 47r, if the 
emission is isotropic: 

es{^,0) = — I N{^)P{^,u,9)dj (4.22) 



mm 



Doing the integral one easily finds that, in an appropriate range of frequen- 
cies: 

es{iy,9) oc i^s{p+i)/2j,-{p-i)/2 (4^23) 

The important thing is that a power law electron distribution produces a 
power law spectrum, and the two spectral indices are related. We tradition- 
ally call Q the spectral index of the radiation, namely eg oc We then 
have 

a = ^ (4.24) 

This result is so important that it is worth to try to derive it in a way as 
simple as possible, even without doing the integral of Eq. 14.221 We can 
in fact use the fact that the synchrotron spectrum emitted by the single 
particle is peaked. We can then say, without being badly wrong, that all 
the power is emitted at the typical synchrotron frequency: 

= rvh; = (4.25) 

zvrmeC 

In other words, there is a tight correspondence between the energy of the 
electron and the frequency it emits. To simplify further, let us assume that 
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the pitch angle is 90° . The emissivity at a given frequency, within an interval 
du, is then the result of the emission of electrons having the appropriate 
energy 7, within the interval d'j 

1 r jyV^^ d-f jy-V2 

= _p.«(-,)d,; y ; _2=_ (4.26) 

L 

we then have 

du 

OC ifB(p+i)/2^-(P-i)/2 (427) 

where we have used uj^ (x. B. 

The synchrotron flux received from a homogeneous and thin source of 
volume y oc i?^, at a distance cZl, is 

F,{v) = ATTe,{u) ^ 



oc ^KB^+'^v-^ 

oc elRKB^+'^v-'^ (4.28) 

where Og is the angular radius of the source (not the pitch angle!). Observing 
the source at two different frequencies allows to determine a, hence the slope 
of the particle energy distribution. Furthermore, if we know the distance 
and -R, the normalization depends on the particle density and the magnetic 
field: two unknowns and only one equation. We need another relation to 
close the system. As we will see in the following, this is provided by the 
self-absorbed flux. 



4.5 Synchrotron absorption: photons 

All emission processes have their absorption counterpart, and the synchrotron 
emission is no exception. What makes synchrotron special is really the fact 
that it is done by relativistic particles, and they are almost never distributed 
in energy as a Maxwellian. If they were, we could use the well known fact 
that the ratio between the emissivity and the absorption coefficient is equal 
to the black body (Kirchhoff law) and then wc could easily find the absorp- 
tion coefficient. But in the case of a non-thermal particle distribution we 
cannot do that. Instead we are obliged to go back to more fundamental 
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relations, the one between the A and B Einstein coefficients relating spon- 
taneous and stimulated emission and "true" absorption (by the way, recall 
that the absorption coefficient is what remains subtracting stimulated emis- 
sion from "true" absorption). But we once again will use some tricks, in 
order to be as simple as possible. These are the steps: 

1. The first trick is to think to our power law energy distribution as a 
superposition of Maxwellians, of different temperatures. So, we will 
relate the energy jrUeC^ of a given electron to the energy kT of a 
Maxwellian. 

2. We have already seen that there is a tight relation between the emitted 
frequency and 7. Since the emission and absorption processes are 
related, we will assume that a particular frequency v is preferentially 
absorbed by those electrons that can emit it. 

3. As a consequence, we can associate our "fake" temperature to the 
frequency: 

/ J, \ 1/2 

kT - jruec^ ~ mec^ ( — j (4.29) 

4. For an absorbed source the brightness temperature Ti„ defined by 

I{v) = 2kn— (4.30) 
must be equal to the kinetic "temperature" of the electrons, and so 

I{u) = '^kT-^ ~ 2m^v^ [ — 

^ ^ (4-31) 



These are the right dependencies. Note that the spectrum is oc i^^/^, not 
i^^, and this is the consequence of having "different temperatures". Note 
also that the particle density disappeared: if you think about it is natural: 
the more electrons you have, the more you emit, but the more you absorb. 
Finally, even the slope of the particle distribution is not important, it controls 
(up to a factor of order unity) only the normalization of I{u) (our ultra- 
simple derivation cannot account for it, see the Appendix). 

The above is valid as long as we can associate a specific 7 to any v. 
This is not always the case. Think for instance to a cut-off distribution, 
with 7min ^1- In this case the electrons with 7min are the most efficient 
emitters and absorbers of all photons with v < I'min = 7^in^'L- So in this case 
we should not associate a different temperature when dealing with different 
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Figure 4.6: The synchrotron spectrum from a partially self absorbed source. 
Observations of the self absorbed part could determine B. Observations of 
the thin part can then determine K and the electron slope p. 



V < ^min- But if do not change T, we recover a self-absorbed intensity 
I{i') oc i^^ (i.e. Raleigh- Jeans like). 

Now, going from the intensity to the flux, we must integrate /(z^) over 
the angular dimension of the source (i.e. 9s), obtaining 

F{v) OC 02 __ (4.32) 

if we could observe a self-absorbed source, of known angular size, we could 
then derive its magnetic field even without knowing its distance. 



4.5.1 Prom thick to thin 

To describe the transition from the self absorbed to the thin regime we have 
to write the radiation transfer equation. The easiest one is for a slab. Calling 
Ky the specific absorption coefficient [cm~^] we have 

/(j,) = !M(i_e--); Ty = RKy (4.33) 

it is instructive to write Eq. 14.331 in the form: 

I{v) = e{v)R (4.34) 

because in this way it is evident that when ^ 1 (self absorbed regime), 
we simply have 

/(.) = ii^ = iM; ry»l (4.35) 
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Figure 4.7: The synchrotron absorption cross section as a function of z//z/l 
for different values of 7, as labeled, assuming a pitch angle oi 9 = 90° and a 
magnetic field of 1 Gauss. 



One can interpret it saying that the intensity is coming from electrons lying 
in a shell within R/ti, from the surface. 

Since we have already obtained /(i/) oc v^l'^B^^I'^ in the absorbed regime, 
we can derive the dependencies of the absorption coefficient: 

- J(j,) ^ 1,5/25-1/2 - ^ ^'^••^'^^ 

Note the rather strong dependence upon frequency: at large frequencies, 
absorption is small. 

The obvious division between the thick and thin regime is when t^, = \. 
We call self- absorption frequency, vt, the frequency when this occurs. We 
then have: 



RKly^ = 1 — )■ f t oc 



2/(p+4) 

(4.37) 



The self-absorption frequency is a crucial quantity for studying synchrotron 
sources: part of the reason is that it can be thought to belong to both regimes 
(thin and thick), the other reason is that the synchrotron spectrum peaks 
very close to ut (see Fig. 14. 6p even if not exactly at Uf (see the Appendix) . 



4.6 Synchrotron absorption: electrons 

In the previous section we have considered what happens to the emitted 
spectrum when photons are emitted and absorbed. This is described by 
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the absorption coefficient. But now imagine to be an electron, that emits 
and absorbs synchrotron photons. You would probably be interested if your 
budget is positive or negative, that is, if you are loosing or gaining energy. 
This is most efficiently described by a cross section, that tells you the prob- 
ability to absorb a photon. Surprisingly, the synchrotron absorption cross 
section has been derived relatively recently (Ghisellini and Svensson 1991), 
and its expression is: 

a^iu, ^,6) = ^^ (^^) (4.38) 

3a/3 ij sm 6^ ' \ Uc sm 6 J 

For frequencies v -^Vc this expression can be approximated by: 

87r2(3sin^)2/3r(5/3) e / v \ "^/^ , 

Note these features: 

• At the fundamental frequency 1^/7, the cross section does not depend 
on 7. 

• The dimensions are given by e/B: this factor is proportional to the 
product of the classical electron radius and the Larmor wavelength (or 
radius). Imagine an electron with 90° pitch angle, and to see its orbit 
from the side: you would see a rectangle of base tl and height tq. The 
area of this rectangle is of the order of e/B. At low frequencies, 
can be orders of magnitudes larger than the Thomson scattering cross 
section. 

• There is no explicit dependence on the particle mass. However, protons 
have much smaller z/l, and the dependence on mass is hidden there. 
Nevertheless, electrons and protons have the same cross section (of 
order e/B) at their respective fundamental frequencies. 

Fig. 14.71 shows function of v/vi^ for different 7. The thing it 

should be noticed is that this cross section is really large. Can we make 
some useful use of it? Well, there are at least two issues, one concerning 
energy, and the other concerning momentum. 

First, electrons emitting and absorbing synchrotron photons do so with 
a large efficiency. They can talk each other by exchanging photons. There- 
for, even if they are distributed as a power law in energy at the beginning, 
they will try to form a Maxwellian. They will form it, as long as other 
competing processes are not important, such as inverse Compton scatter- 
ings. The formation of the Maxwellian will interest only the low energy part 
of the electron distribution, where absorption is important. Note that this 
thermalization process works exactly when Coulomb collisions fail: they are 
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inefficient at low density and high temperature, while synchrotron absorp- 
tion can work for relativistic electrons even if they are not very dense. 

The second issue concerns exchange of momentum between photons and 
electrons. Suppose that a magnetized region with relativistic electrons is 
illuminated by low frequency radiation by another source, located aside. 
The electrons will efficiently absorb this radiation, and thus its momentum. 
The magnetized region will then accelerate. 
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4.7 Appendix: Useful Formulae 

In this section we collect several useful formulae concerning the synchrotron 
emission. When possible, we give also simplified analytical expressions. We 
will often consider that the emitting electrons have a distribution in energy 
which is a power law between some limits 71 and 72. Electrons are assumed 
to be isotropically distributed in the comoving frame of the emitting source. 
Their density is 



N{'j) = 71 < 7 < 72 



(4.40) 



The Larmor frequency is defined as: 



eB 



(4.41) 



'In in, ( 



4.7.1 Emissivity 



The synchrotron emissivity es{u,6) [erg cm 



^ s -"^ sterad is 
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(4.42) 



where Ps(z^, 7, 0) is the power emitted at the frequency u (integrated over 
all directions) by the single electron of energy 'ynieC^ and pitch angle 6. 
For electrons making the same pitch angle 9 with the magnetic field, the 
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emissivity is 



J. ( 3p-A p ( 3p+19 \ 

SaTcKUs ( ^ \ ^ / . ^\2±ig£ W2 y^l^ 12 y 



between z^i ^ 7i^'l and ^^2 ^ l2^L ■ If the distribution of pitch angles is 



isotropic, we must average the (sin 6) 2 term, obtaining 



„2L ^^^2±5^ 

<(sin^)— >= / (sin6')— sin^d^ = ^^^^ ) f (4.44) 

Therefore the pitch angle averaged synchrotron emissivity is 

p-i 



The function /e(p) includes all the products of the F-functions: 

3i r (¥) r (31^) r (^) 



feip) = 



- 3i(p+ 0.105) (4.46) 
where the simplified fitting function is accurate at the per cent level. 

4.7.2 Absorption coefficient 

The absorption coefficient K,y{0) [cm~^] is defined as: 

Written in this way, the above formula is valid even when the electron distri- 
bution is truncated. For our power law electron distribution Ki,{9) becomes: 

Above = ^fi'L, we have: 

""^^'^ = W^'^^""') ^.^32 ri^-^j ^ 

Averaging over the pitch angles we have: 
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resulting in a pitch angle averaged absorption coefficient: 

^e^K £±2 _p+4 
K. = ^.,^ . ^ AW (4.51) 

where the function f^ip) is: 

■p f 3p+22 \ p / 3p+2 \ p / p+6 \ 

r 

^P+i / 1.8 
~ 3~ h — 

The simple fitting function is accurate at the per cent level. 




4.7.3 Specific intensity 

Simple radiative tranfer allows to calculate the specific intensity: 

/(j,) = f^M (1 _ e---) (4.53) 

where the absorption optical depth t^, = k^R and R is the size of the emitting 
region. When Ti, S> 1, the esponential term vanishes, and the intensity is 
simply the ratio between the emissivity and the absorption coefficient. This 
is the self-absorbed, ot thick, regime. In this case, since both es(z^) and 
Ki, depends linearly upon K, the resulting self-absorbed intensity does not 
depend on the normalization of the particle density K: 

m = -^^fiip) {l-e-^-^) (4.54) 

we can thus see that the slope of the self-absorbed intensity does not depend 
on p. Its normalization, however, does (albeit weakly) depend on p through 
the function flip), which in the case of averaging over an isotropic pitch 
angle distribution is given by: 



flip) 



r(¥ r r r 1^ 



p + i r (3^) r (3£±2) r (2±^) r(2+6) 



where again the simple fitting function is accurate at the level of 1 per cent. 
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4.7.4 Self-absorption frequency 

The self-absorption frequency ft is defined by Tl,^ = 1: 



Tie^RK 



Up) 



4 
p+4 



tt-v/tF eRK 



B 



Up) 



2 

p+4 



(4.56) 



Note that the term in parenthesis is adimensional, and since RK has units 
of the inverse of a surface, then e/B has the dimension of a surface. In 
fact we have already discussed that this is the synchrotron absorption cross 
section of a relativistic electron of energy ^rrieC? absorbing photons at the 
fundamental frequency vl/i- 

The random Lorentz factor 7t of the electrons absorbing (and emitting) 
photons with frequency Vi is 7t ~ [3ft /(4z//,)]^/^. 



4.7.5 Synchrotron peak 

In a F{v) plot, the synchrotron spectrum peaks close to ft) at a frequency 
fs,p given by solving 



dl(y) 
dv 



dv 



.5/2 (l-e--") 







which is equivalent to solve the equation: 



e^P V^s,v) ^^^s,P -1 = 



whose solution can be approximated by 

2 



5 



(4.57) 



(4.58) 



(4.59) 



Chapter 5 



Compton scattering 



5.1 Introduction 

The simplest interaction between photons and free electrons is scattering. 
When the energy of the incoming photons (as seen in the comoving frame 
of the electron) is small with respect to the electron rest mass-energy, the 
process is called Thomson scattering, which can be described in terms of 
classical electro-dynamics. As the energy of the incoming photons increases 
and becomes comparable or greater than mgC^, a quantum treatment is 
necessary (Klein-Nishina regime). 

5.2 The Thomson cross section 

Assume an electron at rest, and an electromagnetic wave of frequency u <^ 
rUeC? /h. Assume also that the incoming wave is completely linearly po- 
larized. In order to neglect the magnetic force (e/c)(v x B) we must also 
require that the oscillation velocity v ^ c. This in turn implies that the 
incoming wave has a sufficiently low amplitude. The electron starts to os- 
cillate in response to the varying electric force eE, and the average square 
acceleration during one cycle of duration T = 1 /z/ is 



The emitted power per unit solid angle is given by the Larmor formula 
dP/dCl = e^a^sin^ 0/(47rc^) where 6 is the angle between the acceleration 

vector and the propagation vector of the emitted radiation. Please note 
that Q is not the scattering angle, which is instead the angle between the 
incoming and the scattered wave (or photon). We then have 




(5.1) 



dPe 

dfl 



sin^e 



(5.2) 
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The scattered radiation is completely linearly polarized in the plane defined 
by the incident polarization vector and the scattering direction. The fiux of 
the incoming wave is Si = cEI/{Stx). The differential cross section of the 
process is then 

where ro = /{rrif.cp') is the classic electron radius^ tq = 2.82 x 10""^^ cm. We 
see that the scattered pattern of a completely polarized incoming wave is a 
torus, with axis along the acceleration direction. The total cross section can 
be derived in a similar way, but considering the Larmor formula integrated 
over the solid angle [P = 2e^a^/(3c^)]. In this way the total cross section is 

(^poi = 5: = Y^o (5-4) 

Note that the classical electron radius can also be derived by equating the 
energy of the associated electric field to the electron rest mass-energy: 



/•oo roc 

me(? = / — Airr^dr = / 



roo rp2 rco 2 -1 2 



Why is ao slightly different from ro? Because there is an intrinsic uncertainty 
related to the distribution of the charge within (or throughout the surface 
of) the electron. See the discussion in Vol. 2, chapter 28.3 of "The Feynman 
Lectures on Physics" , about the fascinating idea that the mass of the electron 
is all electromagnetic. 

5.2.1 Why the peanut shape? 

The scattering of a completely unpolarized incoming wave can be derived 
by assuming that the incoming radiation is the sum of two orthogonal com- 
pletely linearly polarized waves, and then summing the associated scattering 
patterns. Since we have the freedom to chose the orientations of the two 
polarization planes, it is convenient to chose one of these planes as the one 
defined by the incident and scattered directions, and the other one perpen- 
dicular to this plane. The scattering can be then regarded as the sum of 
two independent scattering processes, one with emission angle B, the other 
with 7r/2. If we note that the scattering angle (i.e. the angle between the 
scattered wave and the incident wave) is^ = 7r/2 — 0, we have 



1 



do_\ 

^Junpol 2 
1 



^r2(l+sin2e) 
^rlil + cos'^e) (5.6) 




Figure 5.1: Photons are coming along the y-axis. The top panels shows the 
pattern of the scattered radiation for photons completely linearly polarized 
along the z-axis (left) and along the x-axis (right). The sum of the two 
torii corresponds to the pattern for unpolarized radiation (bottom panel). 
This explains why we have a "peanut" shape, elongated along the velocity 
vector of the incoming photons. Courtesy of Davide Lazzati. 
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In this case we see that the cross section depends only on the scattering 
angle 9. The pattern of the scattered radiation is then the superposition 
of two orthogonal "tori" (one for each polarization direction), as illustrated 
in Fig. I5.1[ When scattering completely linearly polarized radiation, only 
one "torus" survives. Instead, when scattering unpolarized radiation, some 
polarization is introduced, because of the difference between the two "tori" 
patterns. Both terms of the RHS of Eq. 15.61 refer to completely polarized 
scattered waves (but in two perpendicular planes). The difference between 
these two terms is then associated to the introduced polarization, which is 
then 

^ 1 - cos^ 9 . ^. 

" = TT^9 ^'-'^ 

The above discussion helps to understand why the scattering process intro- 
duces some polarization, which is maximum (100%) if the angle between the 
incoming and the scattered photons is 90° (only one torus contributes), and 
zero for 0° and 180°, where the two torii give the same contribution. 

The total cross section, integrated over the solid angle, is the same as 
that for polarized incident radiation (Eq. 15. 4p since the electron at rest has 
no preferred defined direction. This is the Thomson cross section: 

aT = [(^) dn=^-^ [il + cos^9)dcos9 ^ 

J V"^'/unpol 2 J 



3 
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6.65 X 10'^^ cm^ (5. 



5.3 Direct Compton scattering 

In the previous section we considered the scattering process as an interaction 
between an electron and an electromagnetic wave. This required hi' <C 
rrif.cp'. In the general case the quantum nature of the radiation must be 
taken into account. We consider then the scattering process as a collision 
between the electron and the photon^ and apply the conservation of energy 
and momentum to derive the energy of the scattered photon. It is convenient 
to measure energies in units of mgC^ and momenta in units of mgC. 

Consider an electron at rest and an incoming photon of energy xq, which 
becomes xi after scattering. Let 9 be the angle between the incoming and 
outgoing photon directions. This defines the scattering plane. Momentum 
conservation dictates that also the momentum vector of the electron, after 
the scattering, lies in the same plane. Conservation of energy and conserva- 
tion of momentum along the x and y axis gives: 

" l + xo(l-cos0) ^^-^^ 

Note that, for xq ^» 1 and cos^ / 1, xi — >• (1 — cos^)^^. In this case the 
scattered photon carries information about the scattering angle, rather than 
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about the initial energy. As an example, for 9 = tt and xq ^ 1, the final 
energy is xi = 0.5 (corresponding to 255 keV) independently of the exact 
value of the initial photon energy. Note that for xq <C 1 the scattered energy 
xi — xq, as assumed in the classical Thomson scattering. The energy shift 
implied by Eq. 15.91 is due to the recoil of the electron originally at rest, and 
becomes significant only when xq becomes comparable with 1 (or more). 
When the energy of the incoming photon is comparable to the electron rest 
mass, another quantum effect appears, namely the energy dependence of the 
cross section. 



5.4 The Klein— Nishina cross section 



The Thomson cross section is the classical limit of the more general Klein- 
Nishina cross section (here we use x as the initial photon energy, instead of 
Xq, for simplicity): 



3 fxi 

lOTT V X 



\ sm^ I 

Xl X 



(5.10) 



This is a compact form, but there appears dependent quantities, as sinO is 
related to x and xi. By inserting Eq. 15.91 we arrive to 



167r [l + x(l -cos6')]^ 



x(l — cos 9) + 



l + x(l - cos 6*) 



+ cos 



(5.11) 

In this form, only independent quantities appear (i.e. there is no xi). Note 
that the cross section becomes smaller for increasing x and that it coincides 
with daT:/dfl for 9 = (for this angle xi = x independently of x). This 
however corresponds to a vanishingly small number of interactions, since 
dn-^Ofor 9-^ 0). 

Integrating Eq. 15.111 over the solid angle, we obtain the total Klein- 
Nishina cross section: 



CRN 



1 + X 



2x(l + x) 
l + 2x 



ln(l + 2x) 



H ln(l + 2x) - 

2x ^ ^ (1 + 2x)2 

(5.12) 



Asymptotic limits are: 



UKN — ctt ^1 - 2x + 
3 (Tt 



26x' 



(TKN — 



ln(2x) + 
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X <C 1 
X > 1 



(5.13) 



The direct Compton process implies a transfer of energy from the photons to 
the electrons. It can then be thought as an heating mechanism. In the next 





Figure 5.3: The differential Klein-Nishina cross section (in units of cjt), 
for different incoming photon energies. Note how the scattering becomes 
preferentiahy forward as the energy of the photon increases. 
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Figure 5.4: Scattered photons energies as a function of the scattering angle, 
for different incoming photon energies. Note that, for x ^ 1 and for large 
scattering angle, the scattered photon energies becomes xi ~ 1/2, indepen- 
dent of the initial photon energy x. 

subsection we discuss the opposite process, called inverse Compton scatter- 
ing, in which hot electrons can transfer energy to low frequency photons. 

We have so far neglected the momentum exchange between radiation 
and the electron. One can see, even classically, that there must be a net 
force acting along the direction of the wave if one considers the action of 
the magnetic field of the wave. In fact the Lorentz force ev x B is always 
directed along the direction of the wave (here v is the velocity along the E 
field). This explains the fact that light can exert a pressure, even classically. 

5.4.1 Another limit 

We have mentioned that, in order for the magnetic Lorentz force to be neg- 
ligible, the electron must have a transverse (perpendicular to the incoming 
wave direction) velocity <C c. Considering a wave of frequency uj and electric 
field E = Eq sm{ujt), this implies that: 



This means that the scattering process can be described by the Thomson 
cross section if the wave have a sufficiently low amplitude and a not too 
small frequency (i.e. for very small frequencies the electric field of the wave 
accelerates the electron for a long time, and then to large velocities). 




(5.14) 
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5.4.2 Pause 

Now pause, and ask if there are some ways to apply what we have done up 
to now to real astrophysical objects. 

• The Eddington luminosity is derived with the Thomson cross section, 

with the thought that it describes the smallest probability of inter- 
action between matter and radiation. But the Klein-Nishina cross 
section can be even smaller, as long as the source of radiation emits 
at high energy. What are the consequences? If you have forgotten the 
definition of the Eddington luminosity, here it is: 

LEdd = ^ = 1.3 X lO^s — erg s'^ (5.15) 

(Tt Mq 

• In Nova Muscae, some years ago a (transient) annihilation line was 
detected, together with another feature (line-like) at 200 keV. What 
can this feature be? 

• It seems that high energy radiation can suffer less scattering and there- 
fore can propagate more freely through the universe. Is that true? Can 
you think to other processes that can kill high energy photons in space? 

• Suppose to have an astrophysical source of radiation very powerful 
above say - 100 MeV. Assume that at some distance there is a very 
efficient "reflector" (i.e. free electrons) and that you can see the scat- 
tered radiation. Can you guess the spectrum you receive? Does it 
contain some sort of "pile— up" or not? Will this depend upon the 
scattering angle? 



5.5 Inverse Compton scattering 

When the electron is not at rest, but has an energy greater that the typical 
photon energy, there can be a transfer of energy from the electron to the 
photon. This process is called inverse Compton to distinguish it from the 
direct Compton scattering, in which the electron is at rest, and it is the 
photon to give part of its energy to the electron. 

We have two regimes, that are called the Thomson and the Klein-Nishina 
regimes. The difference between them is the following: wc go in the frame 
where the electron is at rest, and in that frame we calculate the energy of the 
incoming photon. If the latter is smaller than mgC^ we are in the Thomson 
regime. In this case the recoil of the electron, even if it always exists, is 
small, and can be neglected. In the opposite case (photon energies larger 
than rUeCp), we are in the Klein-Nishina one, and we cannot neglect the 
recoil. As we shall see, in both regimes the typical photon gain energy, even 
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if there will always be some arrangements of angles for which the scattered 
photon looses part of its energy. 

5.5.1 Thomson regime 

Perhaps, a better name should be "inverse Thomson" scattering, as will 
appear clear shortly. 

5.5.2 Typical frequencies 

In the frame K' comoving with the electron, the incoming photon energy is 

x' = X7(l — /3 cos V') (5.16) 

where ip is the angle between the electron velocity and the photon direction 
(see Fig lS.Sp . 



— 

V 

Lab Frame 



Figure 5.5: In the lab frame an electron is moving with velocity v. Its 
velocity makes an angle ip with an incoming photon of frequency z/. In the 
frame where the electron is at rest, the photon is coming from the front, 
with frequency z^', making an angle ip' with the direction of the velocity. 

At first sight this is different from x' = x6 derived in Chapter 3. But 
notice that i) in this case the angle ip is measured in the lab frame; ii) it is 
not the same angle going into the definition of 6 (i.e. in S we use the angle 
between the line of sight and the velocity of the emitter, i.e. 9' = ir — ip'). 



V 



¥ 



Rest Frame K' 



84 



CHAPTER 5. COMPTON SCATTERING 



Going to the rest frame of the electrons we should use (recalling Eq. 3.16 
for the transformation of angles): 

/ ^ + cosV'' 
cosip = (5.17) 

Substituting this into equation 15.161 we have 

x' = — § — (5.18) 

7(1 + /3cos^/;') ^ ' 

Finally, consider that cos^' = cos(7r — = — cos^'', validating x' = x5. 

If x' <C 1, we are in the Thomson regime. In the rest frame of the 
electron the scattered photon will have the same energy x'l as before the 
scattering, independent of the scattering angle. Then 

x[ = x' (5.19) 

This photon will be scattered at an angle ip'i with respect to the electron 
velocity. The pattern of the scattered radiation will follow the pattern of 
the cross section (i.e. a peanut). Think to the scattering in the comoving 
frame as a re-isotropization process: even if the incoming photons are all 
coming from the same direction, after the scattering they are distributed 
quasi-isotropically. Going back to K the observer sees 

xi = x[-f{l + /3cosV'i) (5-20) 

Recalling again Eq. 3.16, for the transformation of angles: 

, cosV'i-/? 
1 — p cos tpi 

we arrive to the final formula: 

,,.,lz|£2Sf (5.22) 

1 — p COS tpi 

Now all quantities are calculated in the lab-frame. 

Let us see the minimum and maximum energies. The maximum is when 
= IT (head on collision), and when ipi = (the photon is scattered along 
the electron velocity vector). In these head-on collisions: 

xi = = 7^(1 + /3)2x ^ 47^; head-on (5.23) 

where the last step is valid if 7 ^ 1. The other extreme is for ■0i = vr and 
^ = 0. In this case the incoming photon "comes from behind" and bounces 
back. In these "tail-on" collisions: 

^^ = ^TT| = 7^(i + /3)2^i^' ^'-''^ 
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Figure 5.6: Maximum and minimum scattered frequencies. The maximum 
occurs for head-on colhsions, the minimum for tail-on ones. These two 
frequencies are one the inverse of the other. 

where again the last step is valid if 7 ^ 1. Another typical angle is s'mipi = 
1/7, corresponding to cos^/^i = /3. This corresponds to the aperture angle 
of the beaming cone. For this angle: 

1 — /3 cos Ip 2/ n :\ ^ 

xi = — - — — X = 7 (1 — p cos V'ja;; beaming cone (5.25) 

which becomes xi = x/{l + /3) for ip = 0, xi = 7^x for ip = 7r/2 and 
xi = 7^(1 + /3)x for tp = IT. 

For an isotropic distribution of incident photons and for 7^1 the 
average photon energy after scattering is (see Eq. I5.47p : 



Total loss rate 

We can simply calculate the rate of scatterings per electron considering all 
quantities in the lab-frame. Let n(e) be the density of photons of energy 
e = hv, V the electron velocity and ip the angle between the electron velocity 
and the incoming photon. For mono-directional photon distributions, we 
have: 



/ \ ^2 



(5.26) 




(5.27) 
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Vrei = c — v COS IS the relative velocity between the electron and the incom- 
ing photons. We then have 

^ = J aTc{l- /3 cos ip)n{€)de (5.28) 

Note that the rate of scatterings in the lab frame, when the electron and/or 
photon are anisotropically distributed, can be described by an effective cross 
section (Teg = / o"t(1 — l3costp)dCl/ATT. For photons and electrons moving 
in the same direction the scattering rate (hence, the effective optical depth) 
can be greatly reduced. 

The power contained in the scattered radiation is then 

dE^ eidN f (1 - j3 cos 4^ f , , 

- - CTTC ^ --^en{e)de (5.29) 



dt dt J 1 — P cos il^i 

Independently of the incoming photon angular distribution, the average 
value of 1 — j3cosij)i can be calculated recalling that, in the rest frame of 
the electron, the scattering has a backward-forward symmetry, and there- 
fore {cosip'i) = 7r/2. The average value of cosV'i is then /3, leading to 
(1 — /3cos^i) = 1/7^. We therefore obtain 

dE^ 2 



dt 



o-TC7^ y (1 - /3 cos ipfen{e)d€ (5.30) 



If the incoming photons are isotropically distributed, we can average out 
(1 — Pcosip)"^ over the solid angle, obtaining 1 + /3^/3. The power produced 
is then 

"^^^ - a^n' (l + ^]u, (5.31) 



dt ^ ' V 3 , 
where 

Ur = j en{e)de (5.32) 

is the energy density of the radiation before scattering. This is the power 
contained in the scattered radiation. To calculate the energy loss rate of the 
electron, we have to subtract the initial power of the radiation eventually 
scattered 

Pc(7) = ^ = ^ - <^TcC/r = -aWP'U, (5.33) 
dt dt 6 

A simple way to remember Eq. 15.331 is: 

n / N f # oi collisions \ , , r \ 

i-c(7j = I 1 (average phot, energy after scatt.) 

U, \ /'4„ , 2 



Note the similarity with the synchrotron energy loss. The two energy loss 
rates are identical, once the radiation energy density is replaced by the 
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magnetic energy density Ub- Therefore, if relativistic electrons are in a 
region with some radiation and magnetic energy densities, they will emit by 
both the synchrotron and the Inverse Compton scattering processes. The 
ratio of the two luminosities will be 



1^ ~ ^ ~ ITr 



(5.35) 



where we have set dEic/dt = dEg/dt. This is true unless one of the two 
processes is inhibited for some reason. For instance: 

• At (relatively) low energies, electrons could emit and absorb synchrotron 
radiation, so the synchrotron cooling is compensated by the heating 
due to the absorption process. 

• At high energies, electrons could scatter in the Klein-Nishina regime: 
in this case, since the cross section is smaller, they will do less scat- 
terings, and cool less. 



A The bowl 




Figure 5.7: In the center of a semi-sphere (the "bowl") we have relativistic 
electrons going down and going up, all with the same 7. Since the seed 
photon distribution is anisotropic, so is the scattered radiation and power. 
The losses of the electron going down are 7 times larger than those of the 
electron going up (if 7 3> 1). Since almost all the radiation is produced 
along the velocity vector of the electrons, also the downward radiation is 7 
times more powerful than the upward radiation. 

But let us go back to Eq. 15.301 that is the starting point when dealing 
with anisotropic seed photon distributions. Think for instance to an accre- 
tion disk as the producer of the seed photons for scattering, and some cloud 
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of relativistic electrons above the disk. If the cloud is not that distant, and 
it is small with respect to the disk size, then this case is completely equal to 
the case of having a little cloud of relativistic electrons located at the center 
of a semi-sphere. That is, we have the "bowl" case illustrated in Fig 15.71 
Just for fun, let us calculate the total power emitted by an electron going 
"up" and by its brother (i.e. it has the same 7) going down. Using Eq. I5.3UI 
we have: 

^up £(l-/3/.)2d^ 1-/3 + /3V3 ' 

where fi = cosip and the last step assumes /3 — t- 1. Since almost all the 
radiation is produced along the velocity vector of the electrons, also the 
downward radiation is more powerful than the upward radiation (i.e. 7 
times more powerful for 7^1). What happens if the cloud of electrons 
is located at some height above the bowl? Will the -Pdown/-fup be more or 
less? 

5.5.3 Cooling time and compactness 

The cooling time due to the inverse Compton process is 

E S'ymeC^ SnieC^ 2 /,- o^T^ 

This equation offers the opportunity to introduce an important quantity, 
namely the compactness of an astrophysical source, that is essentially the 
luminosity L over the size R ratio. Consider in fact how and L are related: 

Although this relation is almost universally used, there are subtleties. It 
is surely valid if we measured C/r outside the source, at a distance R from 
its center. In this case AirR'^c is simply the volume of the shell crossed by 
the source radiation in one second. But if we are inside an homogeneous, 
spherical transparent source, a better way to calculate 11^ is to think to 
the average time needed to the typical photon to exit the source. This is 
tesc = 3i?/(4c). It is less than R/c because the typical photon is not born 
at the center (there is more volume close to the surface). If V = {4tt/3)R^ 
is the volume, we can write: 

jj _ L_ _ 3L 3R _ 9L 

This is greater than Eq. I5.38l bv a factor 9/4. Anyway, let us be conventional 
and insert Eq. 15.381 in Eq. 15.371 

3-KmeC^R^ tiQ 37r rrieC^R Sir 1 

tic = J > = ^ = — 7 (5.40) 

aT^L R/c 7 (TT-t^ 7 t 
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where the dimensionless compactness £ is defined as 

i = ^ (5.41) 

For i close or larger than unity, we have that even low energy electrons cool 
by the Inverse Compton process in less than a light crossing time R/c. 

There is another reason why £ is important, related to the fact that 
it directly measures the optical depth (hence the probability to occur) of 
the photon-photon collisions that lead to the creation of electron-positron 
pairs. The compactness is one of the most important physical parameters 
when studying high energy compact sources (X-ray binaries, AGNs and 
Gamma Ray Bursts). 

5.5.4 Single particle spectrum 

As we did for the synchrotron spectrum, we will not repeat the exact deriva- 
tion of the single particle spectrum, but we try to explain why the typical 
frequency of the scattered photon is a factor 7^ larger than the frequency of 
the incoming photon. Here are the steps to consider: 

1. Assume that the relativistic electron travels in a region where there 

is a radiation energy density Ui- made by photons which we will take, 
for simplicity, monochromatic, therefore all having a dimensionless 
frequency x = hv/nicCp. 

2. In the frame where the electron is at rest, half of the photons appear 
to come from the front, within an angle I/7. 

3. The typical frequency of these photons is x' ~ 7a; (it is twice that for 
photons coming exactly head on). 

4. Assuming that we are in the Thomson regime means that 1) x' < 1; ii) 
the cross section is the Thomson one; iii) the frequency of the scattered 
photon is the same of the incoming one, i.e. x'l = x' ~ 7X, and iv) 
the pattern of the scattered photons follows the angular dependence 
of the cross section, therefore the "peanut" . 

5. Independently of the initial photon direction, and therefore indepen- 
dently of the frequencies seen by the electrons, all photons after scat- 
terings are isotropized. This means that all observers (at any angle 

in this frame sec the same spectrum, and the same typical frequency. 
Half of the photons are in the semi-sphere with ^[ < Tr/2. 

6. Now we go back to the lab-frame. Those photons that had 'ip'i< 11/2 
now have ■i/'i < 1/7- Their typical frequency if another factor 7 greater 
than what they had in the rest frame, therefore 

xi ^ 7^x (5.42) 
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This is the typical Inverse Compton frequency. 

The exact derivation can be found e.g. in Rybicki &; Lightman (1979) 
and in Blumenthal &: Gould (1970). We report here the final result, valid 
for a monochromatic and isotropic seed photons distribution, characterized 
by a specific intensity 

= — d[x — xq) (5.43) 

X X 

Note that I{x)/x is the analog of the normal intensity, but it is associated 
with the number of photons. If we have n electrons per cubic centimeter we 
have: 

The function Fiq contains all the frequency dependence: 

Xl 1 



Fic{xi) 
Fic{xi) 



Xq 
Xl 
Xo 



xo (1 + /3)272 
1 _ ^ 1 



xo (1 + /3)V 



1 Xl 

(1 + /3)V 
1< — < (1 + /3)V (5.45) 

Xo 



The first line corresponds to downs cattering: the scattered photon has less 
energy than the incoming one. Note that in this case -Fic(a^i) oc xf. The 
second line corresponds to upscattering: in this case Fic(xi) oc xi except 
for frequencies close to the maximum ones. The function -Fic(xi) is shown 
in Fig. 15.81 for different values of 7. The figure shows also the spectrum of 
the photons contained in the beaming cone I/7: the corresponding power is 
always 75% of the total. 

The average frequency of Fic(xi) is 

(xi) = 27^x0; energy spectrum (5.46) 

This is the average frequency of the energy spectrum. We sometimes want 
to know the average energy of the photons, i.e. we have to calculate the 
average frequency of the photon spectrum Fic(xi)/xi. This is: 

(xi) = -7^x0; photon spectrum (5-47) 

5.6 Emission from many electrons 

We have seen that the emission spectrum from a single particle is peaked, 
and the typical frequency is boosted by a factor 7^. This is equal to the 
synchrotron case. Therefore we can derive the Inverse Compton emissiv- 
ity as we did for the synchrotron one. Again, assume a power-law energy 
distribution for the relativistic electrons: 

N{^) = Kj-P = N{E) —■ 7^in < 7 < 7max (5.48) 

a7 
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Figure 5.8: Spectrum emitted by the Inverse Compton process by electrons 
of different 7 (as labeled) scattering an isotropic monochromatic radiation 
field of dimensionless frequency xq. The dashed line corresponds to the spec- 
trum emitted within the I/7 beming cone: it always contains the 75% of the 
total power, for any 7. For xi < xq we have downscattering, i.e. the photons 
loose energy in the process. Note also the power law segments arising when 
7 3> 1: Fic{xi) oc xf for the downscattering tail, and Fjc{xi) oc xi for the 
upscattering segment. 




and assume that it describes an isotropic distribution of electrons. For 
simplicity, let us assume that the seed photons are isotropic and monochro- 
matic, with frequency z^o (we now pass to real frequencies, since we are 
getting closer to the real world..). Since there is a strong link between the 
scattered frequency i^c and the electron energy that produced it, we can set: 



:7 ^0 



o \ 1/2 





-1/2 

_ f c / 


dv 


2 \ 



3 



1/2 



(5.49) 



Now, repeating the argument we used for synchrotron emission, we can state 
that the power lost by the electron of energy jmcC^ within nieC^d'y goes into 
the radiation of frequency u within dv. Since we will derive an emissivity 
(i.e. erg cm~'^ sterad^^ Hz~^) we must remember the Air term (if the 
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emission is isotropic). We can set: 



€,{i^,)di^, = -^mec'Pci^)N{j)d-f (5.50) 
47r 



This leads to: 



Again, a power law, as in the case of synchrotron emission by a power law 
energy distribution. Again the same link between a and p: 

a = ^ (5.52) 

Of course, this is not a coincidence: it is because both the Inverse Compton 
and the synchrotron single electron spectra are peaked at a typical frequency 
that is a factor 7^ greater than the starting one. 
Eq. 15.511 becomes a little more clear if 

• we express ec(fc) as a function of the photon energy hi'c- Therefore 
edhuc) = ec(fc)//i; 

• we multiply and divide by the source radius R; 

• we consider a proxy for the scattering optical depth of the relativistic 
electrons setting Tc = a^KR. 

Then we obtain: 

,,(ft.,) = _LM):_^itf;iy" (5.53) 

^ ^ 47r 2 R/c huo \uoJ ^ ' 

In this way: Tc (for < 1) is the fraction of the seed photons U^/huQ 
undergoing scattering in a time R/c, and z/c/fo ~ 7^ is the average gain in 
energy of the scattered photons. 



5.6.1 Non monochromatic seed photons 

It is time to consider the more realistic case in which the seed photons are 
not monochromatic, but are distributed in frequency. This means that we 
have to integrate Eq. l5.51l over the incoming photon frequencies. For clarity, 
let us drop the subscript in i/q- We have 

1 (4/3)" Tc r--u,{u) 

eM = — -^T^c / ^ du (5.54) 

Att 2 R/c J^^^^^ V 

where U-ciy) [erg cm~^ Hz^^] is the specific radiation energy density at 
the frequency v. The only difficulty of this integral is to find the correct 
limit of the integration, that, in general, depend on v^- Note also another 
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interesting thing. We have just derived that if the same electron population 
produces Inverse Compton and synchrotron emission, than the slopes of the 
two spectra are the same. Therefore, when Ur{i^) is made by synchrotron 
photons, then Ur{iy) oc The result of the integral, in this case, will be 

In ( Vraax / ^min ) • 



O 




Figure 5.9: The i^-fc plane. The two diagonal lines delimit the regions of 
the seed photons that can be used to give a given frequency fc- 



Fig. 15.91 helps to understand what are the right h'max and to use. 
On the y-axis we have the frequencies of the seed photon distribution, which 
extend between vi and U2. On the x-axis we have the scattered frequencies, 
which extend between fc,i = (4/3)7 ^j^ z^i and fc,4 = (4/3)7^ax'^i- The 
diagonal lines are the functions 



472. 

'mi 



^ I mi 



(5.55) 



that tell us what are the appropriate v that can give i/c once we change 7. 
There are three zones: 

1. In zone (1), between fc,i and fc,2 = (4/3) 7^^;^ 1^2 the appropriate limits 
of integration are: 

3i/c 



472. 

'mm 



(5.56) 
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2. In zone (2), between i'c,2 and 



(4/3)7^3x^1 the limits are: 



(5.57) 



3. In zone (3), between v^s and 



(4/3)7^3x^2 the hmits are: 




I^max — 1^2 



(5.58) 



We see that only in zone (2) the limits of integration coincide with the exten- 
sion in frequency of the seed photon distribution, and are therefore constant. 
Therefore ec(fc) will be a power law of slope a only in the corresponding 
frequency limits. Note also that for a broad range in [i'i;i'2] or a narrow 
range in [7min;7max] we do not have a power law, since there is no Vc for 
which the limits of integrations are both constants. 

5.7 Thermal Comptonization 

With this term we mean the process of multiple scattering of a photon due 
to a thermal or quasi-thermal distribution of electrons. By quasi-thermal 
we mean a particle distribution that is peaked, even if it is not a perfect 
Maxwellian. Since the resulting spectrum, by definition, is due to the su- 
perposition of many spectra, each corresponding to a single scattering, the 
details of the particle distribution will be lost in the final spectrum, as long 
as the distribution is peaked. The "bible" for an extensive discussion about 
this process is Pozdnyakov, Sobol &: Sunyaev (1983). 

There is one fundamental parameter measuring the importance of the 
Inverse Compton process in general, and of multiple scatterings in partic- 
ular: the Comptonization parameter, usually denoted with the letter y. Its 
definition is: 

y = [average # of scatt.] x [average fractional energy gain for scatt.] 



If y > 1 the Comptonization process is important, because the Comptonized 
spectrum has more energy than the spectrum of the seed photons. 

5.7.1 Average number of scatterings 

This can be calculated thinking that the photon, before leaving the source, 
experience a sort of random walk inside the source. Let us call 



(5.59) 



tt = CF^nR 



(5.60) 



the Thomson scattering optical depth, where n is the electron density and 
R the size of the source. When tt < 1 most of the photons leave the source 
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directly, without any scattering. When tt > 1 then the mean free path is 
d = R/tt and the photon will experience, on average, scatterings before 
leaving the source. Therefore the total path travel by the photon, from the 
time of its birth to the time it leaves the source is: the photon is born, is 



cAt = r| — = ttR (5.61) 

TT 

and At is the corresponding elapsed time. 

5.7.2 Average gain per scattering 
Relativistic case 

If the scattering electrons are relativistic, we have already seen that the 
photon energy is amplified by the factor (4/3)7^ (on average). Therefore 
the problem is to find what is (7^) in the case of a relativistic Maxwellian, 
that has the form 

iV(7) oc 7^6-^/®; e = — ^ (5.62) 

Setting Xo = hvo/^nieC^) we have that the average energy of the photon of 
initial frequency xq after a single scattering with electron belonging to this 
Maxwellian is: 

4 4 /-f7^e-VQrf7 

4 ^2r(5) 

= ^^o| = lee^xo (5.63) 

Non relativistic case 

In this case the average gain is proportional to the electron energy, not 
to its square. The derivation is not immediate, but we must use a trick. 
Also, we have to account that in any Maxwellian, but especially when the 
temperature is not large, there will be electrons that have less energy than 
the incoming photons. In this case it is the photon to give energy to the 
electron: correspondingly, the scattered photon will have less energy than 
the incoming one. Averaging out over a Maxwellian distribution, we will 
have: 

Ax Xi-XQ 

= = aB — X (5.64) 

X Xq 

Where aQ is what the photon gains and the —x term corresponds to the 
downscattering of the photon (i.e. direct Compton). We do not know yet 
the value for the constant a. To determine it we use the following argument. 
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We know (from general and robust arguments) what happens when photons 
and electrons are m equilibrium under the only process of scattering, and 
neglecting absorption (i.e. when the number of photon is conserved). What 
happens is that the photons follow the so-called Wien distribution given by: 

Fw(x) oc x^e-^/® ^ iVw(x) = oc x^e"^/® (5.65) 

where F correspond to the radiation spectrum, N to the photon spectrum, 
and Q is the dimensionless electron temperature. When a Wien distribution 
is established we must have (Ax) = 0, since we are at equilibrium. So we 
require that, on average, gains equal losses: 

{Ax) = ^ a@{x) - {x^) = (5.66) 

Calculating [x) and (x^) for a photon Wien distribution, we have: 

x3e-^/®da; r(4) 3! 

, .r^c-''''^dx Tib) o 4! , , . 

This implies that a = 4 not only at equilibrium, but always, and we finally 
have 

At 

— = 46 - X (5.68) 

X 

Combining the relativistic and the non relativistic cases, we have an expres- 
sion valid for all temperatures: 

A /J. 

= 166^ + 46 - a; (5.69) 

X 

going back to the y parameter we can write: 

y = max(TT, r|) x [166^ + 46 - x] (5.70) 
Going to the differential form, and neglecting downscattering, we have 

^ = [166^ + 49] dK ^xi = xo e(i6e^+40)^ ^ Xf = xo (5.71) 

where now K is the number of scatterings. If we subtract the initial photon 

energy, and consider that the above equation is valid for all the xq of the 
initial seed photon distribution, of luminosity Lq, we have 

^ = ey-1 (5.72) 

Then the importance of y is self evident, and also the fact that it marks the 
importance of the Comptonization process when it is larger than 1. 
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# scatt. Fraction of escaping (x) 



photons 

e^"^ — )• 1 — Tt xq 

1 ~ Tt xqA 

2 ~ r| xqA'^ 

3 ~ r| xqA^ 

4 ~ r| xoA^ 



Table 5.1: When tt < 1, a fraction e~'^'^ of the seed photons escape without 
doing any scattering, and a fraction 1 — e~'^ — )• r undergoes at least one 
scattering. We can then repeat these fractions for all scattering orders. Even 
if a tiny fraction of photons does several scatterings, they can carry a lot of 
energy. 

5.7.3 Comptonization spectra: basics 

We will illustrate why even a thermal (Maxwellian) distribution of electrons 
can produce a power law spectrum. The basic reason is that the total pro- 
duced spectrum is the superposition of many orders of Compton scattering 
spectra: when they are not too much separated in frequency (i.e. for not too 
large temperatures) the sum is a smooth power law. We can distinguish 4 
regimes, according to the values of tt and y. As usual, we set x = hv/{meC^) 
and = kT/{meC^). 

The case tt < 1 

Neglect downscattering for simplicity. The fractional energy gain is Ax/x = 
160^ + 40, so the amplification A of the photon frequency at each scattering 
is 

A = — = 1602 + 40 + 1 ~ ^ (5.73) 

X Tt 

We can then construct Table I^TTl 

A look to Fig. IS.lOl should convince you that the sum of all the scattering 
orders gives a power law, and should also make clear how to find the spectral 
slope. Remember that we are in a log-log plot, so the spectral index is 
simply Ay/ Ax. We can find it considering two successive scattering orders: 
the typical (logarithm of) frequency is enhanced by log^, and the fraction 
of photons doing the scattering is — logrx. Remember also that we use 
F{x) oc as the definition of energy spectral index. 

Therefore 

logTT logTT 
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XAv" A^v A*^v' A'^Y A^Y A®v" A'''v" 

Log X 



Figure 5.10: Multiple Compton scatterings when tt < 1- A fraction tt of 
the photons of the previous scattering order undergoes another scattering, 
and amplifying the frequency by the gain factor A, until the typical photon 
frequency equals the electron temperature G. Then further scatterings leave 
the photon frequency unchanged. 

When y ~ 1, its logarithm is close to zero, and we have a ~ 1. When 
y > 1, then a < 1 (i.e. flat, or hard), and vice-versa, when y < 1, then its 
logarithm will be negative, as the logarithm of tt, and a > 1 (i.e. steep, or 
soft). 

Attention! when tt ^ 1 and A is large (i.e. big frequency jumps be- 
tween one scattering and the next), then the superposition of all scattering 
orders (by the way, there are fewer, in this case) will not guarantee a perfect 
power-law. In the total spectrum we can see the "bumps" corresponding to 
individual scattering orders. 

The case tt>1 

This is the most difficult case, as we should solve a famous equation, the 
equation of Kompaneet. The result is still a power law, whose spectral index 
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is approximately given by 



a 




(5.75) 



The case tt ^ 1: saturation 

In this case the interaction between photons and matters is so intense that 
they go to equihbrium, and they wih have the same temperature. But 
instead of a black-body, the resulting photon spectrum has a Wien shape. 
This is because the photons are conserved (if other scattering processes 
such as induced Compton or two-photon scattering are important, then one 
recovers a black-body, because these processes do not conserve photons). 
The Wien spectrum has the slope: 



At low frequencies this is harder than a black-body. 
The case tt > 1, y > 1: quasi saturation 

Suppose that in a source characterized by a large tt the source of soft 
photons is spread throughout the source. In this case the photons produced 
close to the surface, in a skin of optical depth tt = 1, leave the source 
without doing any scattering (note that having the source of seed photons 
concentrated at the center is a different case). The remaining fraction, 
1 — 1/tt, i.e. almost all photons, remains inside. This can be said for each 
scattering order. This is illustrated in Fig. 15.111 where tt corresponds to 
the ratio between the flux of photon inside the source at a given frequency 
and the flux of photons that escape. If I start with 100 photons, only 1 - 
say - escape, and the other 99 remain inside, and do the first scattering. 
After it, only one escape, and the other 98 remain inside, and so on, until 
the typical photon and electron energies are equal, and the photon therefore 
stays around with the same final frequency until it is its turn to escape. 
This "accumulation" of photons at x ~ 30 gives the Wien bump. Note that 
since at any scattering order only a fixed number of photons escape, always 
the same, then the spectrum in this region will always have a = 0. This is 
a "saturated" index, i.e. one obtains always zero even when changing tt or 
Q. What indeed changes, by increasing tt, is that i) the fiux characterized 
by decreases, ii) the Wien peak will start to dominate earlier (at lower 
frequencies), while nothing happens to the fiux of the Wien peak (it stays 
there). Increasing still tt we fall in the previous case (equilibrium, meaning 
only the Wien spectrum without the part). 




(5.76) 
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Figure 5.11: Multiple Compton scatterings when tt > 1 and y S> 1. For the 
first scattering orders, nearly all photons are scattered: only a fraction I/tt 
can escape. Therefore the number of photons escaping at each scattering 
order is the same. This is the reason of the flat part, where F{x) oc x^. 
When the photon frequency is of the order of 0, photons and electrons are 
in equilibrium, and even if only a small fraction of photons can escape at 
each scattering order, they do not change frequency any longer, and therefore 
they form the Wien bump, with the slope F{x) oc x^e~^/®. If we increase 
Tt, the flux with slope x^ decreases, while the Wien bump stays the same. 
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Chapter 6 

Synchrotron Self— Compton 



Consider a population of relativistic electrons in a magnetized region. They 
will produce synchrotron radiation, and therefore they will fill the region 
with photons. These synchrotron photons will have some probability to 
interact again with the electrons, by the Inverse Compton process. Since the 
electron "work twice" (first making synchrotron radiation, then scattering 
it at higher energies) this particular kind of process is called synchrotron 
self-Compton, or SSC for short. 



6.1 SSC emissivity 

The importance of the scattering will of course be high if the densities of 
electrons and photons are large. If the electron distribution is a power 
law [A^(7) = 1^7"^], then we expect that the SSC flux will be oc K^, i.e. 
quadratic in the electron density. 

We should remember Eq. 15.541 and, instead of a generic Ur{i^), we should 
substitute the appropriate expression for the specific synchrotron radiation 
energy density. We will then set: 

where 3i?/(4c) is the average photon source-crossing time, and V is the 
volume of the source. Now a simple trick: we write the specific synchrotron 
emissivity as 

e,{u) = es,oi^~" (6.2) 

Remember: the a appearing here is the same index in Eq. 15.541 Substituting 
the above equations into Eq. 15.541 we have 

(4/3)"-i dz. 
essc(i^c) = 5 'rces,oz^c / — (6-3) 

^ ./]/ . 
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As you can see, es^o'^cT" = es(t'c) is nothing else than the specific synchrotron 
emissivity calculated at the (Compton) frequency Uq- Furthermore, the in- 
tegral gives a logarithmic term, that we will call In A. We finally have: 



(4/3)""^ 
M = rces(z^c)lnA 



(6.4) 



In this form the ratio between the synchrotron and the SSC flux is clear, it 
is [(4/3)'^~"'^/2]tc In A ~ TclnA. It is also clear that since Tc = a^RK and 
es(fc) oc KB^^", then, as we have guessed, the SSC emissivity essc(i^c) oc 
(i.e. electrons work twice). Fig. 16.11 summarizes the main results. 
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Figure 6.1: Typical example of SSC spectrum, shown in the vF^, vs v rep- 
resentation. The spectral indices instead correspond to the F^ oc z^~° con- 
vention. 
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6.2 Diagnostic 

If we are confident that the spectrum of a particular source is indeed given by 
the SSC process, then we can use our theory to estimate a number of physical 
parameters. We have already stated (see Eq. I4.32p that observations of 
the synchrotron spectrum in its self-absorbed part can yield the value of 
the magnetic field if we also know the angular radius of the source (if it 
is resolved). Observations in the thin part can then give us the product 
RK = Tc/dT (see Eq. I4.28p . But Tc is exactly what we need to predict 
the high energy flux produced by the SSC process. Note that if the source 
is resolved (i.e. we know Bg) we can get these information even without 
knowing the distance of the source. To summarize: 



There is an even simpler case, which for reasons outlined below, is the most 
common case employed when studying radio-loud AGNs. In fact, if you 
imagine to observe the source at the self absorption frequency ft, then you 
are both observing the thick and the thin flux at the same time. Then, let 
us call the flux at vt simply Ft. We can then re- write the equation above: 



Once again: on the basis of a few observations of only the synchrotron flux, 
we can calculate what should be the SSC flux at the frequency Vc- Note the 
rather strong dependencies, particularly for Og, in the sense that the more 
compact the source is, the larger the SSC flux. 

If it happens that we do observe the source at high frequencies, where 
we expect that the SSC flux dominates, then we can check if our model 
works. Does it? For the strongest radio-loud sources, almost never. The 
disagreement between the predicted and the observed flux is really severe, 
we are talking of several orders of magnitude. Then either we are completely 
wrong about the model, or we miss some fundamental ingredient. We go for 
the second option, since, after all, we do not find any mistake in our theory. 

The missing ingredient is relativistic bulk motion. If the source is mov- 
ing towards us at relativistic velocities, we observe an enhanced flux and 
blueshifted frequencies. Not accounting for it, our estimates of the magnetic 




(6.5) 





oc i.2(2+a)^-(5+3a)^-2(3+2«)^^- 



(6.6) 
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field and particles densities are wrong, in the sense that the B field is smaller 
than the real one, and the particle densities are much greater (for smaller B 
we need more particle to produce the same synchrotron flux). So we repeat 
the entire procedure, but this time assuming that F{i^) = (53+"F'(i/), where 
5 = l/[r(l — /3 cos 6)] is the Doppler factor and F^v) is the flux received by 
a comoving observer at the same frequency z^. Then 

^tMn(^) « elRKB^+''u"'5'^+'^ (6.7) 
The predicted SSC flux then becomes 

F,M OC F2(2+a)^-(5+3a)^^„2(3+2a)^-a^-2(2+a) (g_g) 

If we now compare the predicted with the observe SSC flux, we can estimate 
5. And indeed this is one of the most powerful 5-estimators, even if it is not 
the only one. 



6.3 Why it works 

We have insisted on the importance of observing the synchrotron flux both 
in the self-absorbed and in the thin regime, to get B and t^- But the 
self-absorbed part of the synchrotron spectrum, the one oc v^/"^ is very 
rarely observed in general, and never in radio-loud AGNs. So, where is the 
trick? It is the following. In radio-loud AGN the synchrotron emission, at 
radio frequencies, comes partly from the radio lobes (extended structures, 
hundreds of kpc in size, very relaxed, unbeamed, and usually self-absorbing 
at very small frequencies) and from the jet. The emission from the latter is 
beamed, and it is the superposition of the fluxes produced in several regions: 
the most compact ones (closer to the central engine) self-absorb at high radio 
frequencies (say, at 100 GHz), and the bigger they are, the smaller their self- 
absorbed frequency. But what is extraordinary about these jets is that the 
peak flux of each component (i.e. the flux at the self-absorption frequency) 
is approximately constant (in the past, this phenomenon was called cosmic 
conspiracy). Therefore, when we sum up all the components, we have a flat 
radio spectrum, as illustrated by Fig. 16.21 

Of course the emission components of the jet, to behave in such a coher- 
ent way, must have an electron density and a magnetic field that decrease 
with the distance from the central engine in an appropriate way. There is 
a family of solutions, but the most appealing is certainly B{R) oc R~^ and 
K{R) oc R^^. It is appealing because it corresponds to conservation of the 
total number of particles, conservation of the bulk power carried by them 
(if r does not change) and conservation of the Poynting fiux (i.e. the power 
carried by the magnetic field). 
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Figure 6.2: Typical example of the composite spectrum of a flat spectrum 
quasars (FSRQ) shown in the F^, vs representation, to better see the flat 
spectrum in the radio. The reason of the flat spectrum is that different parts 
of the jet contributes at different frequencies, but in a coherent way. The blue 
line is the SSC spectrum. Suppose to observe, with the VLBI, at 22 GHz: 
in this framework we will always observe the jet component peaking at this 
frequency. So you automatically observe at the self-absorption frequency of 
that component (for which you measure the angular size). 



To our aims, the fact that the jet has many radio emission sites self- 
absorbing at different frequency is of great help. In fact suppose to observe 
a jet with the VLBI, at one frequency, say 22 GHz. There is a great chance to 
observe the emission zone which is contributing the most to that frequency, 
i.e. the one which is self-absorbing at 22 GHz. At the same time you 
measure the size. Then, suppose to know the X-ray flux of the source. It will 
be the X-ray flux not only of that component you see with the VLBI, but an 
integrated flux from all the inner jet (with X-ray instruments the maximum 
angular resolution is about 1 arcsec, as in optical). But nevertheless you 
know that your radio blob cannot exceed the measured, total. X-ray flux. 
Therefore you can put a limit on 5 (including constants): 



1 



6 > (0.08a + 0.14) (1 + z) ( 1 



In 



2(2 + a) 



( 20, \--^ 



5 GHz/ V m.a.s 



Jy / VJy / ^1 keV/ 

3 + 2a 1 

)2(2+Q,) 
(6.9) 



For some sources you would find 5 > 10 or 20, i.e. rather large values. 
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Chapter 7 

Pairs 



In the presence of energetic particles and pliotons one has to wonder about 
the possibility that there are collisions between them. One result of these 
collisions is the production of e^. If we have a photon-photon collision, 
then the original photons might disappear, so that this process becomes an 
important absorption process. 

The importance of this process came initially from the realization that 
the virial temperature of protons, in the vicinity of black holes, can be very 
large. As an estimate, the kinetic energy of a proton at 3 Schwarzschild 
radii is 



therefore particles can be very energetic in accreting systems. Be aware, on 
the other hand, that protons can efficiently give their energy to electrons, 
that will emit this energy. This would keep the protons much cooler than 

the above vahic. 

The other key ingredient, realized in the '60s, is that Active Galactic Nu- 
clei (AGNs) vary quickly (with a variability timescale iyar)- By the causality 
argument they cannot be larger than R ~ ctvar- Therefore their emitting 
regions must be small, and yet produce a huge luminosities. Densities are 
thus very large, and the collisions between particles, between particles and 
photons, and between photons must be probable. The discovery that blazars 
(AGNs with jets pointing at us) are very efficient 7 ray producers is now 
a proof that high energy particles indeed exist in these objects, that can 
in turn produce pairs. And finally, also Gamma Ray Bursts (GRBs) are 
emitting 7-rays. 

The importance of electron-positron pair production in the nuclei of 
AGNs was realized quite early (Jelley 1966) and at the beginning of the 
'70s' two seminal papers appeared (one by Bysnovaty-Kogan, Zeldovich 
& Sunyaev 1971, the other by Bonometto & Rees 1971): the first con- 
cerned thermal plasmas, the second was about pairs in non-thermal plas- 




MeV 



(7.1) 
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mas. Thereafter the field was sleepy for about one decade, when the paper 
by Lightman (1982) and Svensson (1982) resurrected the interests of the 
scientific community. 

At those time it was believed that pair processes were fundamental for 
the formation of the spectrum of all AGNs, galactic binaries and GRBs. Now 
we know that they do not play such a key role for AGNs and binaries, but 
they still could for GRBs. Nevertheless the entire subject is worth studying, 
because, at the very least, it allows to pose robust and important limits on 
the physics of all compact objects. 

Wc will divide the study into thermal and non-thermal plasmas. Wc 
will see two fundamental results of pairs in thermal plasmas, and will see 
in more detail what happens when the particles have a relativistic (and 
not-Maxwellian) energy distribution. 



7.1 Thermal pairs 

Consider particles characterized by a Maxwellian distribution and a tem- 
perature T. Since we are dealing with the production of e^ , it is convenient 
to measure all energies in units of mgC^ (=511 keV). Let us define: 

kT hv 
e = 2, X = 2 7.2 

Furthermore, define: 

n± = number density of positrons (+) and electrons (— ) 

= number density of photons 

n± = production rate of pairs 

riA = rate of annihilation of pairs (7.3) 

To these processes we should add the possibility that the leptons escape the 
source, with a rate 

•hesc = rate of escape of leptons (7.4) 
The process that can produce e^ are: 
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while the processes that annihilate pairs or that correspond to the escape of 
pairs and electrons are: 

e^e^ — 77 : h\ = 2rgCn+n_FA 

Escape: hesc = r^c (n+ + n_)^ Fesc (7.6) 

The F-factors are averages of the energy dependent part of the cross sections 
over the Maxwellian distribution. They are dimensionless. The geometrical 
part of the cross section is always of the form of a^r"^: ap = 1/137 is the 
fine structure constant, and is the classical electron radius. For photon- 
photon interaction, a = 0, for particle-photon a = 1, for particle-particle 
a = 2. Consider that pair production processes have an energy threshold: 
there must be enough energy to produce a pair. But, when dealing with 
a Maxwellian distribution, one has particles (even if a few) at all energies, 
so the process can occur even if <C 1. Here is a brief summary of the 
-F-factors: 





e < 1 


e> 1 


Comment 


-^77 
7? 

^ ee 
7? 

^ ep 

Fa 

^ esc 


vr 

— B 


Trine 

2 /0^2a 

3 \ 2j 

#fln'(2e) 

ln^(2e) 

TT 

Trine 
262 

3tt ^esc 


Wien 

power law (a > —1) 



Table 7.1: Summary of the F-factors for the main rates of pair produc- 
tion mechanism. Only photon-photon and particle-particle processes are 
considered. 



When dealing with processes involving photons, we have to know their 
spectral distribution. This is why in Tab. 17. II there are the rates for a Wien 
distribution [L{x) oc exp(— x/B)], and for a power law [L{x) oc x~"]. 
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7.1.1 First important result 

We impose that the source is stationary. This imphes a detailed balance 
between the created and the destroyed pairs: 

n+ = UA (7.7) 

For simplicity, we neglected escape as a way to get rid of pairs (but in reality 
there will always be some escape). To solve the above equation, we have 
to know the i^-terms, that are functions of 0. Now the important point: 
for particle-particle pair production processes, the F~terms are increasing 
functions of the temperature. On the contrary, the pair annihilation rate 
decreases with the temperature. Compare this behavior with the "Klein- 
Nishina" decline of the scattering cross section. 

Therefore, if the temperature is high enough, there is the possibility that 
we produce more pairs than what we can destroy. The source then becomes 
non-steady. As a result, if we want a steady source, we require that the 
temperature cannot be greater than a critical value. This value is: 

e = 24 ^ /cTmax = 12 MeV (7.8) 

To see this, consider a plasma that is pair-dominated (i.e neglect the original 
electrons associated with protons) . Let us set = ?^+ + ?^- . Neglect also 
all the particle-photon processes. Stationarity demands: 



a| rig Fee + n^F-y-y = '^u^Fa 

2Fa - al F^e 

Tie / ^77 



> (7.9) 



Fig. 17.11 shows how the two functions -Fa and Op-Fee behave. We require 
that Fa > Op Fee, since otherwise the annihilation rate cannot keep the pace 
with the production rate, the number of particles increases, and if they are 
at the same temperature, then they continue to produce even more pairs... 
(and this would quickly make an explosion). 

But in reality what happens is that the energy stored in the system 
is limited: if the number of particles increases, then they will share the 
available energy among themselves. There will be less energy per particle — >• 
the temperature will decrease. This is the second important result, that we 
will see now in more detail. 



7.1.2 Second important result 

Usually, when we increase the energy of a system, the temperature increases. 
But when kT approaches nieC^, the pair production processes start to be 
important. New particles are created: the energy per particle decreases 
— )• kT decreases. 
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Figure 7.1: To illustrate pair balance: the annihilation rate can balance 
pair production at low temperatures, but as = kT/nieC^ becomes large, 
the particle-particle processes increase their rate. Above Tmaxj annihilation 
cannot keep up with pair production, pair balance is impossible, and the 
number of pairs increases. Having more particles to share the available 
energy it is very likely that the temperature then decreases, making the 
system to return to pair balance. 



Example: Wien plasmas 

The name "Wien plasma" indicates a plasma in which the photons and 
particles are in equilibrium, reached through scattering. We must have tt = 
aicricR > 1. The spectrum of the radiation is not a black body, because 
we do not have absorption and re-emission, but only scattering. In these 
conditions, whatever the process of photon production (i.e. bremsstrahlung) 
the spectrum of the radiation is 

n^(e) oc exp(— e/ZcT) — )■ n^(x) oc exp(— x/G) (7.10) 

where e = hiy. n-y{e) is the photon production rate. To convert it in a density 
of photons, we use a very simplified equation of radiative transport: 

R 

n^ix) = h^ix) — (1 + Tt) (7.11) 
c 

where R is the size of the region. 
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Pause and observe the above equation: when tt ^ 1, the difference 
between and is the light crossing of the photon inside the source. It is 
the time needed for a photon created in the center to escape. When tt ^ 1, 
we have that the "transit time" is longer by a factor tt. This corresponds 
to have done Tr^ scatterings (remember the random walk) before escaping. 
The distance travelled between two consecutive scatterings is the mean free 
path, namely R/tt- Hence the time needed to escape is Tr^R/ (ttc) = ttR/c. 
With this link between the production rate and the density of photons, we 
can write the luminosity 

L = ^vriJ^y eh^{e)de 

= ^[e^de (7.12) 

3R/cJ l + TT ^ ^ 

When Tt ^» 1, we can set tt + 1 ~ tt = fTT-R?^e• In this case 

L = ^{e)^ (7.13) 
3(Tt rie 

where (e) is the mean energy of the photons. Note these remarkable facts: 

1. L depends on the ratio n^/nc. For pair plasmas (i.e. when pairs 
outnumber original electrons), this ratio is fixed by the F terms. 

2. For a Wien spectrum, (e) = 3kT (and then (x) = 30). 

3. For large luminosities, the dominating pair production process is photon- 
photon. 

4. Therefore n^/uc is fixed by (Fa/F^^)^/^, p^om Tab. EH assummg 
<C 1, we have 



Ue \F^^J Vvr2e-2/e 

n. 
n, 



We can see how the luminosity of a pair dominated, steady Wien source 
depends on the temperature: 

L = ^TT— 3 mec^ 2^27^ 03/2 e^/® 
3 (Tt 

oc ©5/2 gi/e (7 

For a given R, high luminosities are possible only for low temperatures! In- 
troducing the compactness 

i = ^ (7.16) 
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Figure 7.2: Wien plasmas: assuming that the pairs outnumber original 
electrons, there is a fixed relation between the compactness and the tem- 
perature. Large compactnesses (and then luminosities, for a given size), are 
possible only if the temperature is small. 

the luminosity-temperature relation has a simple form: 

^Wien = lOvrv^ gl/0 ^ (7^17) 

Fig. 17.21 shows ^wien as a function of 0. Note the minimum for ~ 1. The 
dashed line corresponds to the region where the asymptotic approximations 
are invalid. 



7.2 Non— thermal pairs 

In non-thermal plasmas the density of the particles is small, and the particle 
energies can be large. If so, the production rate is usually large, correspond- 
ing to a large density of very energetic photons even if the plasma is rarefied. 
In these conditions the main pair production process is photon-photon. 
The threshold energy for pair production is (x = hu/nieC^): 

xiX2 > ^ ^ , (7.18) 
1 — cos 

where 9 is the angle between the two photons: 6 = 27r means head on 
collisions. In this case the threshold is xiX2 > 1- For tail on encounters 
(9 = 0) the process cannot occur. 
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Figure 7.3: The photon-photon cross section, as a function of xia;2(l— cos ip), 
where ^p is the angle between the direction of the two photons (i.e. ^p = 180° 
corresponds to head on colhsions). The cross section is in units of the 
Thomson cross section ut. 



Naively, one would have thought that the threshold condition were xi + 
X2 > 2, namely that the sum of the photon energies is larger than 2meC^. 
But the threshold condition must be calculated in a reference frame where the 
energies of the two photons are equal. To illustrate this point, consider the 
simple case of an head on collision, with xi ^ X2- Let us go in a frame where 
x'l = X2- To do that, let us go in a reference frame having a Lorentz factor 
r such that the two energies become equal. In this frame the angle would 
still be 27r (head on collisions), and the two energies will be x'l = xi/{2T) 
and X2 = 2Tx2- In this frame we require that 



x'l + x'2 > 2 



and 



2r=(^ 

X2 



1/2 



^ + 2rX2 > 2 XiX2 > 1 



(7.19) 



The cross section resembles the Klein-Nishina cross section, and it has a 
peak for xiX2 ^ 2 where its value is 



a. 



77 



{xiX2 ~ 2) 



(7.20) 



The cross section is of course zero below threshold, and decreases after the 
peak. This means that a 7-ray photon with energy x can interact with all 
photons with energy xt > I/2;, but it will interact preferentially with those 
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photons near its threshold, i.e. xt = l/x (here the subscript "T" stands for 
"target"). 

7.2.1 Optical depth and compactness 

We want to calculate the optical depth t^^ for a source of a given luminosity. 
Assume that this luminosity is at high frequencies. A photon with energy 
~0.5-l MeV will interact with photons of the same energy (i.e. 0.5-1 MeV) 
to produce pairs (photons of 5 MeV will have x = 10 and will interact mainly 
with photons of energy xt ~ 1/10, i.e. 50 keV. 

Like all the optical depths, t^^ depends on the cross section, the density 
of the targets, and the size for which the process can occur. The density of 
targets is 

n (IMeV) ~ ^rad(lMeV)xlMeV ^ 

^ r^^(l MeV) c^^R ^ „ = ^ (7.21) 

" 5 AirR^meC^ 207r RrrieC^ 

Therefore T^-y depends (apart from the numerical factor) upon the so-called 
compactness parameter i, defined as: 

£ = — ^ (7.22) 

and we have: 

r^^(lMeV) ^ ^ (7.23) 
If we measure L in Eddington units, and R in Schwarzschild radii, we have 

^ ^ 27r mp L Rs 

3 TTte i^Edd R 

imd = Cax ^ 10^ (7.24) 



7.2.2 Photon-photon absorption 

Let us assume that a compact source produces a spectrum with a power law 
shape: F(x) oc a;~", where a is the energy spectral index. Assume also that 
this spectrum extends above the threshold for pair production. We want to 
calculate the optical depth for the 7^ — >■ e"^ process. 

Consider a 7-ray photon with x > 1. First we need to know the number 
density of the targets. These are all photons of energy xt > 1/x. On the 
other hand, we know that the cross section decreases when increasing the 
energy of the targets. In most cases, furthermore, we deal with spectral 
indices a > 0, meaning that also the number of photons decreases with 
energy (but be aware that this "rule" has an important exception, since 
several Gamma Ray Bursts have very flat spectra). 
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Figure 7.4: Photons of energy x > 1 can interact with photons of energy 
XT > 1/x to produce pairs. Since the cross section is decreasing when 
increasing the energy of the targets, and since for a power law spectrum 
with a > the number of photons is also decreasing when increasing their 
energy, we can approximate the number of target photons as (1/x) • n(l/x). 
Note the "x-l/x" symmetry when we plot F(x) in a logarithmic scale. 



We then adopt the following important simplification: 



T^^ = a^^R / n{x)dx 

Jl/x 

~ — R — n — 

5 X \x 

1 F(l/x) 1 /I 
oc oc — — 

XX X \x 

oc x° 



-a-1 



(7.25) 



Therefore r^^ increases increasing x: there are more targets (if a > 0). 
Consider also that the spectral index a is the spectral index of the targets: 
the same calculation would hold if we had a broken power law, with the 
break at -say- x = 1. Let us calculate the resulting spectrum, recalling the 
simplest equation of radiative transfer through a slab (a slab, instead of a 
sphere, is here used for simplicity): 



/fx) 



e{x)R 
k{x)R 



1 



(7.26) 



where /(x) is the monochromatic intensity; e(x) is the emissivity and k(x) 
is the absorption coefficient. The little trick that we have used here is to 
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multiply both the numerator and the denominator by R, the size of the 
source, because in this way k{x)R = r^^(x), and this allows to write: 



I{x) = e{x)R 



1 



(7.27) 



Written in this way, it is easy to see the behavior when T^-y(x) is much 
smaller or much larger than unity: 



T^7(x) ^ 1 
T^^{x) » 1 



I{x) = e{x)R oc x"" 
I{x) 



T^^[X) 



(7.28) 



There is a well defined break in the spectrum, occurring at the energy x for 
which Tj-y{x) ~ 1... 

BUT 

this is not the end of the story. 




-2-10 1 2 

Log X 

Figure 7.5: The expected spectrum after photon-photon absorption, but 
without reprocessing. The initial spectrum is a single power law of index a. 
Absorption produces a break at the energy x, where the optical depth r-y-y 
is unity. After the break, the spectrum becomes oc x~'^°'. 

In fact the produced pairs emit, and contribute to the spectrum (if they 
remain into the source). In general, each pair of electron and positron 
will share the energy of the photon creating them: (i.e. the electron and 
positron will have ^nieC^ ~ hv/2, i.e. 7 = x/2). Therefore the radiation 
they can produce will be at energies smaller than x. Depending on the 
energies, the secondary photons can themselves pair-produce. Alternatively, 
the secondary photons can be used as targets, enhancing the rate of the 
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process, and the number of pairs that the source can produce. Therefore 
the pair production process is 

highly non hnear 
photons that produce pairs that produce photons that produce pairs... 



7.2.3 An illustrative case: saturated pair production 

If the source is compact, meaning that the compactness parameter £ > 
1, the photon density is huge. This means that the cooling timescale for 
inverse Compton scattering is very short, shorter than the hght crossing 
time [remember that i measures not only the importance of the 7-7 ^ 
process, but also the ratio tcooi/(-R/c), see Eq. 4.38]. Therefore, since all 
leptons cool in a short time, to have a steady source we have to replenish 
the "dead leptons" with new ones. In general, we have two ways to do that. 

The first is to have "heating". With this word we mean that, besides 
cooling, leptons are also heated by some mechanism. Then there is a com- 
petition between energy losses and energy gains. Since the cooling rate for 
inverse Compton is oc 7^, it is very likely that cooling "wins" at high en- 
ergies, and heating at low energies. There will be a specific energy where 
heating and cooling balance. All leptons will quickly move (in the energy 
space) towards this energy. In the absence of diffusion (in the energy space) 
the leptons will form a monoenergetic distribution. Diffusion will broaden 
this, and we may have the formation of a Maxwellian particle distribution 
even if the particles do not directly exchange energy between them. 

The second way to "refill" the source with fresh energy is to continuously 
inject new energetic particles throughout the source. A shock does exactly 
that: at the shock front the particles are accelerated in a time that must be 
shorter than the cooling time. Particles are therefore continuously injected 
downstream, where they cool. Proton-proton collisions is another example 
(but also photon-photon collisions producing pairs can do). There can be 
differences among these processes: for instance, we expect a stratification in 
the downstream region of a shock, with the most energetic particles closer to 
the shock front, while we might expect a more homogeneous distribution of 
particles following the injection through proton-proton or photon-photon. 
We will neglect these differences in the following, and assume, for simplicity, 
that the injection and cooling mechanisms are homogeneous. 

Let us assume that we inject Q{'j) leptons per second and per cm^: 



number of particles 



cm^ s 



injection (7.29) 



This is not the quantity we need to calculate the spectrum. We want the 
density: 

number of particles 



Nil) 



cm-^ 



density (7.30) 
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The injection rate and the density are hnked by the continuity equation. In 
its simplest form it reads: 

^IM. = !_ [^iv(7)] + Q(7) (7.31) 

It describes how A^(7) evolves in time. If we want a steady source, we require 
dN{'y)/dt = 0. This allows to obtain the solution: 

Nh) = '-^ (7.32) 

7 

To go further let us assume that: 

1) the injection rate is a power law in energy: (5(7) = QoT"''; 

2) the cooling rate is quadratic with energy: 7 oc 7^; 

3) all particles cool down to 7 = 1, and then "disappear" (i.e. we neglect 
the accumulation of "dead leptons"). 

We have 2 cases: 

1) s < 1 ^ iV(7) oc 7"^ 

2) s > 1 ^ iV(7) oc 7-(^+i) (7.33) 

If we inject a monoenergetic (^(7) = Qo — 7max) we have A'' (7) oc 7"^. 

We have learnt that if N{'^) is a power law, the emitted spectrum will 
also be a power law, with a well defined link between the slope p of N{'^) 
and the radiation spectrum a. We have a = {p — l)/2. 

Now we have all the ingredients to see the effects of the reprocessing of 
pairs on the emitted spectrum. To fix the ideas, consider a region of radius 
R close to an accretion disk emitting UV radiation. Assume to continuously 
inject, throughout this region, relativistic leptons. Just to make it simpler, 
consider a monoenergetic injection at some 7 = 7max- Assume that the 
inverse Compton process is the dominant one. Let us see what happens in 
a few steps: 

1) After cooling, the particle density iV(7) oc 7"^, and the inverse Compton 
spectrum will have L{x) oc x^^'^ . Assume that the maximum photon energy 
(xmax) is well above the threshold for pair production, and that the com- 
pactness is very large for all photons with a; > 1. In these condition, all high 
energy photons get absorbed and produce pairs {saturated pair cascade). 

2) The generation of pairs means that, besides (5(7), there is another injec- 
tion term. Call it -Pi (7). We know the slope of Pi (7). In fact it will be the 
same of the photon spectrum, whose energy index is a = 0.5, and therefore 
its photon index is a -|- 1 = 1.5. 

Lix) 

Piij) = 2n^(27) oc x>l 

Pi(7) oc 7~""' = 7"'-' (7.34) 
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3) This new injection term originates a corresponding density of the first 
generation of pairs: 



iVi(7) oc 



r-Pi(y)d7' 



r 



oc 7 



-2.5 



(7.35) 



And consequently the radiation spectrum due to this first generation of pairs 
will have a spectral index 



2.5 - 1 



0.75 



(7.36) 




Figure 7.6: Reprocessing due to pairs. A source produces a soft photon dis- 
tribution (here it is a black body, blue solid line). The relativistic primary 
leptons injected throughout the source produce, by inverse Compton, a pri- 
mary high energy spectrum oc x~^'^ (blue short-dashed line). All photons 
above threshold get absorbed and transformed in pairs. They emit, and 
the spectrum of this first generation of pairs is oc x~^''^^ (green long-dashed 
line). Again, the photons above thresold are absorbed and create the second 
generation of pairs, whose spectrum oc x~^'^'^^ (dotted red line). The process 
can continue, and the limiting spectrum has a spectral index close to unity. 
The thick black solid line shows how the final spectrum would look like: it 
is a broken power law with indices ~1 and ~2 below and above x ~ 1. 



4) The cycle repeats. If pairs have enough energy to emit above threshold, 



7.2. NON-THERMAL PAIRS 



121 



their photons will get absorbed, and generate the second generation of pairs. 



P2(7) oc 7""^"' = r'-'' 
A'2(7) oc — 2 ^ 7 



7 _ ' ' ' „,-2.75 

7^ 

2.75 - 1 



a2 = = 0.875 (7.37) 



5) If these second-generation pairs have enough energy to produce photons 
above threshold, we will have a third generation of pairs, and so on. A 
very interesting result of this exercise is that the final spectral index tends 
towards unity. 

= y + ^ (7-38) 

where i stands for the i-th generation of pairs. The equation above is 
valid even if we start with a steep spectrum (q > 1). But in this case the 
reprocessed energy will be a minor fraction of the total: most of the power 
is emitted at low frequencies. So, even if we correctly derive the slope of the 
reprocessed spectrum, this will be energetically unimportant. 

If the initial spectrum is flat (i.e. a < 1), instead, the reprocessing due 
to pairs can be very important, and will have the effect to soften the overall 
spectrum, redistributing photons at lower frequencies. 



7.2.4 Importance in astrophysics 

The pair production processes were intensively studied in the '80s, especially 
in the field of Active Galactic Nuclei (AGNs). Researchers, at that time, 
hoped that the reprocessing due to pairs could explain the "universality" of 
the X~ray spectrum of AGNs, that has ax ~ 0.7. This would have implied a 
non-thermal origin of the spectrum, required to extend much above thresh- 
old. In other words, it was though that the typical AGN spectrum would 
resemble the one shown in Fig. 17.61 It was then predicted that all AGNs, 
i.e. also the radio-quiet ones, emit in the 7-rays, with a relatively steep 
power law shape. The first results obtained at high energies, coming from 
the Compton Gamma Ray Observatory, CGRO satellite showed instead an 
exponential cutoff at a few hundreds keV, and no signs of 7-rays in radio- 
quiet AGNs. It was then concluded that 1) the X-ray spectra of radio-quiet 
AGNs have a thermal origin (i.e. Comptonization) and 2) that pairs have 
no or little role in the formation of this spectrum. 

On the other hand, the main surprise of CGRO was the discovery that 
radio-loud AGNs are strong 7-ray emitters, as a class (i.e. we believe 
that all of them do). The properties of the observed spectra tell us that, 
again, photon-photon absorption inside the source rarely occurs (we do see 
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a lot of 7-rays...). The absence of this effect poses a problem, because the 
high energy flux of these sources is varying quickly, requiring very compact 
sources, and the observed large luminosities imply very large values of the 
compactness parameter. The process should indeed occur, and yet it does 
not. The solution of this puzzle is relativistic motion, shortening variability 
timescales, blue-shifting the observed frequencies, and enhancing the ap- 
parent luminosities. Indeed, requiring that the sources are optically thin for 
photon-photon absorption, we can derive a limit on the bulk Lorentz factor 
of the emitting sources. 

Another interesting field where pairs can be important is Gamma Ray 
Bursts (GRBs). Again, these sources emit high energy radiation that varies 
quickly. Since we now know that they are cosmological, the implied powers 
are huge, as well as the estimated photon densities. As in radio-loud AGNs, 
we can estimate a limit to their bulk Lorentz factor requiring that they are 
optically thin to the photon-photon absorption process. 

At this point you may ask if pairs are indeed produced in some class of 
sources... The answer is yes, even if the reprocessing features outlined above 
do not have (yet) a clear application to specific sources. At the beginning 
of the '90s, it was possible to detect TeV photons from the ground, through 
optical Cherenkov telescopes. Soon, low power blazars (i.e. BL Lac objects) 
were found to strongly emit at these frequencies. Some of them are distant 
enough that the interaction of TeV photons with the cosmic IR background 
photons is not negligible. We do see signs of the corresponding absorption. 
The importance of the latter depends, of course, of the density of the target 
photons, i.e. the level of the cosmic IR background. Since the latter is 
not precisely known, one can use the absorption features of the observed 
spectrum to estimate the background itself. 
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Chapter 8 

Active Galactic Nuclei 



8.1 Introduction 

Up to the seventeenth century, the unaided eye was the only receiver that 
humanity could use to observe the Universe. Evolution was able to apt 
this "instrument" to be sensitive to the light of the star we happen to be 
orbiting, the Sun. The invention of the telescope amplified the sensitivity 
of the eye and its angular resolution, letting humanity discover, less than 
a century ago, that other galaxies exist, far beyond out Milky way, and 
that these galaxies are moving apart: the Universe expands. However, all 
we could discover using the eye and its extension, the optical telescope, 
concerned a tiny, very tiny, part of the entire electromagnetic spectrum. As 
soon as technology enabled humanity to open new windows, we discovered 
other phenomena, other objects, and could push our knowledge farther out 
in space and in time. 

The opening of the radio window made the sixties the golden decade for 
astronomy, with the discovery of the microwave background and of pulsars. 
The third great discovery made in that decade was that of quasars. The 
term "quasar" originally stood for "quasi-stellar radio source", and refers 
to the fact that when an optical telescope is pointed towards the direction of 
some radio source, which can be as extended as minutes of arc in radio maps, 
the resulting optical plate shows a source which looks like a star, i.e. a not 
extended, a "pointlike" object. This apparently innocuous point is instead 
a gigantic energy plant, able to produce much more power that an entire 
galaxy like our own, in a volume which is ridiculously small, if compared 
with the dimension of the Galaxy, and comparable with our Solar System. 
The process that powers the stars, thermonuclear reactions, is not enough to 
power quasars. For this reason we believe that at the core of these sources 
there is a massive black hole, with a mass between a million and a billion 
the mass of the Sun. Matter around the hole is attracted by the black hole 
gravity, it is compressed, heated, and then radiates. This was realized quite 
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soon (Salpeter 1964; Zeldovich 1964; Lynden-Bell 1969; Shakura & Sunjaev 



Another major advance came with the opening of the X-ray window, 
first (in the sixties) with rocket experiments pioneered by Riccardo Giacconi, 
Bruno Rossi and others, and then with the first X-ray satehites, in the early 
seventies. The Uhuru, Ariel 5, HEAO-1 and then Einstein satellites made 
clear that all kinds of quasars were strong X-ray emitters: at the same time, 
people started to believe that quasars were the major contributors to the 
cosmic diffuse X-ray background, already discovered in 1962 (Giacconi et 



A third qualitative "jump" was the improvements of the interferometric 
technique in the radio band, in the early seventies. Radio telescopes in 
different continents, looking at the same source, can resolve details as close 
as a few tenths of a millisecond of arc. This enabled us to discover that 
some radio emitting quasars have spots of radio emission which are observed 
to move. Given the huge distances, this motion corresponds to velocities 
that exceed the speed of light. Far from challenging special relativity, this 
"superluminal" motion, as it is now called, was even predicted before it was 
observed, by Martin Rees in 196(jl, and corresponds to real fast motion (but 
slower than the velocity of light!) at an angle close to our line of sight. 

8.2 The discovery 

8.2.1 Marteen Schmidt discovers the redshift of 3C 273 

Prom Maarteen Schmidt (Schmidt 1990 

"The puzzle was suddenly resolved in the afternoon of February 5, 1963, 
while I was writing a brief article about the optical spectrum of 3C 273. Cyril 
Hazard had written up the occultation results for publication in Nature and 
suggested that the optical observations be published in an adjacent article. 
While writing the manuscript, I took another look at the spectra. I noticed 
that four of the six lines in the photographic spectra showed a pattern of 
decreasing strength and decreasing spacing from red to blue. For some 
reason, I decided to construct an energy-level diagram based on these lines. 
I must have made an error in the process which seemed to contradict the 
regular spacing pattern. Slightly irritated by that, I decided to check the 
regular spacing of the lines by taking the ratio of their wavelengths to that 

^ Lynden-Bell D., 1969, Nature, 223, 690; 
Salpeter E.E., 1964, ApJ, 140, 796; 
Shakura N.I. & Sunjaev R.A., 1973, A&A, 24, 337; 
Zeldovich Ya. B., 1964, Dock. Akad. Nauk SSSR, 157, 67 

^ Giacconi R., Gursky H., Paolini F., Rossi B., 1962 Phys. Rev. Lett., 9, 439. Giacconi 
was awarded the Nobel Prize in Physics in 2002. 

^ Rees M.J., 1966, Nature, 211, 468; Martin Rees was born in 1942. He was 24 in 1966. 
Schmidt M., 1990, in Modern Cosmology in Retrospect; (Cambridge University Press) 
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of the nearest line of the Balmer series. The first ratio, that of the 5630 hne 
to H-/3, was 1.16. The second ratio was also 1.16. When the third ratio was 
1.16 again, it was clear that I was looking at a Balmer spectrum redshifted 
by 0.16. 



I was stunned by this development: stars of magnitude 13 are not sup- 
posed to show large redshift! When I saw Jesse Greenstein minutes later 
in the hallway and told him what had happened, he produced a list of 
wavelengths of emission lines from a just completed manuscript about the 
spectrum of 3C 48. Being prepared to look for large redshift, it took us only 
minutes to derive a redshift of 0.37. 



The interpretation of such large redshifts was an extraordinary challenge. 
Greenstein and I soon found that an explanation in terms of gravitational 
redshift was essentially impossible on the basis of spectroscopic arguments. 
We recognized that the alternative explanation in terms of cosmological 
redshifts, large distances, and enormous luminosities and energies was very 
speculative but could find no strong arguments against it. The results for 
3C 273 and 3C 48 were published six weeks later in four consecutive articles 
in Nature (Hazard, et al. 1963; Schmidt 1963; Oke 1963; Greenstein and 
Matthews 1963 j§' . 

8.3 Basic components 

The first quasars that were discovered were radio-loud, but we now know 
that these radio-loud AGNs are only roughly the 10% of all AGNs. The 
majority are radio-quiet. 

The basic structure of AGNs includes: 

• A black hole, with 10^ < M < W^^Mq. It is probably spinning 
at some level, even if we do not have secure measurements of the spin 
value. 

• An accretion disk. Matter with even a small amount of angular mo- 
mentum, attracted by the black hole gravity, spirals in and forms a 
disk. This is a major source of power. 

• An X-ray corona, sandwiching the accretion disk. It is supposed to 
be a hot layer, or an ensemble of clumpy regions particularly active in 
the inner parts of the disk. 

Greenstein J.L. & Matthews T.A.,1963, Nature, 197, 1041, 
Hazard C, Mackey M.B. & Shimmins A.J., 1963, Nature, 197, 1037; 
Oke J.B., 1963, Nature, 197, 1040; 
Schmidt M., 1963, Nature, 197, 1040. 
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• An obscuring torus located at several parsec from the black hole, in- 
tercepting some fraction of the radiation produced by the disk and 
re-emitting it in the infrared. 

• A region of many small clouds at a distance of ~ 10^''-10^^ cm from 
the hole (i.e. less than a parsec) moving rapidly (~ 3000 km s"^). 
They intercept ~10% of the ionizing radiation of the disk, and re- 
emit it in the form of lines. Doppler shifts broaden the observed lines. 
This is why this region is called Broad Line Region (BLR). 

• At larger distance (~100 pc) there is another region where less dense 
clouds are moving, less rapidly. This is called Narrow Line Region 
(NLR). 

• About 10% of AGNs, besides accreting matter, are able to expel it in 
two oppositely directed jets. Their direction likely traces the rotational 
axis of a spinning black hole. The material inside these jets is moving 
at relativistic speeds. Therefore the jet emission is highly beamed, 
and their appearance depends on the viewing angle. AGNs whose jets 
are pointing at us are called blazars. AGNs whose jets are pointing 
elsewhere are called radio- galaxies. 

8.4 The supermassive black hole 
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Figure 8.1: The mass of the supermassive black hole correlates with the 
luminosity of the bulge of the host galaxy and with the velocity dis- 
persion of the stars of the central regions of the host galaxy. From 
http:/ / chandra.as.utexas.edu/~kormendy/bhsearch.html. 
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The black hole in the center of AGNs has a mass that "knows" of the 
mass of the host galaxy. In fact there is a correlation between the black 
hole mass and the luminosity of the bulge, due to the stars (see Fig. 18.1^ 
left panel). The black hole mass correlates even more with the velocity 
dispersion of these stars (see Fig. 18.11 right panel). The correlation with 
the bulge luminosity is linear: Mbh ~ 2 x 10~^Mbuige) while the correlation 
with the velocity dispersion a is between Mbh and Mbh oc o"^. 



8.5 Accretion disk: Luminosities and spectra 



46 



45.5 




44 I 1 I I I I I I I L I I I I I I 

12 14 16 18 

Log u [Hz] 

Figure 8.2: Spectrum from a standard accretion disk, the infrared torus and 
the X-ray emitting corona. 

The luminosity of the accretion disk depends on the amount of the mass 
accretion rate M. It is usually measured in solar masses per year. The 
process has a very large efficiency rj, defined as the fraction of rest mass 
energy that is converted into radiation, before being swallowed by the black 
hole. We have 

idisk = vMc"" (8.1) 

The efficiency r] depends on the details of the accretion mode, but it is of 
the order of ~ 0.1. To produce a disk luminosity of L^isk = W^^L^q erg 
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s ^ we then require 



M = ~ 1.75 L46 Mqyi-^ (8.2) 



There is a Hmit on the disk luminosity. This is called Eddington luminosity. 
It corresponds to radiation pressure balancing gravity. Radiation pressure 
acts on electrons, gravity acts both on electrons and their companions pro- 
tons, that are heavier. 

^Edd GMmp iircGMmp , o ^ ^38 ^ -1 
-7—52- = — 52 ^Edd = = 1-3 X 10^*— erg s 

(8.3) 

The Eddington luminosity depends only on the black hole mass. If the 
accretion rate for some reason produces a luminosity Ldisk > -^^Edd) then 
the accreting mass feels a radiation pressure that contrasts gravity, and 
the accretion stops, making L^isk to decrease below i^Edd- We have a self- 
regulating process. 
Ask yourself: 

• What implicit assumptions have we made to derive LEdd? 

• Did we neglect some other forces? 

• Why is there (Tt? 

• What happens if we have electron-positron pairs? 

• Can you envisage cases in which the Eddington limit can be avoided? 



8.5.1 Spectrum 

Accreting matter, by spiraling inside, has to dissipate angular momentum 
and potential energy. The potential energy of a proton at distance R (much 
larger than the Schwarzschild radius Rs = 2GM/c^) is 



= — ^ (8.4) 



namely, it increases (or decreases, if you take it negative) as R decreases. 
More energy has to be dissipated close to the black hole than far from it. 
Besides, the surface 2 x 2'KrdR (there are two faces of the disk...) of the ring 
at a distance R and width dR becomes smaller and smaller as you move close 
to the black hole. More energy has to be dissipated from a smaller surface. 
Assuming that the produced disk radiation is a blackbody, this inevitably 
implies that the temperature must increase as R decreases. If all rings of 
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the disk emit black-body radiation, the temperature profile is: 



T = 



1/4 



oc 



3i?S^disk 

IQirrjaMBR^ 
i?-^/^, for R > 3Rs 



3Rs 
R 



1/2- 



1/4 



(8.5) 



where o"mb is the Maxwcll-Boltzmann constant. The term in the second 
square brackets allows for the fact that not all the orbits are stable around 
a black hole, but only those beyond the marginally stable orbit: for a 
Schwarzschild black hole this orbit is at 3i?s- Matter at i? < 3i?s falls 
directly onto the hole in a short time. If we measure R in units of the 
Schwarzschild radius Rs, x = R/Rs we have: 



T = 



647r(7MB2;^G2M2 
MV4 



1/4 



1/2- 



1/4 



OC 



MV2 



X 



-3/4 



for X 3> 3 



(8.6) 



We can now have an idea of the emitted spectrum. It will be the sum 
of black-bodies peaking at different temperatures: each ring of width dR 
emits a luminosity dL: 



dL = 2x 2'KaMBRT^dR 
We can then set: 



dL 
dR 



= 2 X 27r(7MBr-^/^T'' 



dL dLdR rp-4/3rp4dR _ ^1/3 

dT = dRdT ^'"''^^ ^dT^"^ 



(8.7) 



(8.8) 



where we have used R oc T^^/^ and then (dR/dT) oc T"''/^. This is valid as 
long as S> 3Rs. Since there is a one to one correspondence between the 
temperature and the peak frequency of the black-body, we have that 



Ldisk{i^) oc u^^^; hv < feT„ 



(8.9) 



up to a frequency corresponding to the maximum temperature. Beyond 
that the emission is not a superposition of black-bodies anymore, and only 
the region at the maximum temperature contributes. Therefore, at large 
frequencies, we will sec an exponential drop. At the opposite side (small 
frequencies) we see only the black-body produced by the outer radius of the 
disk i?out) and correspondingly, we see the Raleigh-Jeans emission of that 
radius. Therefore we can approximate the disk spectrum as: 



Ldisk oc 



hv < kT{Rout) 

hv > kT{Ront) 



(8.10) 
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Figure 8.3: The spectrum from a standard accretion disk, emitting black- 
body radiation from all its surface, with a temperature profile given by Eq. 
18.51 The disk extends from Ri^ = 3Rs to different values of i?out, as labelled. 
Note that the u^^^ slope is present only for very large values of i?out- 



The peak is determined by Tmax- We can estimate it by assuming that most 
of the luminosity is produced within -say- lORs- 

2 2 4 f Ldisk ^ 

\4U07r/Cg(TA/B 



Ml, 



3.5 X loM ^ ] K (8.111 



where Mg is the black hole mass in units of IO^Mq. This temperature 
corresponds to a frequency 

^peak ~ ~ 7.4 X 10^^ Hz (8.12) 

i.e. in the near ultraviolet. For disks emitting at the same Eddington ra- 
tio (i.e. -^^diskZ-^^Edd =const), we have t'peak oc M~-^/^. Fig. 18.31 shows the 
spectrum calculated numerically according to the assumed temperature dis- 
tribution of Eq. 18.51 for different values of the outer radius of the disk i?out • 
Note that the v^l"^ behaviour can be seen only if i?out is very large. 
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8.6 Emission lines 

A distinctive characteristic of AGNs is the presence, in their optical spectra, 
of emission Hnes. They comes in two flavors: the broad and the narrow. 

8.6.1 Broad emission lines 

Fig. 18.41 shows an atlas (i.e. the average spectrum of many sources) of the 
optical-UV spectrum of typical AGNs. Typical lines are labelled. The most 
prominent are the Hydrogen Lya (transitions from n=2 to n=l), the lines 
from partially ionized Carbon, Magnesium, Oxygen and then the Hydrogen 
Ha (transitions from n=3 to n=2). If you look carefully, you should notice 
that most of these lines are broad, while the [OIII] line is narrow. The 
square brackets indicate that this is a "forbidden" line. So why do we see 
it? The reason to call this line "forbidden" is that when the density is larger 
than some critical value, the de-excitation of the excited state is made by 
collisions, occurring faster than the radiative ones. In our labs the densities 
are always larger than the critical ones, so in the lab we do not see the line. 
This gives some hints of the density where the narrow lines are made. 
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Figure 8.4: Atlas optical/UV spectrum for optically-selected quasars in the 
Sloan Digital Sky Survey, from the work of Vanden Berk et al. (2001, AJ, 
122, 549). Notice the presence of broad lines and also of the forbidden [OIII] 
narrow line. 
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The most used way to quantify the width of a line is to measure it at a 
flux level which is half of the peak flux. In this case we speak of Full Width 
at Half Maximum (FWHM). Full width means that it is the entire width, 
not only one half as in the case of a cr of a Gaussian. The units can be A or 
km s~^. In fact the width can be related immediately with a Doppler shift, 
and thus to a velocity. 

The typical values for the broad emission lines are: 

• FWHM: from 1000 to 10,000 km s"^ They arc too large to be due 
to thermal motions. In fact, with a temperature of T ~ 10^ K, we 
have that thermal broadening implies: 

/LT-X 1/2 

V ^ [ ~ 10^ cm s-^ = 10 km s'^ (8.13) 

Therefore we need bulk motions (in different directions with respect 
to our line of sight) of the material emitting the lines. 

• Density: between 10^ and lO^"*^ cm~^. This comes from the presence 
or absence of specific lines. 

• Temperature: around 10^ K. Again, this comes from the presence 
or absence of specific lines. 

• Covering factor: around 0.1. A simple estimate comes from the 
luminosity of all the broad lines with respect to the continuum. 

• Distance: depends on the luminosity of the ionizing continuum. Re- 
cent studies use the technique of "reverberation mapping": taking 
spectra of a source regularly, one can see the time delay between a 
variation of the continuum and the variations of the lines. Prom these 
studies one gets a size: 

i?BLR ~ W^'Ll^l.^ ~ 10^'4(sM5 cm (8.14) 

Where we have assumed that most of the disk radiation can ionize the 
atoms in the BLR. Please note that this is only approximately true. 

One useful quantity when studying lines is the "ionization parameter 
^" . ^ can be defined as the ratio of the densities of the ionizing photons 
and the electrons: 



1 f°° L(u) , i Ionization rate ,^ , ^, 

^ = — J; \ du = ^ oc : (8.15) 

Ue J Hi ^irR'^chv rip. Recombination rate 



So, what is then the Broad Line Region? We do not know exactly yet. 
It can be an ensemble of clouds in Keplerian motions around the black hole, 
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or material that is infalling, or outflowing, with a range of velocities. Broad 
lines could also be produced within the accretion disk itself, but in this case 
we should see "double horned" lines, because for a typical inclination of the 
disk, the opposite sides of the disk produce redshited and blueshifted lines. 
They are seen sometimes, but not often. It has been proposed also that 
broad line clouds are stars with bloated envelopes, coming close to the black 
hole. But a large number of them is required, possibly too large. If broad 
lines are due to individual clouds, we must have many of them, because each 
one contributes for only 10 km s~^ to the broadening, and if they were a few, 
we should see a not completely smooth profile of the line. Estimates range 
from 10^ to billions. But the total required mass is tens of solar masses, a 
very small quantity. 

8.6.2 Narrow emission lines 

The main properties of the narrow emission lines are: 

• FWHM: around 300-500 km s'^. These hues are "narrow", but in 
any case much larger than thermal broadening. 

• Density: around 10^ cm^'^. This comes from the presence of forbid- 
den lines, requiring small densities. 

• Temperature: around 10^ K. This comes from the presence/absence 
of specific lines. 

• Covering factor: around 0.01. This comes from the entire NLR 
luminosity divided by the continuum luminosity. 

• Distance: around 100 parsec. They can be resolved in nearby objects. 
Sometimes they appear contained in a double cone. 

Since the NLR is located at greater distances than the BLR, they are un- 
affected by the possible presence, in a region close to the accretion disk, of 
absorbing material. This is a crucial fact, as we will soon see. 

8.7 The X— ray corona 

From Eq. I8.5l it is clear that the standard accretion disk is relatively "cold", 
in the sense that even the inner regions cannot reach ~keV temperatures (1 
keV ~ 10^ K). On the other hand, AGNs are powerful X-ray emitters, and 
their X-ray spectrum extends out to hundreds of keV. We must assume that 
there is an additional component, besides the accretion disk, responsible for 
this emission. Observed in the energy range 0.1-10 keV, with low spec- 
tral resolution, the X-ray spectrum can be approximated by a power law: 
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Fx{v) ~ with ax ~ 0.7-0.9. At first sight, it seems a non-thermal 

spectrum, because it is characterized by a power law. But then more accu- 
rate observations were made, that revealed three important features: 

• When high energy observations in the 10-200 keV were available, a 
cut-off was observed, of an exponential type. The spectrum was thus 
a power law ending with an exponential cut: 



where Vc is somewhat different for different objects, and it is in the 
range 40-300 keV. 

• When spectral resolution was improved in the soft X-rays (soft means 
0.1-10 keV), an emission line appeared at ~ 6.4 keV (see Fig. 18. 5p . 
This is the Ka line from cold Iron. It is a fluorescent line. It means 
that a photon of enough energy interacts with one electron of the inner 
shell of the Iron atom, and ionizes it. There is then a "vacant place", 
triggering an = 2— )-n = l transition. The energy of this transition 
depends on the ionization state of the Iron. If it is partially ionized, 
the energy of the jump is more. So the energy of the line tells us about 
the ionization state of the Iron. For almost completely ionized Iron, we 
should observe a line at 6.7 keV (and sometimes we do). A line at 6.4 
keV means that it is not largely ionized, and therefore it is relatively 
cold (i.e. less than 10^ K). 

• At the same time, the improved spectral resolution revealed that the 
slope of the spectrum was more complex than a simple power law. 
There is the presence of a hump, peaking at roughly 30 keV, super- 
posed to the power law. So the previously determined slope of the 
power law was not correct, because it included the hump. Depurating 
from this extra component, the spectral index of the power law be- 
comes ax ~ 0.9-1. The right panel of Fig. 18.51 illustrates this point. 
The nature of this 30 keV "hump" was soon interpreted as the emission 
coming from an irradiated disk, and was baptized "Compton hump" 
or "Compton reflection" . 

We will now discuss in somewhat more details these three components: i) 
the power law; ii) the Compton reflection, and iii) the iron line. 

8.7.1 The power law: Thermal Comptonization 

We have seen in the previous chapters that a power law does not necessarily 
imply a non-thermal emission. Bremsstrahlung and thermal Comptoniza- 
tion are made by thermal electrons, but their shape is a power law (ending 
in an exponential cut). Following this idea we will now postulate that there 
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Figure 8.5: Left: The "Compton reflection hump", together with the pro- 
duced emission hnes. Right: X-ray spectrum and its components: a power 
law continuum, a soft excess, and a "Compton reflection" , together with the 
Iron Ka hne. From Fabian & Miniutti 2005 (astro-ph/0507409). 

is a region above and below the disk, where the electrons are hot, much 
hotter than the disk. It is something similar to the corona of our Sun, and 
therefore this region is called X-ray corona. It cannot be located far from 
the black hole for several reasons: 

1. Its luminosity is comparable (even if smaller) than the luminosity of 
the accretion disk. If its power depends on the release of the grav- 
itational energy, it must be close to the black hole (i.e. remember 
Eg oc l/R). 

2. The X-ray flux varies rapidly. 

3. It must illuminate the disk and produce the ionization flux for the 
production of the fluorescent Iron line. As we shall see, the Iron line 
is sometimes very broad, as a result of Doppler broadening due to the 
high Keplerian velocities in the inner regions of the disk. So also the 
"illuminator" must be close. 

Can we derive some more information from the spectral shape? We can 
exclude bremsstrahlung as the main radiation process because it is too inef- 
ficient, meaning that we would require too many electrons to produce what 
we see. Variability tells us that the emitting volume is small, so the elec- 
tron density would be large, and the Thomson optical depth would be much 
larger than unity. In these conditions, Compton scatterings would be more 
important. Therefore let us assume that the spectrum is due to thermal 
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Comptonization. The temperature can be inferred from the exponential 
cut: 

For each scattering, the energy gain A = xi/x of the photon is 

A = lee^ + 49 + 1 ~ 1.9 (for 9 = 0.2) (8.18) 

If tt < 1, to be checked a posteriori, we have a relation between the spectral 
shape, the amplification factor A and the optical depth: 

ax = log^T = -ax log A — > tt = = 0.5 (8.19) 

where the last equality assumes ax = 1 and A ^ 2. The Comptonization 
parameter y must be of the order of unity, because ^ ~ y/rx. What derived 
here are typical values. There are sources where we need tt > 2: for these 
sources the simple equation we have used to relate ax, tt and A is not 
valid, and another relation must be used. 

Why the corona is so hot, even if the disk is cold? What are the mecha- 
nisms able to extract energy from the accreting matter and to deposit it in 
the corona? We do not know exactly, even if we suspect that the magnetic 
field plays a crucial role. 



8.7.2 Compton Reflection 

If the corona emits isotropically, half of its flux will be intercepted by the 
disk. We have seen that the disk in AGNs should be relatively cold, in the 
sense that most of its material is not ionized. X-ray photons will interact 
with the material in the disk in the following ways: 

1. at low energies {hv <10 keV) they will be photoelectrically absorbed 
by the metals present in the disk; 

2. at intermediate energies (10 < hiy < AO keV) they will be Thomson 
scattered and part of them will be scattered in the upward direction; 

3. at large energies {hu > 40 keV) the scattering will not be in the Thom- 
son regime, but Klein-Nishina effects will start to be important: the 
scattering is preferentially forward directed (and therefore the scat- 
tered photons will penetrate more deeply into the disk) and the pho- 
ton energy will be reduced. The decrease of the energy means that 
photoelectric absorption, initially negligible, is again important, and 
the photon will be absorbed. 

The result is that some of the incoming radiation will be scattered back, 
with a modified shape. This Compton reflection component will be peaked 
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at ~30 keV, where photoelectric absorption and Klein-Nishina effects are 
not important. See Fig. l8.5l to see the shape of this component, and how it 
modifies the original power law spectrum. 

If the disk is hot, and partially ionized, the amount of free electrons will 
increase the importance of Thomson scattering, and the "left shoulder" of 
the Compton reflection increases. Increasing the ionization state of the disk 
implies that the Compton reflection is not peaked any longer, but it retains, 
at low energies, the slope of the original power law. 

8.7.3 The Iron Line 

Among the metals in the disk irradiated by the hot corona. Iron is the one 
producing the most prominent emission line. This is due to a combination 
of two factors: it is one of the most abundant elements (but not the most 
abundant), and it is the one "suffering" less from the Auger effect. This 
is, in very simple terms, a non-radiative transition: the energy of the n = 
2 — )• n = 1 transition is used to expel an electron from the atom instead 
of producing a photon. See the left panel of Fig. 18.51 to see the relative 
strength of the emission lines. The energy, luminosity, width, and profile of 
the iron line are a powerful diagnostic of the condition of the inner parts of 
the accretion disk. 

The energy of the Ka Fe line (from 6.4 to 6.7) tells about the ionization 
state of the iron, and thus the temperature of the disk. 

The luminosity of the line tells about the amount of iron, and thus about 
the abundance of metals of the disk. The ratio about the line to contin- 
uum luminosity tells about the geometry of the emitting region and possible 
anisotropy of the continuum. 

The width of the line tells about the velocities of the irradiated material 
forming the line. 

The profile (symmetric, double horned, skewed) tells about Doppler boost- 
ing and gravitational redshift. 

The studies of iron lines in AGNs received a lot of attention in the recent 
past. The most intriguing observation concerns the finding of relativistically 
broadened lines. Fig. 18.61 shows in a schematic way what happens. Suppose 
to see an accretion disk from an intermediate angle (i.e. neither face on or 
edge on). If the emitting material is moving with a non relativistic veloc- 
ity (Newtonian case, upper panel) we see a symmetric double peaked line, 
because the parts of the disk that are approaching us emit a line that we 
observe blueshifted, while the parts of the disk moving away from us emit 
a line that we observe redshifted. The symmetry refers to the fact that the 
blue and red lines have equal fiuxes. But if the material is orbiting very close 
to the black hole, the Doppler enhancement of the fiux becomes important, 
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Figure 8.6: The profile of an intrinsically narrow emission line is modified 
by the interplay of Doppler/gravitational energy shifts, relativistic beam- 
ing, and gravitational light bending occurring in the accretion disc. Upper 
panel: the symmetric double-peaked profile from two annuli of a nonrela- 
tivistic Newtonian disc. Second panel: the transverse Doppler shifts make 
the profiles redder; the relativistic beaming enhances the blue peak with 
respect to the red. Third panel: the gravitational redshift shifts the overall 
profile to the red side and reduces the blue peak strength. Bottom panel: 
The combination of all the effects gives rise to a broad, skewed line profile. 
From Fabian & Miniutti (2005) (astro-ph/0507409). 



and the blue peak will have more fiux than the red one (second panel from 
the top). If the emission takes place very close to the black hole (within a 
few Schwarzschild radii) also the gravitational redshift is important: then 
all frequencies are redshifted by an amount that is larger as the emitting 
regions are closer to the black hole (third panel). Furthermore, we have to 
account for gravitational light bending, changing the emission pattern and 
thus the received flux. 

All these effects have opened the way to the possibility to test general 
relativity in the strong field case. 

Besides this, there is another, crucial effect, related to the spin of the 
black hole. This is conceptually simple: for a Kerr black hole fastly rotating 
the innermost stable orbit is not 3Rs (i.e. 6 gravitation radii = 6Rg), but 
it moves closer to the hole, becoming equal to one gravitational radius for 
a maximally rotating hole. The closer to the hole the emitting material is. 
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Figure 8.7: The Ka Fe line from a Schwarzschild (no spin) and a max- 
imally rotating Kerr black hole. From Fabian & Miniutti (2005) (astro- 
ph/0507409). 

the larger the velocities, and the stronger the gravitational effects. So the 
broadening of the line becomes larger, as it is illustrated in Fig. 18.71 We 
can tell the spin of the black hole by observing the profile of the iron line. 
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8.8 The torus and the Seyfert unification scheme 

Seyfert galaxies are spiral galaxies with a bright core, that host an Active 
Galactic Nucleus. They come in two flavors: the Type 1 and the Type 2. 
Type 1 Seyferts have both broad and narrow emission lines, while Type 2 
have only narrow lines. See Fig. 18.81 




4000 5000 6000 

Figure 8.8: The optical spectrum of a Seyfert Type 1 and Type 2. Note that 
Type 1 Seyferts have both broad and narrow lines, while Type 2 Seyferts 
have only narrow lines. 

Antonucci & Miller (1985 j£] made spectropolarimetry observations of 
NGC 1068, a nearby and bright Seyfert 2 galaxy. In total light the broad 
lines were not seen, but in polarized light they emerged. Therefore the broad 
line clouds are there, but for some reason their emission is overwhelmed by 
the unpolarized continuum, and in total light their emission is not seen. 
Antonucci & Miller then argued that this is due to obscuration. Fig. 18.91 
illustrates the point. If there is a big and dusty torus surrounding the accre- 
tion disk, and the observer is looking from the side, the emission from the 
disk and the BLR is absorbed by the dust in the torus. Instead the narrow 
lines, coming from an extended region ~100 pc in size, are not intercepted 
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by the torus, and can reach observers located at any viewing angle. 

Furthermore, if the "funnel" of the torus is filled with free electrons, they 
can scatter the emission of the disk and the broad lines also into the direction 
of side observers. The amount of this scattered radiation can be tiny, but 
it is polarized (remember: for scattering angles close to 90° the scattering 
photons are maximally polarized). Observing the polarized light we select 
these photons and discard the much more intense unpolarized light: as a 
result the broad lines now appear. The free electrons act as a periscope. 




Figure 8.9: Schematic illustration of the disk+torus system, with the broad 
line clouds that are close to the disk, and the narrow line clouds external to 
the torus. Narrow lines are then observable at all viewing angles, while side 
observers cannot see the disk and the BLR. However, if a populations of 
free electrons fill the region (black little dots), they can scatter the emission 
line photons into the line of sight. Since the scattering angle is close to 90°, 
the scattered photons are polarized. Observing in polarized light the broad 
emission lines emerge. 

This simple idea implies that all Seyferts are intrinsically equal, but 
their appearance depends on the viewing angle: observers looking face on 
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can see all the components (and thus also the broad lines): they see a Type 
1 object. Observers looking from the side see only the narrow lines and 
the emission of the torus. Being at some distance from the accretion disk, 
the torus in fact intercepts part of the disk emission and re-emits it in the 
infrared. We can even predict approximately at which frequencies the torus 
emission will start to be important. This is because the dust cannot exist at 
temperatures hotter than ~2000 K. Taking this temperature as reference, 
we have ~ 3kT/h ~ 10^'* Hz, which corresponds to a wavelength of 3 
/J,. If we are correct, then we predict that all the radiation that the torus 
intercepts from the disk comes out at frequencies smaller than 10^^ Hz. 

8.8.1 The X— ray background 

Another very important consequence of the presence of the torus concerns 
the X-ray emission, because also X-rays are absorbed by the torus, espe- 
cially in the soft energy range (soft X-rays are much more absorbed than 
hard X-rays). In addition, the walls of the torus will also "reflect" X-rays, 
producing a "Compton reflection hump" , like the disk. The main difference 
is that the fluorescent iron line will in this case be extremely narrow (since 
the typical velocities of the torus material, located far from the black hole, 
are very small). 

100 II I I I I 1 1 — I — I — I I I I I 1 1 — I — I — I I I I I 1 1 — I — 



Figure 8.10: The X-ray background. The lines are synthesis models as- 
suming some prescriptions about the intrinsic slope of the X-ray spectrum 
and the ratio of "thick to thin" AGNs. From Gilh et al. (2007) (astro- 
ph/0704.1376). 
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The importance of the absorption depends on the metals present along 
the line of sight. For solar abundances, there is a well defined relation be- 
tween the amount of metals and the amount of hydrogen. So it is customary 
to express the amount of absorption in terms of the "column of hydrogen" 
along the line of sight. It is A'^h = ^H^i where R is the size of the absorbing 
region, and nn is the hydrogen number density. The units of N}^ are thus 
cm~^. When A^h > 10^'' cm~^ we have that tt = ctt-^Vh becomes close or 
larger than unity, and we speak of Compton thick sources. 

With these properties, a mixture of Type 1 and Type 2 objects can 
account for the X-ray background^ namely, the integrated X-ray emission 
from all sources. It was a mystery for decades, because its spectrum in the 
0.5-10 keV energy range is very hard {ax ~ 0.4), unlike most of the known 
extragalactic X-ray sources (that have ax ~ 1). Fig. 18.101 shows a modern 
collection of data defining the X-ray background from 0.5 to 300 keV. 

To explain it, we need a population of very hard X-ray sources and the 
Type 2 ones have indeed the right properties: since soft X-rays are more 
absorbed than hard ones, their spectrum is very hard. Since the X-ray 
continuum is depressed, but the reflection hump (from the torus) is not, they 
can contribute significantly at ~30 keV, just where the X-ray background 
peaks. 
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8.9 The jet 

About 10% of AGNs, besides accreting matter, are able to expel part of it at 
relativistic speeds in two oppositely directed jets. These radio-loud AGNs 
were the first to be discovered, due to their relevant radio emission. Radio 
jets can be spectacular, since they can reach a few Mpc in size, tens of times 
the radius of their host galaxies. AGNs with relativistic jets were believed 
to always have giant ellipticals as host galaxies, but now this paradigm is 
challenged by the observations of Narrow Line Seyfert 1 Galaxies (beware 
that the term "narrow" here indicates lines that are in any case broad, but 
"less" broad than usual, namely with FWHM< 2000 km s^^). 

The radio emission produced by jets is only a small fraction of the en- 
tire electro-magnetic power they emit. We now know that most of it is 
produced in the range mm-optical and in the 7-ray band. In turn, the 
electro-magnetic output is only a small fraction of the total power carried 
by the jet. Most of it is carried in the form of bulk kinetic energy of the mat- 
ter flowing relativistically, and by the moving magnetic field (i.e. Poynting 
flux); in powerful sources it reaches the large radio structure: the hot spots 
and the radio lobes. In less powerful sources, however, these structures are 
absent. We call the powerful sources FR //radio-galaxies, and the weaker 
sources FR I radio-galaxies. The name FR comes from Fanaroff & Riley, 
that classified radio-galaxies in this way. 

8.9.1 Flat and steep radio AGNs 

Since the emitting material is moving with a bulk Lorentz factor F, its 
radiation is beamed. What we see is therefore amplified if the jet points at 
us, and depressed if the jet points elsewhere. The large structures (hot spots 
and radio-haloes), instead, are not moving, and their emission is isotropic. 
The ratio between the jet emission and the radio-lobe emission is thus a 
strong function of the viewing angle. 

Fig. 18.111 illustrates this point: the jet is producing a very fiat spectrum 
[i.e. F{ij) oc u^] up to 10^^ Hz. The observed intensity is I^et{i^) 6^ 
and changes a lot changing the viewing angle. The emission coming from the 
extended, unbeamed, region is characterized by a steeper slope [for instance: 
hohei'^) Of '^^^]- Therefore: 

• At very low frequencies we expect that the emission is always domi- 
nated by the lobe, for all viewing angles. 

• Lobe and jet emission are equal at some frequency that is larger as 
the viewing angle increases. 

• For misaligned jets the depressing factor due to beaming is so large 
that we do not see the jet, we see only the lobe emission. 
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Figure 8.11: The emission of the jet is beamed, and depends on the viewing 
angle. Since it is a superposition of emission from several blobs of different 
sizes, the spectrum is flat up to the self-absorption frequency of the most 
compact component (hundreds of GHz). The emission from the radio-lobe 
is instead isotropic, and is optically thin. Since it is not beamed, it does not 
depend on the viewing angle. The overall appearance of the radio spectrum 
depends therefore on the viewing angle. Steep radio spectra correspond to 
large 0^ (and small 6), flat radio spectra to small 9^ (and large 5). 

• Aligned sources show a flat radio spectrum. Misaligned sources show 
a steep radio spectrum. The slope of the radio spectrum is thus an 
indication of the viewing angle. 

• If we perform a radio survey with some limiting flux, the number of 
aligned and misaligned sources changes according to the observing fre- 
quency: if we use 178 MHz, we preferentially look at the lobe emission, 
and we can pick up aligned and misaligned sources; if we observe at 
30 GHz we pick up preferentially the flat, aligned sources. 

To summarize, and simplifying a bit: 

flat radio spectrum — > aligned, beamed 

steep radio spectrum — > misaligned, de beamed. 
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8.9.2 Blazars 

Sources whose jet is pointing at us are called blazars. The sub-classification 
of blazars is rather complex, reflecting either the way they were discovered 
(radio or X~rays) or the location of their synchrotron peak (in vFy), or the 
presence or absence of broad emission lines. Particularly relevant in this 
respect is the concept of equivalent width of a line (EW). It is defined as 



EW = 




(8.20) 



where F\ is the total flux (line + continuum) and Fq is the flux of the 
continuum. This describes the EW of both absorption and emission lines, 
but in the case of emission lines the result is negative, so sometimes the 
absolute value is given. The units of the EW are A. 

The classical subdivision of blazars is 
BL Lac objects (BL Lacs) — They have weak or absent lines. Their 
equivalent width is EW< 5 A. 

Flat Spectrum Radio Quasars (FSRQs) — They have strong emission 
lines, whose EW> 5 A. 

Since we are considering very variable objects, this subdivision depends on 
the state of the source when it is observed: this is because the EW is a ratio 
between the line and (beamed) optical flux, and the latter varies. Now we 
believe that the absence or weakness of the lines in BL Lacs is not due to 
a particularly amplified continuum, but it is an intrinsic property. In other 
words, the luminosity of the lines in BL Lacs is less than in FSRQs. 



The overall spectrum of blazars 

Fig. 18.121 shows the source of 3C 454.3, a powerful FSRQ, that was the 
brightest blazar detected in the 7~ray band for a few years. One can see the 
extraordinary variability, encompassing 2 orders of magnitude in flux, both 
in the 7-ray band and in the optical. When the optical flux is in the low 
state, one can also see the contribution of the accretion disk (see the little 
upturn of the optical-UV data at 10^^ Hz). The dashed black line illustrates 
a model for the disk emission (10^^ Hz)and the torus emission (at 10^^ Hz). 

Radio sources in general, and blazars in particular, emit over the entire 
electromagnetic spectrum, from the radio (down to 10^ Hz) to the TeV band 
(up to 10^^ Hz). The main characteristics are: 

• The overall spectral energy distribution (SED), once plotted in uF^, 
shows two broad peaks. The location of the peak frequencies varies 
from object to object, but in general the first peak is between the mm 
and the soft -X-rays, while the high energy peak is in the MeV-GeV 
band. 




Figure 8.12: The overall electromagnetic spectrum of 3C 454.3, the most 
luminous 7-ray source up to now. Note the large amplitude variability, 
even day-by-day. Dates refer to the year 2009. Lines correspond to fitting 
models. See BonnoH et al. (2011, MNRAS, 410, 368). 



• They are variable, at all frequencies, but especially at high energies. 
Minimum variability timescales range between weeks and tens of min- 
utes. 

• In restricted frequency ranges, their spectrum is a power law. 

• The variability is often (even if not always) coordinated and simulta- 
neous in different energy bands (excluding the radio). 

• They are often polarized, in the radio and in the optical. 

• The high energy hump often dominates the power output. 

Besides these common features, it seems that blazars form a sequence 
of SEDs, according to their observed bolometric luminosity. Fig. 18.131 il- 
lustrates this point. It was constructed taking the average luminosity in 
selected bands, and considering ~100 blazars coming from radio and X-ray 
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Figure 8.13: The blazar sequence. As the bolometric luminosity increases, 
the peak of the two humps shifts to smaller frequencies, and the high energy 
hump becomes more important. From Fossati et al. (1998, MNRAS, 299, 
433); Donato et al. (2001 A&A, 375, 739). 

complete samples (complete here means: it contains all blazar whose flux is 
greater than a limiting one). Note that: 

• Low power blazars, that are BL Lacs (i.e. no lines), are bluer than 
powerful blazars, that are FSRQs (with lines). Bluer means that the 
peak frequencies of both peaks are larger. 

• The high energy hump increases its relevance as we increase the bolo- 
metric luminosity. At low luminosities both humps have the same 
power, while the most powerful FSRQ have an high energy hump that 
is ~ lOx the low energy one. 

Interpretation of the SED of blazars 

First hump — It is rather straightforward to interpret the first hump of 
the SED as synchrotron emission. Since it sometimes varies rapidly and 
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simultaneously (at least in optical-UV-soft-X~rays bands) we believe that 
it comes from a single region of the jet. 

High energy hump — It is believed (but with no unanimity) that it is pro- 
duced by the Inverse Compton process, by the same electrons producing 
the synchrotron. The seed photons could be the synchrotron photons them- 
selves (especially in low power blazars) and/or the emission line photons 
(when the lines are present, so only in powerful objects...). Remember that 
in the comoving frame of the jet, the lines are seen blueshifted (factor F) 
and the arrival time is contracted (another factor T) so that in the comoving 
frame the radiation energy density of the lines is seen enhanced by T^. 



Simple estimates 

Suppose to know the luminosity of both peaks of a blazar, call them Lg 
and Lq, and the location of their peak frequency ^ind vq- You also have 
information about the minimum variability timescale tvar- Assume that the 
emitting region is only one, and approximate it with a sphere of radius R. 

Let us consider two cases. For the first we will assume that the high 
energy emission is due to the synchrotron self-Compton process (SSC), in 
the second we will assume that it is due to inverse Compton on the photons 
of the broad line region. For both cases the size R is given by 

R ~ Ctvar (8.21) 

SSC — In this case the ratio between Lq and Ls is given by: 

_ ?^rad _ Ls 1 



Lg U'^ 4ttR^c6^ U'j^ 



1.22) 



where primed quantities refer to the comoving frame. Note the 5 factor 
to pass from the observed to the comoving luminosity. We have assumed 
that C/^g^j is entirely due to the synchrotron luminosity, with no contribution 
from the SSC one. Is that correct? It is, if most of the SSC radiation is at 
frequencies so large that the corresponding scattering process occurs in the 
Klein-Nishina regime, and thus can be neglected. From Eq. 18.221 we can 
derive {U'^ = B'^/Stt): 

B'S' = ^ (^) = {l + z)^ (^) (8.23) 

R \cLcJ Ctvar \cLcJ 

Now let us use the peak frequencies: the ratio between them is 

= Tpeak (8.24) 

Us 
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Figure 8.14: The SED of the BL Lac 0851+202, together with a fitting 
model. Red points are simultaneous, grey points are archival data. 



Note that this is a ratio, therefore there is no 6, no (1 + z). The observed 
synchrotron peak frequency is given by 

- 3-B7peak = 3 ^^7peak (8-25) 

Inserting Eq. 18.241 we arrive to 

B'6='-^^il + z) = (1 + z) (8.26) 

2e i/c 3.7 X 106 z/c 

We have two equations (Eq. 18.231 and Eq. I8.26P for the two unknowns B' 
and 6. 

Example SSC — Fig. KW shows the SED of the BL Lac 0851+201 (=0J 
287), at z = 0.306. For this source we have: Lq ~ Lg ~ 10^^ erg s~^. We 
can guess uq ~ lO^'' Hz extrapolating the X-ray and the 7-ray spectra. Let 
us take z^s ~ 5 x 10^^ Hz, and assume that t^ar = 10^ s. 

From Eq. [823] we have = 3.5 x 10^ G and from Eq. K26\ we derive 
B'6 = 8.8 G. Therefore 6 = 20 and B' = 0.44 G. Notice that if uc were 10^^ 
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Hz instead of 10^*^ Hz, than we would have found 6 = 63 and B' = 0.014 
G. These estimates are very sensitive to the I'c/'^s ratio. The larger the 
separation between the two peaks, the larger 5 and the smaller B'. 

Inverse Compton with broad line photons — In this case let us assume that 
most of the radiation energy density in the comoving frame comes from the 
BLR. We have 

f^r'ad - (8.27) 

We may remember that Lblr ~ O-lL^isk, and that Rblr ~ 10^''ivion,45 cm. 
This assumptions lead to (assuming that Lion ~ Ldisk): 

UL. - ^ (8.28) 

1/2 

a remarkable result, due to the -Rblr oc L^^^y_ dependence. As long as the 
emitting region of the jet is within the BLR, then U^^^ depends only from 
r, but not on the power of the disk or the actual value of the size of the 
BLR. By making the ratio with U'^ we have 

^ _ Kd _ ^ ^ _ /^2Ls y/^ 

The synchrotron peak frequency is still given by Eq. 18.251 if z/s is larger 
than the self absorption frequency The peak frequency of the Compton 
hump is now given by: 

= 3 Tpeak^^Lya (^.30) 

where we have taken the frequency of the Lya line as the characteristic 
frequency of the seed photons (since the Lya line is the most intense). The 
factor T5 comes from considering the scattering in the comoving frame, 
where the typical frequency is 1/^^^, ~ Li/Lya, and then blueshifting the 
frequency after the scattering in the comoving frame by the factor 5. Making 
the ratio vc/^s we have 

UC ^ (4/3)7peak^Ly»'5r/(l+z) ^ Lz^Lyg , . 

(4/3)7^eak^V(l + ^) ^ ^ • ^ 

From this equation we have: 

— ~ 8.8 X 10^ — (8.32) 
1 uc 

Unfortunately, Eq. 18.321 and Eq. 18.291 give two relations for the same quan- 
tity: B' /T. We cannot derive B' and T separately. But we can use them 
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to derive a relation between the peak frequencies and the synchrotron and 
inverse Compton luminosities: 



This can be thought as a consistency check, to see if our assumption of 
inverse Compton on the BLR photons is reasonable. Instead, to the aim to 
derive the parameters, we have to use additional information. One important 
information comes from the superluminal motion of the sources. Usually, it 
is observed at 5-22 GHz. At these frequencies the flux is due to a region of 
the jet self-absorbing there. This region is not the same as the one producing 
most of the emission (we know it because of variability: most of the emission 
is produced in a much more compact region, self-absorbing at ft ~ 100-1000 
GHz). Therefore we must make the assumption that both the viewing angle 
and the bulk Lorentz factor do not change from the two regions. If this is 
true, and assuming that the viewing angle is 9^ ~ l/T, we can set 



and thus solve for B' . Additional information may come from the details of 
the SED, if we see some signs of the SSC emission: even if not dominating 
the bolometric luminosity, it may contribute, especially in the X-rays. 

Example EC — By EC we mean "External Compton", namely that the 
relevant seed photons for scattering are produced externally to the jet, as 
in the case of the BLR photons. Suppose that we observe a source as the 
one shown in Fig. 18.151 Assume /3app = 10. From the figure Lq/Ls = 10. 
Assume uq ~ 10^^ Hz. We do not observe directly the synchrotron peak, 
but the radio spectrum is flat (cx z^'^) and the optical emission is steep [i.e. 
F^v) oc in the optical, once the contribution from the accretion is 

subtracted] . 

From Eq. 18.291 we have B' /T = 0.26. Assuming F = /3app; we derive 
B' = 2.6 G. From Eq. 102] we have us = uq (B'/T) /8.8 x 10^ = 3 x 
10^^ Hz. Is that consistent with what we observe? Yes: the radio and 
optical spectrum indicate a peak (in uF^) between the last observable radio 
frequencies (usually 100 GHz =10^^ Hz) and the optical spectrum. The 
value of z^s = 3 X 10^^ Hz is larger than the self-absorption frequency of our 
emitting region. 

8.9.3 The power of the jet 

It is not easy to estimate the total power of the jet. There are a few ways 
to do it, including studying the energy content of the radio lobes, thought 
to be the calorimeter of the source, and, more recently, the energy required 




(8.33) 



r ~ A 



app 



(8.34) 



8.9. THE JET 



153 




Log u [Hz] 



Figure 8.15: The SED of the blazar 1622-253, together with a fitting model. 
The peak of the Compton emission uq = 10^^ Hz, and the ratio Lq/Ls = 10. 
The figure shows also the contribution from the accretion disk, its X-ray 
corona and the IR torus (short dashed line) and SSC (long dashed). 

for the jet to produce an "X-ray cavity" in cluster of galaxies, around the 
central giant galaxies that are also radio-loud sources. 
The jet power has different components: 

• The jet emits, and the produced radiation must be at the expenses of 
the jet power. Calculating the power (in all directions) in the form of 
radiation sets a lower limit to the jet power. 

• The jet transports matter. If anything, we must have the electrons 
responsible for the radiation we see. They are relativistic in the co- 
moving frame (lower case 7), and have also a relativistic bulk motion 
(capital r). 

• The jet transports also protons. How many? This is controversial, 
because if the emission is totally produced by electron-positrons, there 
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are no protons. If instead the emitting plasma a is a normal electron- 
proton one, than there is one proton per emitting electron. 

• The jet transports a magnetic field. There is a Poynting flux. 

All these components can be formally accounted for by writing 

Pi = ttB?T'^I3cU[ (8.35) 

where U[ is the comoving energy density of the \th component: when dealing 
with the emitting electrons we will have 

U'^ = rriec^ J jN{j)dj = {j)n'^meC^ (8.36) 

When dealing with the magnetic field we have 

Ui, = ^ (8.37) 

and when dealing with the radiation we have 

L' L 

For the protons, assuming them cold. 



%=<—T^ (8-39) 
P me (tX 

The logic is the following: we sit at some place along the jet, cut a cross 
sectional surface of it, and count how many photons, protons, electrons and 
magnetic field are passing by in one second. One factor T accounts for the 
different mass or energy as seen in the lab frame (i.e. a proton will have a 
mass Fmp) and the other factor F comes because the density we see is more 
than the comoving density (lengths are smaller in the direction of motion). 
This is a fiux of energy, i.e. a power. 

It is found that the jet is powerful, often as powerful as the accretion 
disk, and sometimes more. Furthermore, the jet of powerful blazars can- 
not he dominated by the Poynting flux. This is to be expected, since the 
synchrotron component in these sources is only a minor fraction of the bolo- 
metric luminosity, dominated by the inverse Compton component. This 
limits the possible values of the magnetic field. 



8.9.4 Lobes as energy reservoirs 

It is also found that the energy dissipated into radiation is only a minor 
fraction of the jet power, at least in powerful FSRQs. Most of the power 
goes to energize the lobes. There, we have a relaxed region, in the sense 
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that the dimensions are huge, and the magnetic field must be small (i.e. 
B ^ G). Large dimensions also means that the radiation energy density 
is small. Therefore the radiative cooling timescales are very long. Particles 
have not enough time to emit their energy, that instead goes to increase the 
dimension of the region. Extended radio lobes are therefore a reservoir of 
the energy that the jet has provided through the years. If we could measure 
the energy content of the lobe and its lifetime, than we could measure (E/t): 
the average power of the jet. 

Lobes emit in the radio, through synchrotron radiation. They are so 
big that we can easily resolve them. So we know their size. To produce 
synchrotron radiation, there must be some magnetic field and relativistic 
particles. For a given observed synchrotron luminosity, we could have rel- 
atively large magnetic fields and fewer particles, or, on the contrary, less 
magnetic field and more particles. These solutions are not energetically 
equivalent, as we shall see. When the energy content in the magnetic field 
and particles balance, then we have the minimum total energy. To see this, 
consider a lobe of size -Riobe emitting Liobe by synchrotron: 

Aobe = Y^lohe'^eC^ j N {j) js d'J 

= ^i?3^,^laTc|^(72)ne (8.40) 



This gives Ue- 



The energy contained in the magnetic field is 

Eb = Y^obef^ (8.42) 

The energy contained in the particles depends if we have cold or hot protons. 
We do not know. We can parametrize our ignorance saying that the proton 
energy is k times the energy contained in the relativistic electrons. Therefore 
the particle energy in the lobe is 

Air 

Ee,p = Y^k.be(7) "-e^eC^ (1 + /c) 

Note the dependence on B^: Eb oc B^ , while £'e,p oc l/B"^. 

Due to these scalings, the total energy E'tot = Es + Ee^p will have a min- 
imum for Eb = Ee.p. This condition of minimum energy also corresponds 
to equipartition. Fig. 18.161 shows an example, for which the equipartition 
magnetic field is B^^ ~ 10~^ G: with this field we minimize the total energy 
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Figure 8.16: The energy contained in the lobe in magnetic field (E^) and in 
particles (E'e.p)- The black line is the sum. We have assumed that the lobe 
of 100 kpc in size emits synchrotron luminosity -Liobe = 10^^ erg s^^ and 
that protons have the same energy of the emitting electrons (i.e. A; = 1). 
For the particle distribution, we have assumed 7min = 1> 7max = 10^ and 
p = 2.5 [i.e. A^(7) oc 7-2-5]. 



requirement. In the specific case illustrated in Fig. 18.161 the minimum en- 
ergy contained in the lobe is E'tot ~ 10^*^ erg. If we knew its lifetime, we 
could calculate the average jet power. Assuming 10 million years ~ 3 x 10^^ 
s, we have (-Pjet) ~ 3 x 10^^ erg s~^. In this case the lobe synchrotron 
luminosity is only a tiny fraction (i.e. ~ 3 x 10^'^) of the received power. 
This estimates neglects adiabatic losses, likely to be important. Suppose to 
include them: how they modify the estimate on the average jet power? 



8.10 Open issues 

Research in the AGN field is still active, because the open issues are many. 
The previous notes tried to give a glimpse of the "common wisdom" picture 
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of the mainframe about AGNs, but we should avoid to take it as absolutely 
not controversial. The following is an incomplete list of open problems. 

• Black hole masses. We have several ways to estimate the black hole 
masses, but the uncertainties are still large. Knowing the mass more 
precisely, we could know all the quantities in Eddington units, which 
is more physical. 

• Black hole growth. When did the supermassive black hole form? 
It takes time to build them up. Can we live with black hole seeds of 
stellar size or do we need something more exotic? 

• Accretions modes. When the density is small, protons cannot inter- 
act efficiently with electrons (responsible for the emitted radiation and 
the cooling). The energy remain in the protons, that do not cool. The 
pressure increases, and the disk puffs up, and the density becomes 
smaller still. At what accretion rate is there the transition? Are 
Coulomb collisions the only way for protons and electrons to interact? 

• What is the broad line region? We still do not know if it is an 
ensemble of clouds in Keplerian motion, or some form of inflow or 
outflow. 

• The infrared "torus". It should be unstable. Is it formed by a 
clumpy medium? Can it be, instead, a tilted disk? Where does it 
begin: is it a separate structure or just a continuation of the disk? 

• The acceleration mechanisms in the jet. Shocks are the best 
bet, but shocks accelerate protons more than electrons. Accounting 
for hot (and relativistic) protons would increase the estimates for the 
jet power. What is the role of magnetic reconnection? 

• The power of jets. We just started to estimate them in a reliable 
way, but we are still far from robust conclusions. Main uncertainties: 
the amount of electron-positron pairs, and the contribution of hot 
protons. 

• The spin of the black hole. The rotational energy of the black hole 
is a huge energy reservoir that can be extracted. Does the black hole 
spin rapidly? Is it the spin energy that is used to accelerate relativistic 
jets? Is the spin at the root of the radio-loud/radio-quiet dicotomy? 



